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The  development  of  methods  applicable  to  the  X-ray  photoelectron 
spectroscopic  (XPS)  analysis  of  solid  surfaces,  in  general,  and  solid 
state  lead  ion-selective  electrode  (ISE)  membranes,  in  particular,  is 
the  focus  of  the  present  dissertation.  These  methods  fall  into  two 
categories.  The  potential  effects  of  surface  roughness  on  XPS 
analysis  have  been  investigated.   In  addition,  a  method  has  been 
developed  which  allows  the  composition  of  the  near-surface  region  to 
be  investigated  in  a  nondestructive  fashion. 

The  study  of  surface  roughness  effects  has  been  carried  out  from 
both  a  theoretical  and  an  experimental  perspective.  A  program  has 
been  developed  to  determine  the  nature  and  magnitude  of  effects 
attributable  to  surface  topography.  Its  predictions  have  been 
compared  to  experimental  data;  the  agreement  obtained  is  quite  good. 
The  effects  of  the  pressure  used  in  sample  preparation,  which  is 
related  to  the  sample  topography,  on  the  XPS  signal  have  been 


investigated.  Finally,  the  effect  of  increasing  surface  abrasion  (and 
thus  increasing  surface  roughness)  on  the  relative  makeup  of  the 
signal  from  the  various  species  present  in  the  near-surface  region  has 
been  studied. 

In  order  to  obtain  depth  profile  information  nondestructively 
from  a  solid  sample,  a  technique  known  as  variable  angle  XPS  (VAXPS) 
can  be  utilized.  In  this  technique,  the  photoelectrons  are  collected 
at  a  series  of  take-off  angles  with  respect  to  the  mean  plane  of  the 
surface.  These  VAXPS  results  contain  the  desired  depth  profile 
data.  However,  the  data  must  be  mathematically  treated  in  order  to 
obtain  the  depth  profile.  A  method  of  performing  this  data 
transformation  has  been  developed.  Its  description  is  presented, 
along  with  testing  of  the  theory  behind  it. 

Finally,  these  methods,  although  developed  for  general  usage, 
have  been  applied  to  the  study  of  the  surface  composition  of  lead  ISEs 
consisting  of  a  mixture  of  lead  sulfide  and  silver  sulfide.  Several 
studies  have  been  carried  out,  investigating  the  effects  of  membrane 
composition,  preparation,  cleaning,  and  poisoning  of  the  response  on 
the  surface  composition.  The  results  are  presented,  along  with 
several  conclusions  related  to  the  optimization  of  lead  ISE  behavior. 
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CHAPTER  1 
INTRODUCTION 


X-Ray  Photoelectron  Spectroaoopy 
The  roots  of  X-ray  photoelectron  spectroscopy  (XPS)  may  be  traced 
far  back.  One  of  the  four  fundamentally  important  papers  published  by 
Einstein  in  1905  concerned  an  explanation  for  the  photoelectric  effect 
(1).  This  effect  had  long  been  known  to  occur;  however,  Einstein  was 
the  first  person  to  correctly  describe  the  nature  of  the  photoelectric 
effect.  The  phenomenon,  which  is  the  physical  basis  for  XPS,  can  be 
described  simply  as  an  excitation  process,  whereby  a  photon  of  energy 
hv  impinges  on  a  solid,  interacts  with  an  atom's  electron  density 
cloud,  and  causes  the  emission  of  a  so-called  photoelectron,  whose 
kinetic  energy  Ej^  relative  to  the  measuring  spectrometer  is  given  in 
the  absence  of  surface  electric  double  layers  by  the  expression 


Ej^  =  hv  -E^   -^^, 

where  E^^  is  the  binding  energy  of  the  electron  in  the  atomic  or 
molecular  system  measured  relative  to  the  Fermi  level  and  (j)-  is  the 
instrumental  work  function,  all  typically  measured  in  eV  (2).  Thus, 
the  electron's  kinetic  energy  is  independent  of  the  incident  photon 
intensity.  Early  attempts  to  construct  a  high-resolution  electron 
spectrometer  made  use  of  magnetic  focusing  and  thus  were  susceptible 
to  stray  magnetic  fields.  The  theory  behind  such  instrumentation  can 


be  traced  back  to  1956,  when  the  first  double-focusing  magnetic 
spectrometer  was  described  in  the  literature  by  Siegbahn  and  co- 
workers (3).  However,  the  potential  usages  of  such  an  instrument  were 
vastly  increased  with  the  report  by  the  same  group  that  not  only  did 
the  instrument  yield  elemental  information  about  the  surface,  but 
chemical  state  information  as  well.  It  was  discovered  that  Cu  2p 
electrons  from  copper  and  copper  oxide  samples  had  different  kinetic 
energies  and  thus  different  binding  energies.  Combined  with  its 
sensitivity  to  species  at  the  surface  and  in  the  near-surface  region 
(typically,  down  to  approximately  10  nm),  this  factor  has  been 
responsible  for  the  development  of  XPS  as  an  analytical  technique. 

The  next  major  step  in  the  evolution  of  XPS  was  the  elimination 
of  the  problems  caused  by  external  stray  magnetic  fields.  Carlson  and 
Krause  were  the  first  researchers  to  design  an  instrument  using 
electrostatic  focusing.  This  development  was  the  crucial  step  needed 
to  induce  several  instrument  manufacturers  to  construct  and  market 
photoelectron  spectrometers.  The  field  experienced  a  dramatic 
increase  in  popularity  in  the  early  1970s,  as  the  first  generation  of 
commercial  instruments  were  introduced.  Since  then,  such  instrumental 
advances  as  X-ray  monochromators  (to  reduce  the  X-ray  linewidth  and 
thus  enhance  resolution),  multi-channel  detectors  (to  increase  the 
collection  efficiency  and  thus  reduce  run  times),  elaborate  lens 
systems  (to  further  improve  resolution),  various  in  situ  sample 
treatment  methods,  and  small-spot  XPS  instruments  (to  enhance  the 
spatial  resolution)  have  been  introduced.    '   •.- ■  ■ 
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Quantification  in  XPS 
The  achievement  of  quantitative  XPS  results  may  be  approached 
from  two  directions.  One  may  make  use  of  a  number  of  standards  to 
construct  a  series  of  empirical  sensitivity  factors  based  on  a 
comparative  analysis  of  peak  areas.  This  method  was  applied  early  on 
to  collect  a  limited  number  of  sensitivity  factors  (H).     However,  the 
difficulty  in  preparing  suitable  standards  has  limited  the 
applicability  of  this  technique  primarily  to  the  study  of  homologous 
series  of  compounds.  The  other  approach  involves  a  theoretical 
analysis  of  the  chain  of  events  leading  from  the  excitation  of  the 
photoelectron  to  its  energy  determination  and  subsequent  detection.  A 
fundamental  work  by  Fadley  et  al .  (5)  laid  the  foundation  for 
quantification  in  XPS.  This  paper  outlined  for  the  first  time  the 
basic  principles  of  angular  distribution  experiments  in  XPS,  the 
methodology  for  analysis  of  thin  surface  films  and  their  effect  on 
quantitative  analysis  of  the  bulk  sample,  and  the  possible  effects  of 
surface  topography  on  the  XPS  signal.  Since  then,  much  progress  has 
been  made  in  characterizing  the  various  contributing  factors  in  the 
XPS  peak  intensity  equation.  Reports  of  methods  for  quantitative 
analysis  without  standards  began  appearing  around  1977  (6,7).  A  more 
comprehensive  approach  taking  a  larger  number  of  possible  factors  into 
account  was  introduced  in  1982  (8)  and  has  been  revised  since  then 
(9).  The  early  models  typically  made  the  implicit  assumptions  that 
the  sample  composition  was  homogeneous  throughout  the  volume  of  sample 
analyzed,  and  the  surface  topography  had  a  negligible  effect  on  the 
XPS  signal.  Later  modifications  have  made  an  attempt  to  address  these 
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two  factors  in  a  simplistic  fashion.  They  allow  for  the  possibility 
of  an  impurity  overlayer's  interfering  with  the  signal.  The  problem 
of  surface  roughness  in  quantification  has  not  been  addressed  in 
detail;  however,  some  experimental  results  have  been  reported  (10) 
indicating  its  effect  on  the  determination  of  overlayer  thicknesses. 

Measurement  in  XPS 
The  basic  measurement  process  in  XPS  is  outlined  in  Figure  1-1. 
The  technique  makes  use  of  a  photon  source  of  constant  energy, 
typically  the  Ka  X-ray  lines  of  aluminum  or  magnesium,  aimed  at  a 
surface  from  which  photoelectrons  are  emitted.  The  commonly  used 
electron  energy  analyzer  for  XPS,  the  hemispherical  analyzer  (HSA), 
consists  essentially  of  two  closely  spaced  hemispherical  plates, 
across  which  a  potential  can  be  applied,  so  that  only  electrons  of  a 
specified  energy  may  pass  through  to  the  detector.  The  entire 
measurement  process  must  take  place  under  higher  vacuum,  e.g.,  10"^ 
Torr,  in  order  to  slow  the  buildup  of  ubiquitous  hydrocarbon 
contamination  layers,  keeping  the  surface  as  clean  as  possible  and 
allowing  for  greater  sensitivity  toward  the  true  surface  species.  The 
surface  sensitivity  arises  not  from  the  depth  of  penetration  of  the 
incident  X-rays—which  may  travel  100-1000  nm  before  being  appreciably 
attenuated— but  from  inelastic  scattering  of  the  photoelectrons.  This 
process  accounts  for  the  exponential  decay  of  the  photoelectron 
population  as  a  function  of  depth,  with  the  characteristic  inelastic 
mean  free  path  (IMFP)  defined  as  the  straight-line  distance  an 
electron  population  travels  before  its  number  drops  to  1/e  (37$)  of 
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its  original  value  due  to  inelastic  collisions.  For  photoelectrons 
with  typical  kinetic  energies  of  500-1500  eV,  the  IMFP  ranges  from 
0.1-10  nm.      ■■  :  ^  -  -'  -.        :■::' V   ^  ' 


Instrumentation 


The  instrument  used  to  collect  the  majority  of  the  XPS  data  is  a 
Kratos  XSAM  800  X-ray  photoelectron  spectrometer.  A  block  diagram  of 
the  spectrometer  is  shown  in  Figure  1-2.   This  unit  consists  of  an 
ultra-high  vacuum  (UHV)  spectrometer  console,  a  sample  introduction 
system,  an  electronics  console,  and  a  computer  for  instrument  control 
and  data  manipulation.  The  spectrometer  console  comprises  a  spherical 
sample  analysis  chamber,  with  a  water  cooled  dual  anode  X-ray  gun, 
electron  flood  gun,  macrobeam  ion  gun,  electron  energy  analyzer, 
detector  system,  and  a  chamber  pumping  system.  The  X-ray  gun 
furnishes  either  Al  Ka  (1486.6  eV)  or  Mg  Ka  (1253.6  eV)  radiation, 
allowing  differentiation  between  kinetic  energy- invariant  processes, 
such  as  Auger  electron  spectroscopy  (AES),  and  binding  energy- 
invariant  processes,  such  as  XPS.  For  reasons  mentioned  below,  the 
bulk  of  the  experimental  work  was  carried  out  using  the  Al  anode.  The 
X-ray  gun  may  be  isolated  from  the  analysis  chamber  and  separately 
pumped  with  a  small  dedicated  ion  pump,  to  reduce  main  chamber  vacuum 
degradation.  The  electron  flood  gun  is  used  to  neutralize  the  charge 
buildup  on  the  surface  of  insulating  samples  due  to  electron 
emission.  The  ion  gun  is  differentially  pumped  using  a  Leybold- 
Heraeus  TMP  150  turboraolecular  pump,  which  has  a  nominal  pumping  speed 
of  150  L/s.  The  ion  gun  is  typically  used  to  sputter  clean  surfaces 
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for  analysis  using  kr*   ions,  thus  reducing  the  intensity  of  the 
ubiquitous  hydrocarbon  signal  to  a  manageable  level;  it  may  also  be 
used  as  an  excitation  source  in  performing  ion  scattering  spectroscopy 
(ISS),  or  as  an  etching  source  when  carrying  out  destructive  depth 
profiling  of  a  sample  surface.  The  XSAM  800  uses  a  180°  deflection 
double  focusing  hemispherical  electrostatic  analyzer  with  a  mean 
radius  of  127  mm  and  an  aberration  compensated  input  lens  (ACIL) 
system.  A  choice  of  two  analyzer  modes  is  given.  In  the  fixed 
retarding  ratio  (FRR)  mode,  the  kinetic  energy  of  each  electron 
entering  the  ACIL  is  reduced  by  a  fixed  fraction  of  its  initial 
kinetic  energy  as  it  moves  through  the  HSA.  In  the  fixed  analyzer 
transmission  mode  (FAT),  the  electron's  kinetic  energy  is  reduced  to  a 
constant  value  for  all  electrons  passed  to  the  detector.  The  relative 
kinetic  energy  information  necessary  to  distinguish  the  electrons 
energetically  is  thus  contained  in  the  ACIL  focusing  program  to  reduce 
electron  energy.  The  FRR  mode  gives  higher  count  rates  for  peaks  with 
low  binding  energies  (high  kinetic  energies);  this  region  being 
predominantly  the  region  of  study,  the  FRR  mode  was  used  in  nearly  all 
experiments.  The  detector  system  consists  of  a  fast  single  channel 
multiplier  and  a  fast  response  head  amplifier.  It  can  be  operated  in 
either  a  pulse  counting  mode,  which  gives  a  true  electron  count  rate, 
or  a  current  mode,  which  acts  as  a  pulse  count  integrator.  The  former 
mode  is  more  sensitive;  however,  it  leaves  the  electron  multiplier 
susceptible  to  saturation.  For  this  reason,  the  current  mode  was  used 
universally  in  these  studies.  The  chamber  pumping  system  comprises  a 
rotary  vane  roughing  pump,  the  dedicated  X-ray  ion  pump,  the 
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turboraolecular  pump  (which  can  also  be  valved  to  pump  on  the  main 
chamber),  a  titanium  sublimation  pump,  a  200  L/s  ion  pump,  associated 
pumping  manifolds  and  valves,  and  several  vacuum  measuring  devices. 
Pressure  in  the  roughing  vacuum  lines  is  monitored  by  two  thermocouple 
gauges.  An  ionization  gauge  is  attached  to  the  sample  chamber;  it 
serves  to  monitor  pressure  predominantly  in  the  10~^  -  io~^  Torr 
range.  For  UHV  pressure  measurements,  a  cold  cathode  Penning  trigger 
gauge  is  provided,  with  an  operating  range  of  10   -  10~^^  Torr.  In 
order  to  prevent  inadvertent  instrument  operation,  an  interlock  system 
has  been  provided.  If  the  ionization  gauge  cannot  be  energized,  i.e., 
the  pressure  is  >10~^  Torr,  the  analysis  and  detection  electronics 
cannot  be  turned  on.  Also,  if  the  X-ray  gun  ion  pump  is  not  on,  or  if 
the  anode  cooling  water  flow  is  not  sufficient,  the  X-ray  source 
cannot  be  operated.  The  system  has  a  base  operating  pressure  of 
5  X  10    Torr;  typical  operating  pressures  are  in  the  low  10~  to 
high  10    Torr  range.   -'"  .  : 

The  sample  introduction  system  allows  a  sample  to  be  entered  into 
the  UHV  chamber  from  atmospheric  pressure  in  a  matter  of  minutes.   It 
comprises  a  direct  insertion  lock  with  an  isolation  ball  valve  and  a 
fast  insertion  sample  probe.  The  insertion  lock  is  pumped  by  the 
rotary  vane  roughing  pump.  The  sample  probe  can  accommodate  a  variety 
of  probe  tips  for  analyzing  different  types  of  samples.   It  can  also 
be  heated  or  cooled  over  a  temperature  range  of  -ISQOC  to  +600 °C.  The 
probe  can  be  rotated  on  its  axis,  thus  altering  the  orientation  of  the 
sample  surface  relative  to  the  X-ray  gun  and  ACIL  entrance  slit. 
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The  electronics  console  contains  the  control  units  for  pressure 
monitoring  and  sample  temperature  setting,  X-ray  power  supply,  X-ray 
filament  emission  current  stabilizer,  ion  gun  control  and  raster 
units,  electron  flood  gun  control  unit,  and  control  units  for  both  ion 
pumps  and  the  titanium  sublimation  pump.  It  also  contains  the  scan 
control  unit  and  a  display  control  unit,  both  of  whose  functions  have 
been  largely  usurped  by  the  computer  control  and  display  system. 

The  instrument  is  controlled  by  a  Digital  Equipment  Corporation 

Micro  PDP-11  computer.  A  Tektronix  4105A  color  computer  display 

terminal  serves  as  a  user  interface,  and  a  Tektronix  4695  seven-color 

graphics  copier  provides  hardcopy  output.  The  computer  has  256  Kbytes 

"J  f 
of  random  access  memory,  a  dual  5  1/4"  floppy  diskette  drive  system, 

and  a  Winchester  hard  disk  system  with  a  formatted  capacity  of  11.059 
Mbytes.  The  DS800  data  system  included  with  the  system  resides  on  the 
hard  disk.  Backup  files  can  be  stored  on  floppy  diskettes,  each 
having  a  capacity  of  400  kbytes.  The  data  system  allows  computer 
control  of  data  acquisition,  including  spectrometer  control  and 
extensive  data  display  and  manipulation  features,  such  as  peak  area 
determination,  background  subtraction,  peak  deconvolution,  and  curve 
fitting. 

■         Lead  Ion-Selective  Electrodes 
Electrodes  selective  to  lead  ions  date  from  1969  (11).  The  first 
such  electrodes  utilized  a  homogeneous  membrane  consisting  of  a  co- 
precipitate  of  lead  sulfide  (PbS)  and  silver  sulfide  (Ag2S)  pressed 
into  a  pellet,  the  potential  across  which  can  be  related  to  the  lead 
concentration.  Other  lead  ion-selective  electrodes  (ISEs)  were  soon 
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introduced,  including  a  liquid  membrane  ISE  (12);  a  homogeneous 
membrane  using  lead  selenide  (PbSe)  or  lead  telluride  (PbTe)  in  place 
of  PbS  (13);  another  homogeneous  system  consisting  of  PbS,  AgpS,  and 
CU2S  (14);  and  several  homogeneous  membranes  with  PbS  or  PbS/AgjS 
present  at  the  electrode  surface  (15,16).  Since  then,  other  lead  ISEs 
have  been  proposed,  utilizing  PbS/AggS  on  a  PTFE/graphite  support 
(Selectrode),  PbS  deposited  on  an  ionic  conducting  support  material 
(17),  PbOg  on  a  Selectrode  body  (18),  coated  wire  type  Pb  ISEs  (19), 
and  PbS/AgpS/ASpS,  glass  membranes  (20). 

However,  by  far  the  most  popular  designs  have  been  the 
homogeneous  PbS/AgjS  membrane  electrode  and,  to  a  lesser  extent,  the 
Selectrode  body  with  PbS/AgjS  impregnated  on  the  surface.  For  this 
reason,  we  have  concentrated  our  investigations  on  this  powder  mixture 
pressed  into  a  homogeneous  pellet.  The  number  of  studies  carried  out 
on  such  membranes  is  too  great  to  fully  document.  Several  detailed 
electrochemical  studies  have  been  carried  out  (21-23).  The  electrodes 
typically  show  a  linear  response  in  a  potential  versus  log  apj^  plot, 
with  a  response  sensitivity  of  27-30  mV/decade ,  close  to  the  Nernstian 
value  for  Pb   of  29.5  mV/decade.  Several  common  ions  interfere  with 
the  response  of  the  ISE,  including  Ag*,  Cu^"*",  and  Hg^"",  all  of  which 
cause  surface  species  deposition.  The  ions  Mno"  and  Cr  0^~  act  as 
oxidizing  agents,  thus  interfering  irreversibly  with  the  electrode 
response.  The  response  tends  to  degrade  over  a  period  of  several 
weeks'  usage;  however,  the  electrode  can  be  rejuvenated  by  polishing 
its  surface,  or  exposing  it  to  either  ammoniacal  EDTA  or  HClOi, 
solution.  ■  V'  ■■  ■:  V 
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Relatively  few  studies  have  been  carried  out  to  investigate  the 
surface  behavior  of  Pb  ISEs.  Chaudhari  et  al .  (2M)  were  the  first 
investigators  to  use  surface  sensitive  techniques  to  analyze  this 
system.  They  investigated  the  effects  of  various  cleaning  agents  on 
the  surface  composition.  Fixed  angle  XPS  was  also  used  to  examine  the 
effects  of  Cu  *  and  Zn^*  on  these  membranes  (25).  Scanning  electron 
microscopy  (SEM)  with  an  energy-dispersive  X-ray  spectrometer  has  also 
been  applied  to  the  investigation  of  these  ISEs  (26).  A  more 
sophisticated  XPS  study  involving  a  rudimentary  approach  to  resolving 
the  species  concentrations  as  a  function  of  depth  below  the  surface 
has  been  carried  out  (27). 

Thus,  precedent  exists  for  using  surface  techniques  for 
investigating  Pb  ISE  membranes.  The  choice  of  XPS  as  the  analytical 
tool  of  choice  can  be  justified  by  the  nature  of  such  electrode 
membranes.  Their  response  is  governed  chiefly  by  the  interaction  of 
the  solution  being  monitored  with  the  electrode  surface.  An 
understanding  of  the  surface  chemistry  is  critical  to  elucidating  the 
mechanism  of  Pb  ISE  behavior,  particularly  with  regard  to  the 
destruction  of  this  behavior  by  various  agents  present  in  solution. 
As  opposed  to  other  surface  techniques,  XPS  stands  out  as  the  method 
of  choice  because  of  its  ability  to  resolve  chemical  state  information 
of  surface  species.  The  PbS/Ag2S  system  is  a  complex  one,  with 
several  different  molecular  species  present  at  the  surface.  A  mere 
knowledge  of  the  relative  amounts  of  Pb,  Ag,  and  S  is  not  sufficient 
to  provide  a  description  of  the  method  of  membrane  operation;  one  must 
know  the  nature  of  the  species  which  interact  with  the  electrolyte 
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solution.  However,  the  concentration  gradients  of  the  important 
species  in  the  near-surface  region  have  not  been  investigated  in 
detail.  In  addition,  the  possibility  that  the  sample's  surface 
topography  may  affect  the  XPS  results  must  be  considered  in  order  for 
the  results  to  be  treated  with  confidence.  These  factors  have  led  to 
the  present  study.  A  model  for  determining  the  possible  effects  of 
surface  topography  on  the  nature  and  magnitude  of  the  XPS  signal  has 
been  developed  and  applied  to  the  study  of  Pb  ISEs.  In  addition, 
several  more  general  ramifications  of  the  surface  roughness  effect 
have  been  elucidated  by  a  combination  of  theoretical  and  experimental 
investigation.  These  results  are  reported  in  the  following  chapters. 
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■  "  .:   CHAPTER  2 
THEORETICAL  BACKGROUND 


In  order  to  present  an  equation  for  the  XPS  intensity,  certain 
assumptions  must  be  made.  We  initially  assume  a  surface  in  which 
elastic  scattering  and  electron  transport  anisotropies  are 
neglected.  Elastic  scattering  refers  to  the  process  where  an 
electron's  directional  vector  is  altered  with  a  conservation  of 
energy.  Thus,  the  electron,  if  detected,  will  appear  in  the  species 
peak  of  interest.  More  will  be  said  about  this  factor  below.  The 
condition  of  negligible  electron  transport  anisotropies  implies  that 
no  preferences  are  observable  for  electron  motion  along  a  particular 
direction.  This  condition  is  generally  met  in  practice  when  dealing 
with  polycrystalline  or  amorphous  solids;  however,  single  crystals 
show  evidence  of  anisotropies  which  can  be  used  to  obtain  information 
about  surface  bond  orientation  (28).  Also,  X-ray  reflection  and 
refraction  are  assumed  to  be  negligible.  This  condition  is  met  for 
all  geometries  except  grazing  photoelectron  exit,  which  is  not  of 
interest  in  the  present  study.  Finally,  the  X-rays  and  electrons  are 
assumed  to  be  attenuated  exponentially,  the  decay  being  characterized 
by  a  mean  free  path  (MFP).  For  the  electrons,  the  characteristic 
decay  length  is  the  inelastic  mean  free  path  (IMFP  or  A).  The  X-ray 
MFP  is  typically  much  larger  than  the  electron  IMFP  for  the  electron 
kinetic  energies  of  interest  in  XPS.  Given  these  assumptions,  the 
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intensity  I  (E  ,y)  of  photoelectrons  from  level  y  of  species  c  in 
sample  s  is  given  as  (29)  .,      ''     ' 

I  (E  ,y)  =  I  (hv)  daj[/dn(hv,(|))  F(E^,E.)  T(E^,E.,n  ) 

CC         U         w  CACAO 

f 

•  D(E  ,EJ  G(n  ,9)  f"  N  (z)  exp(-z/X  (E  )'Sin  9)  dz. 
cA    0''oc  so 

Here,  I  (hv)  corresponds  to  the  incident  X-ray  flux  of  energy  hv; 

y 

da  /dn(hv,(j))  is  the  differential  photoionization  cross-section  of  sub- 
shell  y  of  species  c,  also  dependent  on  X-ray  energy  and  the  angle  (^ 
between  the  incident  X-rays  and  the  emergent  photoelectrons  which  can 
reach  the  analyzer.     It  is  given  by   (30)    *  •  '  '^^ 

rv  ?  V  "  -V-'  -  -:.  . 

da^/dn  =  (a^ATr)Cl  -  B^,(E^)/4  (3  cos^<«)  -  1)], 

where  a  is  the  total  cross-section  of  sub-shell  y  of  species  c, 

typically  measured  at  1487  eV  in  units  of  the  C  Is  cross-section  of 

13600  barns,  and  3  . (E  )  is  the  angular  asymmetry  factor,  a  function 

of  the  electron  energy  E„.  The  values  of  a^  have  been  theoretically 

*-  c 

determined  and  reported  as  a  function  of  atomic  number  and  core  level 

(31),  and  they  agree  well  with  experimental  results  (32-34).  Although 

the  photoelectron  energy  does  not  depend  on  the  angle  between  the 

incident  X-ray  excitation  beam  and  the  detected  photoelectron  stream, 

the  intensity  does.  The  angular  asymmetry  terra  takes  this  variation 

into  account.  On  the  Kratos  XSAM  8OO,  this  angle  is  fixed  at  80°; 

thus,  this  factor  becomes  a  constant  for  a  given  peak.  Theoretical 

values  for  3  ,(E  )  have  been  reliably  calculated  for  all  core  levels 
nl  c 

and  tabulated  as  a  function  of  atomic  number  (35),  and  the  effects  of 
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elastic  scattering  on  the  angular  distribution  of  photoelectrons  have 
been  studied  as  well  (36,37).  The  electron-optical  factor  F(E  .E  ) 
allows  for  the  effects  of  deceleration  of  the  photoelectrons  from 
energy  E^  to  analyzer  energy  E^.   In  order  to  explain  why  the 
electrons  must  be  decelerated  before  being  energy  analyzed,  it  is 
necessary  to  examine  the  resolving  power  of  the  photoelectron  energy 
analyzer.  For  a  point  source  imaged  into  a  point  image,  the  resolving 
power  is  given  by  -;*'  '"  '^  '""  "•  "^  '" -^  '"• 

A  V ,  I  » 1  V  H 
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Resolving  power  =  AE/E  =  (W/2R)  +  Ko  , 

where  W  is  the  slit  width,  R  is  the  mean  analyzer  radius,  a  is  the 
serai-angle  in  the  dispersive  plane,  and  K  is  a  constant.  For  a  slit 
width  of  8  mm  and  a  mean  analyzer  radius  of  127  mm,  the  resolving 
power  is  >0.012.  Thus,  for  an  electron  with  1000  eV  kinetic  energy, 
resolution  is  about  12  eV,  far  too  poor  to  distinguish  between 
chemical  states  from  a  particular  core  level.  By  removing  energy  from 
the  electron  before  it  is  energy  analyzed,  one  can  improve  the 
resolving  power,  at  the  expense  of  a  loss  in  sensitivity.  Two  methods 
of  energy  removal  can  be  applied.  In  fixed  analyzer  transmission 
(FAT)  mode,  the  transmission  energy  of  the  analyzer  is  held  constant, 
and  the  electrons  are  retarded  by  a  varying  proportion.  In  fixed 
retarding  ratio  (FRR)  mode,  electrons  entering  the  analyzer  have  their 
energy  retarded  by  a  fixed  ratio  of  their  initial  kinetic  energy.  The 
analyzer  transmission  function  T(E^,E^,n^)  is  dependent  on  electron 
kinetic  energy,  the  analyzer  energy,  and  the  solid  angle  of 
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acceptance  Q^  for  the  analyzer;  D(E  ,E,)  is  the  detector  efficiency; 
o  c  A 

and  G(n  ,0)  is  a  geometric  factor  allowing  for  the  change  in  effective 
area  of  sample  analyzed.  For  the  Kratos  XSAM  800  instrument,  the 
geometric  factor  is  not  a  function  of  9  when  operating  in  the  FRR 
mode.  N  (z)  refers  to  the  atom  concentration  (atoms/cm^)  of  species  c 
at  depth  z;  A  (E  )  is  the  electron  IMFP,  dependent  on  the  sample  and 
electron  energy;  and  9  is  the  photoelectron  escape  angle.  The 
product  X    (E  )«sin  9  denotes  the  effective  IMFP,  the  effective  depth 
from  which  photoelectrons  emerge  (Figure  2-1).  For  a  series  of  peaks 
collected  at  various  values  of  9,  the  average  distsince  an   electron 
travels  before  being  inelastically  scattered  remains  constant,  even 
though  the  effective  depth  below  the  surface  from  which  the  electron 
emerged  varies.  The  lower  the  angle  of  photoelectron  exit  from  the 
surface  with  respect  to  the  tangent,  or  slope,  at  the  point  of  exit, 
the  smaller  will  be  the  depth  from  which  the  electron  emerged.  As 
will  be  seen  in  later  chapters,  this  effect  has  a  profound  impact  on 
the  nature  of  XPS  depth  profiling  experiments  and  the  determination  of 
quantitative  information  from  them.  '-:. 

The  actual  determination  of  the  IMFP  is  difficult  to  achieve  from 
either  an  experimental  or  a  theoretical  approach.  Experimentally, 
several  researchers  have  published  universal  mean-free  path  curves  as 
a  function  of  electron  kinetic  energy  (38-41).  A  power  law  dependence 
with  an  exponent  in  the  range  of  0.5-0.75  is  typically  reported, 
although  more  sophisticated  methods  of  determination  exist  (41).  From 
a  theoretical  viewpoint,  several  workers  (42-44)  have  modeled  the  mean 
free  path  of  a  free-electron-like  material  and  shown  agreement  between 
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their  results  and  experimental  measurements.  Neither  approach  yields 

results  of  sufficient  accuracy  to  preclude  further  research  in  this 

area,  which  is  ongoing.  The  uncertainty  in  the  IMFP  determination  is 

a  chief  complicating  factor  in  performing  absolute  quantitative  XPS 

completely  from  basic  principles.  At  least  part  of  the  uncertainty 

arises  from  the  unaccounted-for  contribution  of  the  elastic  mean  free 

path  to  the  overall  mean  free  path,  as  well  as  its  effect  on  the  XPS 

signal  intensity;  the  latter  effect  is  related  to  the  variation  in  the 

angular  asymmetry  term  discussed  above.  The  general  consensus  among 

researchers  is  that  there  is  a  significant  elastic  scattering  effect, 

but  it  does  not  in  general  affect  the  exponential  decrease  of  the 

signal.  Rather,  it  acts  to  effectively  increase  the  electron  mean 

free  path  (45,46).  The  elastic  scattering  contribution  can  be  reduced 

by  using  small  solid  angles  of  acceptance  at  the  energy  analyzer;  it 

1 
is  also  more  pronounced  at  smaller  values  of  the  photoelectron  escape 

angle  (47). 

The  equation  shown  above  indicates  that  the  intensity  depends  on 

many  parameters  and  points  out  what  are  likely  to  be  some  major 

limitations  of  quantitative  XPS.  First  of  all,  a  reliable  method  of 

estimating  the  asymmetry  parameter,  photoemission  cross-section,  and 

IMFP  must  be  found,  in  order  to  account  for  these  factors'  effects  on 

the  XPS  intensity.  As  discussed  above,  extensive  work  has  been 

directed  toward  this  end.  The  first  two  factors  are  generally 

recognized  to  be  well  characterized;  however,  research  is  still 

ongoing  into  the  nature  of  the  IMFP.   Provision  has  been  made  for  the 

uncertainty  in  this  parameter  in  the  development  of  the  methods  to  be 
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discussed  in  this  work.  Secondly,  surface  roughness  may  have  a 

significant  effect  on  0  and  on  G(Q  ,9).  Quantitative  VAXPS  will  be 

o 

generally  useful  only  when  the  effect  of  surface  roughness  may  be 
adequately  taken  into  consideration.  One  of  the  purposes  of  this  work 
is  to  clarify  the  surface  topography  effect,  both  from  a  theoretical 
and  an  experimental  standpoint.  The  various  instrument-dependent 
functions  remain  constant  in  a  VAXPS  experiment,  where  only  the  escape 
angle  is  changed.  Thus,  keeping  in  mind  the  possible  surface 
topography  effects,  as  well  as  the  uncertainty  in  the  IMF?,  these 
factors  can  be  collected  for  a  VAXPS  experiment  into  an  equation  of 
the  form 

■  :  r-n      *  t 

I  (9)  =  K  (0)  r  N^(z)  exp  [-Z/X  (E  )-3in  9]  dz. 

O        C     ^  0   C  3   C 

The  angular  dependence  of  K  (9),  if  present,  can  typically  be 

determined  empirically.  Results  from  VAXPS  experiments  may  be  used  to 

characterize  the  nature  of  the  intensity  I  (9)  as  a  function  of  the 

c 

photoelectron  take-off  angle  9.  The  use  of  such  VAXPS  data  in  this 
equation  provides  an  opportunity  for  obtaining  information  about  the 
depth  profile  N^Cz)  of  the  species  under  investigation.  The 
methodology  necessary  for  extracting  depth  profiles  from  VAXPS  data  is 
presented  in  a  later  chapter,  as  are  several  applications. 


CHAPTER  3 
MODELING  OF  SURFACE  TOPOGRAPHY  EFFECTS  IN  XPS 


Introduction 

In  order  to  achieve  the  desirable  goal  of  quantitative  analysis 
of  surface  species  by  electron  microscopy,  one  important  prerequisite 
is  the  characterization  of  the  surface  topography  of  the  sample  under 
study.  Surface  roughness  may,  depending  on  its  magnitude,  affect  the 
relationship  between  the  instrumental  setting  of  the  angle  of 
photoelectron  escape  and  the  nature  and  magnitude  of  the  resultant  XPS 
signal. 

Previous  studies  of  the  effects  of  surface  topography  on  the 
magnitude  and  nature  of  photoelectron  signal  almost  unanimously 
involved  attempts  to  model  the  surface  roughness  as  a  periodic 
function.  Ebel  et  al .  (48)  made  use  of  a  model  based  on  close-packed 
cubes  to  approximate  the  effects  of  surface  roughness  on  the 
photoelectron  signal  magnitude  of  gold  films  on  abraded  surfaces. 
Their  rudimentary  model  took  account  of  photoelectron  shading 
(discussed  below),  but  did  not  examine  the  changes  in  sampling  depth 
occasioned  by  the  topography  of  a  sample  surface.  A  more 
comprehensive  study  by  Fadley  et  al .  (5)  made  use  of  a  one-dimensional 
sinusoidal  model  to  predict  changes  in  the  average  emission  angle, 
photoelectron  intensity  reduction  effects,  and  the  alteration  of  the 
average  photoelectron  escape  depth.  Although  several  samples  were 
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analyzed,  no  systematic  experimental  investigation  of  surface 
roughness  effects  was  undertaken.  The  change  in  emission  angle  was 
focused  upon.   Baird  et  al .  (49,50)  experimentally  verified  the 
presence  of  surface  roughness  effects  in  the  XPS  analysis  of 
triangularly  periodic  aluminum  diffraction  gratings  and 
unidirectionally  polished  aluminum  foil  samples.  On  the  basis  of 
their  findings,  they  stated  that  the  effects  of  surface  topography 
would  be  nonnegligible  if  the  roughness  was  at  least  of  the  same 
magnitude  as  the  inelastic  mean  free  path  \. 

More  recently,  Wagner  and  BrOmmer  (51 )  investigated  four  one- 
dimensional  models  of  surface  roughness,  the  two  models  discussed 
above  as  well  as  a  close-packed  semi -circular  model  and  a  rectangular 
roughness  model.  No  experimental  comparison  was  made.  Wu  et  al.  (52) 
carried  out  a  purely  experimental  analysis  of  surface  roughness 
effects  on  gold  films,  concentrating  on  the  reduction  in  photoelectron 
intensity  with  increasing  surface  area,  as  measured  by  the  Brunauer, 
Emmett,  and  Teller  (BET)  technique.  They  made  no  attempt  to 
investigate  the  variation  of  the  XPS  signal  as  a  function  of  escape 
angle.  De  Bernardez  et  al .  (53)  constructed  curves  showing  the 
effects  of  surface  roughness  on  the  number  of  photoelectrons  as  the 
escape  angle  was  varied.  These  curves  were  constructed  using  a  two- 
dimensional  conical  surface  roughness  model,  whose  photoelectron 
emission  was  computed  using  a  Monte  Carlo  electron  trajectory 
simulation.  They  compared  their  theoretical  results  to  the 
experimental  work  of  Wu  et  al .  (52).  The  effect  of  electron  shading 
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on  the  determination  of  thin  overlayer  thicknesses  has  also  received 
some  attention  (10,54). 

Surface  roughness  effects  are  important  in  Auger  electron 
spectroscopy  (AES)  as  well.  The  importance  of  induced  surface 
roughness  due  to  the  use  of  Ar"^  ion  sputtering  was  recognized  early  on 
(55,56)  and  has  generated  extensive  interest.  Surface  roughness 
affects  Ar  ion  sputtering  by  decreasing  depth  resolution  (57), 
reducing  ion  sputtering  yield  (58),  and  changing  the  evolution  of 
topography  during  sputtering  (59).  Holloway  (60)  investigated  the 
surface  roughness  effect  on  the  AES  signal,  concluding  that  its 
presence  tends  to  reduce  signal  magnitudes  and  affect  quantification, 
regardless  of  whether  absolute  or  relative  intensities  are  used. 

In  a  recent  paper  (27),  a  sinusoidal  model  of  surface  roughness 
using  data  derived  from  profilometer  tracings  showed  a  negligible 
roughness  effect  for  the  particular  samples  studied  at  the 
experimental  resolution.  However,  the  results  suggested  that  a  more 
detailed  analysis  of  the  effects  of  surface  topography  on  the  XPS 
signal  was  necessary.  It  has  been  shown  that  the  depth  resolution  is 
lowered  by  the  presence  of  surface  roughness  (50).  In  addition,- other 
effects  related  to  the  change  in  the  angle  of  photoelectron  escape  are 
present,  and  they  affect  the  makeup  of  the  XPS  signal.   Also,  the 
presence  of  photoelectron  shading  acts  to  reduce  the  overall  area 
sampled  by  XPS;  one  must  account  for  its  effects.  In  this  chapter, 
the  theory  and  methodology  of  a  technique  developed  for  the  analysis 
of  surface  roughness  is  presented.  Applications  of  the  technique  for 
the  analysis  of  variable  angle  X-ray  photoelectron  spectral  data,  as 
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well  as  a  comparison  to  such  experimental  data  and  a  discussion  of  the 
results,  are  presented  in  the  next  chapter. 

Theory  of  Modeling 
In  order  to  present  an  equation  for  the  XPS  intensity,  certain 
assumptions  must  be  made.  We  initially  assume  a  surface  in  which 
elastic  scattering  and  electron  transport  anisotropies  are 
neglected.  Also,  X-ray  reflection  and  refraction  are  assumed  to  be 
negligible.  Finally,  the  X-rays  and  electrons  are  assumed  to  be 
attenuated  exponentially,  the  decay  being  characterized  by  a  mean  free 
path  (MFP).  For  the  electrons,  the  characteristic  decay  length  is  the 
inelastic  mean  free  path  (IMFP  or  X)  described  above.  The  X-ray  MFP 
is  typically  to  be  much  larger  than  the  electron  IMFP.  In  this  case, 
the  intensity  I  (E  ,y)  of  photoelectrons  from  level  y  of  species  c  in 
sample  s  has  been  shown  to  follow  the  form 

I^(E^,y)  =  I^Chv)  da^/dQ(hv,<|))  F(E^,,E^)  T(E^,E^,Q^) 

•  D(E  ,EJ  G(fl  ,9)  f"  N  (z)  exp(-z/X  (E  )-sin  9)  dz. 

CA      O-'OC  SC 

Of  interest  in  this  chapter  are  the  geometric  factor  G(Q  ,9)  and 
the  integral  term.  A  surface  with  appreciable  topography  possesses  a 
surface  area  greater  than  the  apparent  area  defined  by  the  projection 
of  the  aperture  onto  the  surface,  due  to  the  presence  of  peaks  and 
valleys  rather  than  a  planar  surface.  However,  only  a  fraction  of 
this  rough  surface  is  typically  capable  of  emitting  photoelectrons 
which  may  reach  the  analyzer  aperture,  due  to  electron  shading  of 
portions  of  the  surface  (i.e.,  valleys)  by  other  areas  of  the  surface 
(i.e.,  peaks).  This  electron  shading  is  an  important  factor  to 
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consider  in  quantitative  analysis,  since  its  magnitude  is  related  to 
the  actual  amount  of  surface  area  "visible"  to  the  analyzer.  The 
method  of  determination  of  the  shading  fraction  assumes  that  the 
electrons  act  as  particles  traveling  in  straight  lines  after  exiting 
the  surface.  A  geometrical  approach  to  the  shading  phenomenon  is 
utilized,  wherein  electrons  free  of  the  surface  and  travelling  toward 
the  analyzer  entrance  slits  are  inelastically  scattered  if  their  path 
causes  them  to  re-enter  the  surface.  Within  the  integral  term,  the 
angle  of  photoelectron  escape  is  affected  by  surface  roughness  as  well 
(Figure  3-1).  Note  that  the  photoelectron  escape  angle  is  defined  for 
a  point  on  a  surface  relative  to  a  tangent  drawn  at  that  point.  Thus, 
at  any  point  on  a  rough  surface,  the  actual  photoelectron  escape  angle 
differs  from  the  instrumental  setting—which  is  based  on  the 
assumption  of  an  atomically  flat  surface— by  a  factor  related  to  the 
tangent,  or  slope,  of  the  surface  at  that  point.  As  mentioned  above, 
the  actual  value  for  the  sine  of  the  angle  of  photoelectron  escape  is 
directly  related  to  the  effective  sampling  depth.  If  the  species  c 
has  a  nonconstant  depth  profile  N^Cz),  then  the  surface  roughness  will 
affect  the  XPS  signal  by  changing  the  depth  from  which  the 
characteristic  photoelectrons  are  emerging.  Since  the  concentration 
of  atoms  may  vary  with  depth,  this  change  can  affect  the  total  number 
of  photoelectrons  generated.  To  summarize,  among  the  effects  which 
surface  roughness  may  exert  on  the  XPS  signal  intensity,  the  change  in 
surface  area,  the  presence  of  electron  shading,  and  the  change  in  the 
actual  range  of  photoelectron  escape  depths  are  prominent. 
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•  ';•"  :   .  ■-.  Methodology   -  ;  .'  .  ..'^  ■ 

A  program  to  estimate  the  magnitude  of  these  factors  has  been 
developed.  A  flow  chart  Is  shown  in  Figure  3-2.  As  input,  it  relies 
on  a  digitized  vector  or,  in  the  more  general  case,  a  matrix  of  values 
representing  the  height  above  or  below  some  defined  surface  level  as  a 
function  of  displacement  along  the  surface.  From  these  data,  a 
shading  fraction  Fg,  which  represents  the  amount  of  surface  free  from 
shading,  is  calculated  by  comparing  the  height  of  a  point  on  the 
surface  with  the  heights  of  the  other  successive  points  along  the 
direction  of  photoelectron  emission,  as  shown  in  Figure  3-3. 
Specifically,  if      ■  ^^^  -  *  '•  • 


h,  +nb<h.      n=1,2,..., 
i        i+n 

then  that  particular  point  i  will  be  shaded  from  emitting  photo- 
electrons  capable  of  reaching  the  analyzer  aperture.  Thus,  the 
program  works  under  the  assumption  that  adjacent,  shading  areas  of  the 
surface  do  so  completely,  with  no  probability  of  electron  trans- 
mission. Here,  h^  represents  the  height  of  point  i,  b  =  x  tane, 
where  t^^  is  the  distance  between  successive  height  values  and  9  is  the 
instrumental  setting  for  the  photoelectron  escape  angle  with  respect 
to  the  plane  of  the  surface,  and  h^^^^  refers  to  the  heights  of  points 
along  the  analyzer  axis.  A  "length"  of  surface  along  the  direction  of 
analysis  is  computed  as  the  vector  distance  between  points  h.  and 
^i+n'  ^^^^  surface  length  value  is  used  to  properly  weight  the 
parameters  associated  with  point  i.  The  vector  distance  L  is 
calculated  using  the  equation 


30 


^T„  +  (h.  ,  -  h,  ) 


2 

—  n     1 

i+1 


This  normalization  process  rests  on  the  assumption  that  a  larger 
surface  area  pixel  (greater  value  of  L)  emits  a  proportionately  larger 
number  of  electrons.  All  the  points  in  the  vector  or  matrix  are 
scanned  in  this  manner,  and  Eg  is  found  by  carrying  out  a  weighted  sum 
over  all  unblocked  points  and  dividing  by  the  sum  of  surface  length 
values,  representing  the  total  surface  length  along  that  particular 
vector . 

In  order  to  account  for  the  change  in  the  angle  of  photoelectron 
escape  from  the  instrumental  setting  9  to  the  average  actual  value 
<e'>,  the  value  of  9^'  at  each  point  is  approximated  (Figure  3-4). 
The  slope  at  a  point  on  the  surface  is  estimated  to  be  the  difference 
in  two  successive  points  divided  by  the  distance  x  separating  them. 


X 

In  this  manner. 


9^'    =  9   +   <j)., 


where 


*.    =  Arctan(h^  -  h^^/t^). 

A  weighted  average  over  all  points  previously  found  to  be  unshaded  is 
carried  out  for  9^'   and  sin  9^',   to  yield  <9'>,   <sin  9'>,   the 
distribution  of  Oin  9'>  over  the  range   [0,1],  and  the  value  of  the 
ratio  Oin  9'>/sin  9.      It  should  be  noted  that  sin  <9'>  is  not  in 
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Figure  3-4.   Illustration  of  the  relationship  between  9,  '^^,   and  9[ 
for  the  three  possible  cases. 
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general  equal  to  Oin  e'>.  The  assumption  is  made  that  a  general  idea 
of  the  average  depth  analyzed,  which  has  been  shown  to  be  directly 
related  to  sin  9  for  an  ideally  flat  surface,  may  be  obtained  by 
looking  at  <sin  9'>,  which  represents  the  average  taken  over  the 
entire  rough  surface.  The  program  does  output  the  distribution  of 
sin  9|  values  as  well,  in  order  that  a  more  complete  analysis  may  be 
carried  out  if  called  for.       '?'  \ 

In  order  to  determine  how  the  population  of  electrons  which  can 
reach  the  analyzer  entrance  slits  varies  as  a  function  of  take-off 
angle  and  surface  roughness,  it  is  necessary  to  take  two  factors  into 
account.  First,  the  area  of  each  pixel  which  is  unshaded  must  be 
determined.  Second,  the  depth  below  the  surface  from  which  the 
photoelectrons  may  emerge  without  being  inelastically  scattered  must 
be  found  for  each  pixel.  These  two  values  may  then  be  multiplied 
together  for  each  pixel  to  yield  its  effective  sampling  volume. 
Finally,  the  volume  values  over  all  the  unshaded  pixels  are  summed. 
The  surface  roughness  modeling  program  then  reports  these  values 
normalized  so  that  near  normal  exit  at  9  =  89°  yields  an  electron 
population  of  unity.  The  routine  for  determining  the  normalized 
electron  population  of  a  surface  while  taking  surface  topography  into 
account  has  little  precedent.  It  provides  a  valuable  means  of 
comparing  the  normalized  peak  areas  obtained  in  XPS  as  a  function  of 
escape  angle  with  expected  results  determined  from  the  surface 
roughness  model.  A  literature  review  indicates  that  little  progress 
has  been  reported  along  these  lines  of  investigation.  However,  in 
order  to  adequately  account  for  the  effects  of  surface  topography,  and 
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thus  come  one  step  closer  to  a  truly  quantitative  XPS  methodology, 
such  investigations  are  called  for  to  provide  critical  information 
about  electron  populations  and  the  factors  which  go  into  their  makeup. 

A  means  of  presenting  an  alternative  method  of  characterizing  the 
nature  of  the  surface  roughness  is  provided  as  well.  The  program  has 
the  capability  of  calculating  the  autocorrelation  vector  for  each 
height  vector  taken  along  the  axis  of  analysis.  These  vectors  are 
then  summed  over  the  entire  matrix  to  give  an  effective 
autocorrelation  vector  for  the  height  matrix  in  the  direction  which 
the  photoelectrons  must  travel  to  reach  the  analyzer  slits.  This 
approach  has  been  applied  to  modeling  the  surface  roughness  effect  on 
the  AES  signal  (60),  where  it  affects  both  the  excitation  and 
detection  processes.  Other  researchers  (61,52)  have  also  made  use  of 
autocorrelation  functions  to  characterize  the  surface  roughness 
profile.  If  one  were  to  input  an  experimental  roughness  vector,  this 
function  would  provide  information  about  any  periodicities  in  the 
data.  Also,  the  autocorrelation  length — the  average  distance  over 
which  the  structure  is  correlated — can  be  used  as  a  measure  of  the 
effective  magnitude  of  the  topography;  rough  surfaces  are  in  general 
less  correlated  over  large  distances  than  smooth  surfaces. 

An  IMSL  subroutine  is  used  in  the  program  to  perform  the  actual 
autocorrelation  calculation.  For  a  given  height  vector  in  the 
x-direction  (i.e.,  along  the  analyzer  axis),  the  average  height  <h>  is 
first  calculated.  Then  the  autocovariance  vector  A  (j )  is  calculated 
for  a  specified  number  of  point  offsets  K,  which  is  typically  about 
10$  of  the  number  of  points  in  the  height  vector.  The  equation  used 
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to  calculate  A^^Cj)  is 

A^^(j)  =  (n)  '  Z      (h  -  <h>)  (h    -  <h>)     j=1 K. 

i  =  l  ^   J 

The  variance  s  of  the  vector  is  calculated  by 

2      -1  "  2 

s'^  =  (n)  '  E  (h  -  <h>)'^. 

i  =  1 

2 

The  autocorrelation  vector  is  then  given  by  A  (j)/s  . 

Finally,  the  program  attempts  to  estimate  the  total  surface  area 
for  a  matrix  of  digitized  height  values.  The  approach  is  to  divide 
the  surface  into  a  series  of  triangles  joining  three  near-spaced 
points,  calculating  the  area  of  these  triangles,  and  summing.  Better 
surface  area  measurement  techniques  are  readily  available  in  cases  for 
which  the  quality  of  the  data  justifies  a  more  quantitative  measure. 
It  is  simply  designed  to  give  another,  qualitative  measure  of  the 
magnitude  of  surface  roughness. 

The  program  thus  allows  an  evaluation  of  the  two  major  factors 
related  to  a  sample's  surface  topography  which  are  responsible  for 
varying  the  XPS  signal.  The  signal  can  be  decreased  in  magnitude  by  a 
certain  region  of  the  sample  blocking  electrons  emanating  from  another 
region  in  the  direction  of  the  analyzer  entrance  slits.  In  addition, 
electrons  from  unshaded  regions  may  leave  the  surface  at  an  angle 
quite  different  from  the  instrumental  photoelectron  escape  angle. 
This  phenomenon  leads  to  a  change  in  the  effective  sampling  depth  of 
the  photoelectrons,  which  is  related  to  <sin  e'>.  The  interplay  of 
these  two  factors  to  create  the  overall  electron  population  as 
affected  by  the  surface  topography  is  determined  by  the  electron 
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population  modeling  subroutine.  The  auto-correlation  behavior  of  the 
surface  matrix  along  the  axis  of  analysis  is  also  available  as  a  means 
of  characterizing  the  surface  topography.  Finally,  provision  is  made 
for  approximating  the  total  unshaded  surface  area  represented  by  the 
input  matrix  of  height  values. 
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CHAPTER  4 
TESTING  AND  APPLICATIONS  OF  THE  SURFACE  ROUGHNESS  MODEL 


In  this  chapter,  the  various  analytical  roughness  patterns  used 
to  evaluate  the  performance  of  the  surface  roughness  model  are 
discussed  and  the  results  obtained  are  compared  to  the  expected 
values.  In  addition,  several  applications  of  the  program  output  to 
experimental  data  are  presented.  Lastly,  the  experimental 
investigations  into  the  magnitude  and  nature  of  surface  topography 
effects  are  presented  and  discussed  in  light  of  the  theoretical 
results.  ^  "■  '  ^  ^  - 

Testing  of  Roughness  Program 

Two  functions  which  have  been  applied  several  times  previously  to 
the  analysis  of  possible  surface  topography  effects  are  the  sinusoidal 
function  and  a  generalized  grating  function.  Both  functions  are  one- 
dimensional  in  their  surface  height  variation;  in  fact,  all  known 
previous  studies  have  treated  the  surface  roughness  problem  as 
effectively  a  unidirectional  one,  with  one  exception.  De  Bernardez 
et  al .  (53)  were  interested  in  modeling  the  effects  of  damage  from  a 
static  electron  or  ion  beam,  and  so  they  used  an  axially  symmetric 
cone-shaped  surface.  The  rationale  behind  one-dimensional  modeling 
arises  from  the  geometric  relationship  between  the  sample  surface  and 
the  electron  energy  analyzer  entrance  slits.  The  technique  of 
variable  angle  XPS  (VAXPS)  does  not,  and  in  fact  cannot,  make  use  of 
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an  attractive  potential  between  the  sample  and  the  analyzer  entrance 
3lit3.  If  such  a  potential  were  to  exist,  then  photoelectrons  from  a 
large  range  of  escape  angles,  and  thus  a  large  range  of  depths,  would 
be  energy-analyzed  and  detected  in  the  species  peak.  In  effect,  depth 
and  kinetic  energy  resolution  would  be  severely  degraded.  Thus,  since 
the  axis  of  sample  rotation  in  VAXPS  is  perpendicular  to  the  direction 
which  the  emergent  photoelectrons  must  possess  to  enter  the  energy 
analyzer,  essentially  only  one  dimension  of  the  surface  roughness— 
that  dimension  perpendicular  to  the  axis  of  rotation  and  in  the  plane 
of  the  sample  surface-has  any  effect  on  the  resultant  XPS  signal  in 
the  limit  of  an  inf initesimally  narrow  analyzer  entrance  slit. 

For  both  the  sinusoidal  and  grating  functions,  the  shading 
fraction,  vector  length,  and  values  related  to  the  photoelectron 
escape  angle  can  be  predicted  algebraically  as  a  function  of  the  ratio 
of  amplitude  to  wavelength  for  the  periodic  function.  A  definition  of 
these  parameters,  as  well  as  their  method  of  calculation,  was 
presented  in  the  previous  chapter.  The  ability  to  compare  program 
results  with  the  analytic  solutions  made  these  functions  excellent 
candidates  for  testing  the  accuracy  of  the  surface  roughness  model's 
output. 

Sinusoidal  Function  Results      f f?\  f'-^ 

For  the  sinusoidal  function  S(x)  given  by 

■  *.  *   S(x)  »  A  sin  (2Trx/A),  .■  /■■ . 

where  A  is  the  amplitude  and  A  is  the  wavelength,  the  shading  fraction 
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Fg  is  given  by 


F  = 


i"x^     *  '^^2TrA/A)  G03  (2ttx/A)]^  dx 


Jl""    v/l  +  [(2TrA/A)  cos  (2Trx/A)]^  dx 


where 


x^  -  (A/2Tr)  cos~^  [(A/2irA)  tan  9] 


and 


A  cot  9  sin  (2TrX2/A)  -  x  =  A  cot  9  sin  (2irx  /A)  -  x  -  A. 

The  instrumental  setting  for  the  photoelectron  take-off  angle  is  given 
by  9.  It  can  be  deduced  from  the  above  equations  that  no 
photoelectron  shading  occurs  if  (A/2irA)«tan  9  >  1.0.  Once  x,  is 
obtained  from  the  first  expression,  it  can  be  substituted  into  the 
second  expression  to  obtain  Xj.  These  values  can  in  turn  be 
substituted  into  the  integral  equation  to  yield  the  shading 
fraction.  Although  this  integral  form  does  not  possess  an  analytic 
solution,  it  can  be  numerically  integrated  to  any  desired  precision  to 
give  values  for  Fg  algebraically.  Calculation  of  the  parameters 
related  to  the  photoelectron  escape  angle  9  from  first  principles  is 
more  complex.  Essentially,  the  procedure  involves  the  determination 
of  the  surface  tangent  for  a  vector  element  dx,  multiplying  it  by  the 
length  of  dx,  summing  over  all  the  unshaded  vector  elements,  and 
dividing  by  the  overall  unshaded  length.  Again,  numerical  integration 
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techniques  are  required  to  obtain  an  algebraic  solution  to  the 
problem. 

Comparison  between  the  analytic  results  and  the  program  output 
for  the  one-dimensional  sinusoidal  function  can  be  seen  in  Figures  4-1 
and  4-2  for  several  different  wavelengths  (holding  the  amplitude 
constant)  over  the  range  of  escape  angles  such  that  10°  £  Q  <.  89°. 
The  curves  in  Figure  4-1  represent  the  true,  algebraically  arrived-at 
values  for  F_,  while  the  points  represent  the  roughness  program 
results  for  a  sinusoidal  vector  consisting  of  50  wavelengths  digitized 
to  5000  points.  As  A  is  increased  relative  to  the  amplitude,  the 
surface  topography  becomes  less  pronounced.  Thus,  A  =  0.1  represents 
an  analogy  to  a  very  rough  surface,  while  A  =  20.0  corresponds  to  a 
smooth,  essentially  flat  surface.   In  general,  the  agreement  is  quite 
satisfactory;  relative  errors  are  greatest  at  very  low  values  of  Fg, 
where  the  program  must  make  use  of  a  small  number  of  points  in  its 
calculation.  This  error  can  in  principle  be  reduced  by  increasing  the 
digitized  resolution.  However,  constraints  on  the  memory  available 
for  array  storage  during  execution  of  the  program  make  such  error 
reduction  difficult  to  achieve.  The  general  trend  followed  by  Fg  as  a 
function  of  escape  angle  can  be  seen.  As  9  is  increased,  more  of  the 
surface  has  the  direct  line-of-sight  to  the  analyzer  entrance  slits 
necessary  for  emitted  photoelectrons  to  be  collected.  Thus  Fg 
increases  as  0  is  increased.  At  some  cutoff  angle  9^,  the  surface  is 
fully  visible  to  the  analyzer,  with  no  shaded  areas.  At  this  and  all 
higher  escape  angles,  the  shading  fraction  is  unity.  For  the 
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sinusoidal  function,  the  cutoff  angle  9^  increases  as  the  wavelength 
decreases;  for  a  rougher  surface,  the  shading  effect  is  apparent  over 
a  greater  range  of  9  values.  Finally,  the  curves  demonstrate  that  not 
only  is  there  shading  over  a  greater  escape  angle  range  for  a  rougher 
surface,  but  the  shading  is  greater,  apparently  allowing  less  of  the 
actual  emitted  photoelectrons  to  make  their  way  to  the  analyzer 
entrance  slits.  ""    ^  ^..W-^y- 

In  Figure  4-2  are  several  curves  showing  the  variation  in  the 
ratio  R  =■  <3in  9'>/sin  9  as  a  function  of  escape  angle  9.  Here, 
<3in  9'>  represents  the  average  value  for  the  sine  of  the  true 
photoelectron  take-off  angle,  which  differs  from  the  instrumental 
setting  sin  9  for  rough  surfaces.  Results  for  the  same  four 
sinusoidal  vectors  as  those  shown  in  Figure  4-1  are  illustrated  here; 
again,  the  curves  represent  an  analytic  determination  of  R,  while  the 
points  correspond  to  program  output.  Agreement  is  in  general  better 
between  output  and  true  values.  This  value  is  not  quite  as  sensitive 
to  the  use  of  small  numbers  of  points  in  its  calculation  as  the 
shading  fraction.  For  a  perfectly  smooth  surface,  R  is  equal  to 
unity.-  Over  much  of  the  range  of  escape  angles,  R  is  close  to  unity 
for  the  analogous  smooth  vector  where  A  =■  20.0.  Only  below  the  cutoff 
angle  9-  does  it  begin  to  show  an  increase.  For  the  lower  values  of  A 
corresponding  to  rougher  surfaces,  a  general  decrease  in  R  is  observed 
as  9  is  increased,  with  the  curve  leveling  out  beyond  9^.  Note  that 
at  low  values  of  9,  the  ratio  is  in  general  greater  than  unity,  while 
the  reverse  is  true  at  high  values  of  9.  This  observation  has 
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ramifications  for  the  determination  of  the  overall  sampling  volume  to 
be  discussed  below.  '/•.'■; 

Grating  Function  Results 

The  second  analytic  test  function,  a  grating  function,  was  chosen 
for  several  reasons.     It  has  been  previously  studied  by  several 
researchers.     Also,   its  use  addresses  the  ability  of  the  surface 
roughness  program  to  deal  with  discontinuities — points  of 
indeterminate  slope—in  the  topography  of  a  surface.     Finally,  this 
function  is  not  in  general  isotropic  in  its  roughness  along  the 
direction  of  measurable  photoelectron  exit.     It  has  a  different 
surface  area  facing  toward  the  analyzer  from  the  side  facing  away  from 
the  analyzer,  a  category  into  which  certain  real  world  surfaces  may 
fall.     For  the  one-dimensional  grating  function,  analytic  solution  can 
be  accomplished  as  well.     The  variable  parameters  are  the  primary  and 
secondary  blaze  angles  p  and  Y,  respectively;   the  baseline-to-peak 
distance  is  fixed  at  2.0.     The  wavelength  A  is  then  given  by 

A  =  2.0   (sin  p   +  sin  Y). 

If  the  photoelectron  escape  angle  is  greater  than  the  primary  blaze 
angle  p,   then  the  surface  is  fully  illuminated  and  Fg  =  1 .     Otherwise, 
one  must  use  a  geometric  approach  to  the  determination  of  Fg. 
Recognizing  that  the  function  has  a  peak-to-trough  amplitude  of  2.0, 
the  initial  cutoff  point  x^   is  given  by  , 


X     =«  2.0/tan  p, 


■.'■e- 
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By  drawing  a  series  of  right  triangles  and  using  elementary  analytic 

geometry,  the  second  cutoff  point  X2  is  determined  as 

X  +  A  +  2  cot  9  (tan  Y/tan  p) 
^2  1  +  cot  9  tan  Y      * 


The  equation  for  Fg  turns  out  to  be 


F^  -  (1/L)   (x_  -  X,  )^  +  (h^  -  2.0)^  , 

s  a      I         X 


where  the  surface  length  L  across  one  wavelength  is  given  by 


L  =  2.0  (1/3in  p  +  1/sin  Y) 


and  the  height  h^^  at  point  X2  equals 


h     =  2.0   +  2.0   (tan  Y/tan   p)   -  x     tan  Y. 


The  same  basic  procedure  used  for  calculating  the  parameters  related 
to  the  angle  of  photoelectron  escape  in  the  case  of  the  sinusoidal 
function  is  again  used  for  the  grating  function. 

Comparison  between  analytic  results  (curve)  and  program  output 
(points)  is  shown  for  the  grating  function  in  Figures  ^-3  and  4-4. 
For  the  program's  input  height  vector,  20  wavelengths  were  digitized 
to  give  1000  points.  As  the  blaze  angles  increase,  the  effects 
associated  with  increased  surface  roughness  increase  as  well.  The 
same  general  trends  observed  for  Fg  (Figure  4-3)  and  R  (Figure  4-4)  as 
a  function  of  photoelectron  escape  angle  in  the  sinusoidal  functional 
case  are  reproduced  here  as  well,  lending  an  air  of  generality  to  such 
behavior.  In  this  case,  the  cutoff  angle  is  easily  seen  to  equal  the 
primary  blaze  angle;  only  below  this  angle  can  shading  occur. 
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Agreement  between  the  program  output  and  the  analytic  values  is  better 
for  the  grating  function  than  the  sinusoidal  function,  primarily 
because  the  pixel  slope,  upon  which  several  key  calculations  are 
based,  does  not  vary  continuously  in  this  case,  but  can  take  at  most 
only  one  of  two  discrete  values.  Note  also  that  the  slope  of  the 
curves  near  9^  is  sharper  in  this  function.  This  factor  arises  from 
the  periodic  point  discontinuities  in  the  roughness  vector.  The 
greater  the  blaze  angles  p  or  Y.  the  more  the  surface  vector  behavior 
deviates  from  that  for  a  smooth  surface;  increasing  p  causes  an 
increase  in  9^,  while  increasing  Y  reduces  Fg  and  increases  R. 
Particulate  Surface  Model  Results  and  Applications 
The  surface  roughness  program  was  first  used  to  test  for  any 
difference  in  the  effects  of  surface  topography  on  two  different 
surfaces.  As  experimental  data,  powder  size  distributions  derived 
from  scanning  electron  micrographs  were  compiled.  In  order  to 
transfer  these  data  to  a  form  suitable  for  input  into  the  program,  the 
particulate  surface  model  (PSM)  was  developed.  The  PSM  is  a  simple 
approximation  to  the  actual  surface  topography  vector.  At  present,  it 
is  strictly  one-dimensional  in  its  approach.  The  inherent  assumptions 
are  as  follows:  '  .-  ,> 

1)  Surface  roughness  effects  on  the  determination  of  the  shading 
fraction  and  the  effective  escape  angle  parameters  are 
predominant  primarily  along  the  analyzer  axis. 

2)  Given  a  particle  size  distribution,  it  is  assumed  that  the 
particles  are  close-packed  along  the  analyzer  axis. 


59 


3)   The  maximum  height  of  each  particle  above  a  surface  baseline 

is  modeled  by  taking  a  random  fraction  of  its  radius. 
^O   The  particles  are  assumed  to  be  spherical,  or—to  be  more 
precise— to  give  a  circular  cross-section  along  the  "slice" 
of  surface  which  is  analyzed  by  the  program. 
The  resulting  model  surface  resembles  that  shown  in  Figure  H-5   (a). 
This  surface  is  then  digitized  to  give  a  vector  of  height  values  in 
Figure  4-5  (b).  Given  a  vector  representing  the  randomized  fraction 
of  each  particle's  radius,  the  height  h  as  a  function  of  surface 
displacement  p  for  a  particle  with  diameter  R  is  given  by 


h  =  0.5  Ad  •  (R  -  d)  -  m. 


where 


d  -  R/2  -  F 


q  =  C/2  -  p 


C  -  AF  •  (R  -  F) 


m  =  R/2  -  F 


F  »  (R/2)  •  RND  (0,1), 

as  seen  in  Figure  4-6.  The  function  RND  (0.1)  represents  a  random 
number  generator  with  output  s  such  that  0  <  s  <  1.  The  displacement 
P„  from  the  edge  of  particle  n  is  related  to  the  total  displacement  x 


n 


n-1 
X  -  E  F  . 

i-1   ^ 


■a 


Figure  4-5. 


Illustration  of  a  portion  of  surface  height  vector 
obtained  from  the  PSM:  a)  construction  of  the  vector 
from  a  close  packing  of  particles;  b)  result  of 
digitization,  yielding  a  suitable  input  vector. 
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Particulate   Surface  Model 


Digitized    Height    Values 
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The  sample  powders  whose  surface  topography  effects  are  to  be 

compared  are  homogeneous  Pb  ISEs  manufactured  from  a  co-precipitation 

of  PbS  and  AggS.     Two  slightly  different  powder  preparation  techniques 

have  been  used.     The  general  experimental  procedure  is  to  dissolve 

Pb(N0o)2   (MCB  Chemicals,   Reagent  grade)  and  AgNOo    (MCB  Chemicals, 

Reagent  grade)  in  de ionized  water  in  amounts  such  that  2n„     =  n.    .     A 

Pb    Ag 

0.1  M  NagS  solution  is  then  prepared.   In  the  first  procedure,  the 
NagS  solution  is  added  in  excess  to  the  mixture  of  Pb^"^  and  Ag"^  salts, 
while  stirring  at  room  temperature.  The  second  procedure  is  similar; 
the  mixture  of  Pb^*  and  Ag"*"  salts  is  added  to  an  excess  of  Na2S 
solution.  The  resultant  precipitates  are  isolated  and  washed  several 
times  with  0.1  M  HNOo  to  remove  any  excess  sulfide.  The  powders  are 
finally  washed  with  deionized  water,  filter  dried  overnight,  dried  in 
a  IIO"  oven  for  24  hours,  and  stored  in  a  dessicator. 

The  scanning  electron  microscopy  (SEM)  photographs  shown  in 
Figure  4-7  were  taken  on  a  JEOL  JSM-35CF  microscope  operating  at  25  kV 
accelerating  potential  and  20,000X  magnification.  These  micrographs, 
which  are  representative  of  those  for  each  powder,  correspond  to 
powders  prepared  using  procedure  1  (Figure  4-7  (a))  amd  procedure  2 
(Figure  4-7  (b))  described  above.     -  .. 

Certain  general  differences  between  the  two  powders  may  be 
noted.  The  particles  from  the  first  procedure  appear  to  be  larger  as 
a  rule  and  more  widely  scattered  in  size.  These  indications  are  borne 
out  by  an  analysis  of  the  particle  size  distributions,  shown  in 
Figure  4-8.  The  skewedness  of  the  distribution  in  Figure  4-8  (b )  may 
be  due  in  part  to  the  lower  resolution  of  the  SEM  photo  in 
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Figure  4-7. 


Scanning  electron  photomicrographs  of  two  PbS/Ag2S 
powders:  a)  preparation  by  procedure  1,  and  (b) 
preparation  by  procedure  2;  both  of  which  are  described 
in  the  text.  The  bar  at  the  bottom  of  each  photograph 
corresponds  to  1  pm. 
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Figure  4-8. 


Particle  size  histograms,  showing  the  number  of 
particles  as  a  function  of  particle  diameter,  derived 
from  the  SEM  photomicrographs  shown  in  Figure  4-7:  a) 
histogram  of  particle  sizes  from  the  SEM  photo  in  Figure 
4-7  (a);  b)  histogram  of  particle  sizes  from  the  SEM 
photo  in  Figure  4-7  (b). 
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Figure  4-7  (b).  The  presence  of  any  smaller  particles  would  only 
strengthen  the  argument  that  the  average  particle  size  differs  in  the 
two  photomicrographs.  In  the  particle  size  distributions  determined 
from  the  two  micrographs,  photo  0005  (procedure  2)  shows  a  mean 
diameter  of  0.13  ym,  while  the  particles  in  photo  0009  (procedure  1) 
possess  a  mean  diameter  of  0.15  ym.  If  we  make  the  assumption  that 
each  particle  population  follows  a  Gaussian  distribution  in  particle 
diameter,  then  the  hypothesis  that  the  second  mean  diameter  is  not 
greater  than  the  first  mean  diameter  can  be  rejected  at  a  99.5? 
confidence  level.  The  histogram  shown  for  photo  0005  probably 
underestimates  the  sample  standard  deviation  and  overestimates  the 
sample  mean.  However,  this  hypothesis  rejection  holds  even  if  we  use 
the  larger  of  the  two  sample  standard  deviations  as  a  measure  of  the 
population  standard  deviation,  as  opposed  to  pooling,  or  if  we  lower 
our  estimate  of  the  first  sample  mean. 

These  powder  distributions  were  then  used  as  the  input  vector  for 
the  PSM  and  an  analysis  of  the  comparative  surface  roughness  effects 
was  carried  out.  The  comparison  results  are  shown  for  the  shading 
fraction  Fg  in  Figure  4-9,  and  for  the  ratio  R  =  <sin  e'>/sin  e  in 
Figure  4-10.  It  is  immediately  apparent  that  the  difference  in 
particle  diameters  is  not  sufficient  to  cause  a  noticeable  difference 
in  their  XPS  signals  attributable  to  the  differing  surface 
topography.  A  statistical  comparison  between  the  results  for  both 
powders  supports  this  statement.  However,  both  powders  are  expected 
to  show  appreciable  surface  roughness  effects  in  their  respective 
VAXPS  signals,  especially  at  low  values  for  the  photoelectron  escape 
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Fijmre  4-9.   Plots  of  the  shading  fraction  F  as  a  function  of 

photoelectron  take-off  angle  9  (deg),  derived  from  the 
SEM  photomicrographs  shown  in  Figure  4-7:  a)  F  versus 
9  plot  for  the  powder  depicted  in  Figure  4-7  (a);  b)  F 
versus  9  plot  for  the  powder  depicted  in  Figure  4-7  (b;. 
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Figure  4-10. 


Plots  of  the  ratio  R  =  Oin  9'>/sin  9  as  a  function  of 
take-off  angle  9  (deg),  derived  from  the  SEM 
photomicrographs  shown  in  Figure  M-7:  a)  plot  of  R 
versus  9  for  the  powder  depicted  in  Figure  4-7  (a);  b) 
plot  of  R  versus  9  for  the  powder  depicted  in  Figure  4-7 
(b). 


■'■mm-- 


73 


2.5  ■ 
2- 

1.5- 

I- 

.5- 


PHOTO  8009 


2Q  THETA  60 


90 


2-5" 
2- 

1.5  • 

!•• 

-5- 


PHOTO  0005 


'■•-3.  t''i   '•"*^ 

i::3  * 


••^. 


I  ^^  L'  r^ 


■+ 1- 


3®  THETA  60 


9^ 


7M 


angle.  Also,  note  that  no  sharp  transitions  can  be  observed  in  the 
curves.  As  has  been  previously  seen,  such  transitions  arise  from 
periodicities  in  the  surface  roughness  vector;  since  the  PSM  produces 
an  essentially  aperiodic  roughness  vector,  as  most  randomly  rough 
surfaces  would  be  expected  to  produce,  this  observation  is  expected. 

Thus,  the  PSM  results  provide  several  useful  pieces  of 
information.  They  more  closely  approximate  the  expected  behavior  of 
Fg  and  R  of  randomly  rough  surfaces  as  a  function  of  9.  In  the 
present  case,  the  results  allow  the  conclusion  that  the  two  methods  of 
Pb  ISE  preparation  do  not  produce  different  surface  topographical 
effects.  This  factor  is  of  importance  to  the  discussion  of  Pb  ISEs  in 
a  later  chapter.  Finally,  and  of  more  general  interest,  the  PSM 
provides  a  means  of  estimating  the  surface  roughness  effect  of  a 
powder  sample  using  particle  size  distributions  readily  available  from 
SEM  photomicrographs. 

Experimental  Surface  Roughness  Studies 

An  analysis  of  the  effects  of  varying  surface  roughness  has  been 
carried  out  on  two  surface  chemical  systems  of  relative  simplicity. 
The  first  system  consists  of  sputter-cleaned  gold  surfaces.  These 
surfaces  are  essentially  free  of  significant  levels  of  contamina- 
tion. Therefore,  their  normalized  peak  areas  as  a  function  of 
photoelectron  escape  angle  can  be  used  to  characterize  the  instru- 
mental response  of  the  XPS  instrument.  The  Au  ^fj/2   Peak  at  83. 80  eV 
is  also  commonly  used  to  steuidardize  the  instrumental  binding  energy 
response.  The  second  system  is  a  carbon  sample  with  measurable 
surface  oxidation.  Such  samples  are  more  readily  available  than  gold 
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samples,  and  they  provide  easily  observable  evidence  of  the  contribu- 
tion of  oxidized  species  to  the  overall  carbon  concentration.  Gold  is 
relatively  inert,  and  the  Au  ^'^'t/2   P®^*^  does  not  show  the  presence  of 
oxides.  The  value  of  being  able  to  resolve  oxidized  and  unoxidized 
carbon  arises  from  the  information  about  the  electron  sampling  depth 
available  from  such  measurements.  If  the  sampling  depth  is  relatively 
small,  then  the  carbon  peak  envelope  shows  a  significant  contribution 
from  oxidized  carbon  species.  For  greater  sampling  depths,  the  bulk 
graphitic  carbon  contribution  is  predominant. 
Results  of  Gold  Studies 

Initial  studies  using  gold  samples  focused  on  the  determination 
of  the  response  function  of  the  XPS  instrument.  In  order  to  minimize 
surface  topography  and  vertical  inhomogeneity  effects,  a  gold  foil 
sample  (Alfa  Chemicals,  99.9+$  purity,  0.25  mm  thick)  was  sputter 
cleaned  until  the  residual  carbon  contamination  peak  was  minimized  in 
intensity.  Source  power  was  300  W  (15  kV,  20  mA),  and  a  small  solid 
angle  of  detection  (LO  magnification  mode)  was  used.  The  peak  area  of 
the  Au  ^f'7/2  ^'^'^   '*^5/2  P®^'^  doublet  was  determined  after  subtraction 
of  the  background  signal  (Figure  4-11)  using  an  iterative  technique 
developed  by  Shirley.  Peak  area  is  calculated  using  the  equation 

A  (count)  •  Energy  step  size  (eV/ch) 

A  =  _i _^ 

N  ( count -eV/sec)  Dwell  time  (sec/ch)         ' 

where  A-j.  represents  the  sum  of  counts  in  all  the  channels  in  the  peaks 
after  background  subtraction,  the  energy  step  size  is  an  instrumental 
setting  defining  the  ultimate  energy  resolution  of  the 
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spectrum,  and  the  dwell  time  is  given  by 

Dwell  time  _  Scan  time  (sec)  »  #  of  sweeps  •  Step  size  (eV/ch) 
(sec/ch)  Energy  scan  range  (eV) 

Since  the  X-rays  used  are  not  strictly  monochromatic,  satellite  lines 
are  visible  in  this  spectrum.  The  satellites  do  not  interfere  with 
the  peak  area  determination  in  this  instance.  The  area  was  determined 
for  a  series  of  photoelectron  escape  angle  settings  ranging  from  90° 
to  5°.  This  area  was  then  normalized  by  dividing  by  the  largest  peak 
area  value  for  the  range  of  escape  angle  settings.  The  results  of 
these  variable  angle  plots  are  shown  in  Figure  4-12  for  both  Mg  Ko  and 
Al  Ka  radiation.  The  decrease  in  peak  area  at  9  =  90°  when  using  the 
Mg  anode  results  from  the  arrangement  of  the  anode  faces.  The  anode 
tip  is  wedge-shaped,  with  the  upper  side  aluminum-coated  and  the  lower 
side  magnesium-coated.  For  9  =  90°,  the  sample  surface  is 
approximately  even  with  the  tip  of  the  wedge  and,  thus,  does  not 
receive  the  total  X-ray  flux  when  using  the  Mg  anode.  Since  data 
taken  at  90°  were  to  be  used  in  the  surface  roughness  studies,  the  Al 
anode  was  used  throughout  the  investigation. 

In  order  to  ascertain  the  basic  effects  of  increasing  surface 
topography  on  the  results  of  VAXPS  experiments,  two  other  gold  samples 
were  prepared  and  analyzed.  As  a  model  for  the  behavior  of  a  very 
rough  surface,  a  gold  powder  (Alfa  Chemicals,  99.995J  purity,  -20  mesh 
powder  size)  was  sprinkled  onto  the  surface  of  conducting  copper  tape, 
which  was  then  affixed  to  the  sample  holder  using  colloidal  silver 
adhesive.  In  this  manner,  a  rough  gold  surface  could  be  examined 
without  any  of  the  sample  charging  effects  due  to  the  use  of 
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common  adhesive  tape.  In  addition,  this  powder  was  pressed  at  7  psi 
into  a  pellet  to  be  analyzed  as  well.  This  pellet  was  analyzed  to 
determine  the  general  surface  roughness  behavior  of  a  sample  prepared 
using  a  common,  simple  hand  press  which  is  often  used  in  sample 
preparation  for  routine  XPS.  Operating  conditions  were  the  same  as 
above,  and  the  VAXPS  curves—plots  of  total  doublet  peak  area  as  a 
function  of  photoelectron  escape  angle  9~for  the  foil,  pellet,  and 
powder  samples  are  compared  in  Figure  4-13.  In  this  graph,  all  the 
peak  areas  were  normalized  only  to  account  for  any  differences  in  the 
peak  collection  times.  Thus,  the  overall  peak  areas  as  a  function  of 
escape  angle  are  compared  for  the  three  samples.  Several  noteworthy 
changes  are  apparent  as  the  surface  roughness  increases  in  the  order 
foil  <  pellet  <  powder.  First,  the  doublet  peak  area  drops  off 
significantly  at  high  values  of  9,  then  actually  increases  as  9  is 
decreased.  As  the  surface  roughness  increases,  the  escape  angle  where 
the  maximum  peak  area  is  obtained  decreases.  In  addition,  in  the 
range  of  low  9,  the  peak  area  is  greater  for  rougher  surfaces.  These 
two  factors  combine  to  lower  the  range  of  peak  area  values  over  the 
range  of  escape  angle  settings.  Whereas  the  peak  area  for  the  foil 
increases  by  2.5  orders  of  magnitude  over  the  range  of  angles,  only  an 
order  of  magnitude  increase  is  seen  for  the  pellet  sample,  followed  by 
a  decrease  at  higher  take-off  angle  settings.  For  the  powder  sample, 
the  peak  area  increases  only  by  a  factor  of  3  before  decreasing  at 
higher  9.  Obviously,  this  observation  has  important  ramifications 
when  one  is  varying  the  angle  setting  in  order  to  maximize  the  XPS 
signal.  In  addition,  the  tacit  assumption  that  surface  roughness 
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effects  are  negligible  with  pressed  pellet  samples  is  shown  to  be 
false,  due  to  the  difference  in  VAXPS  signal  between  gold  foil  and 
compressed  powder  pellet.  This  point  forms  the  basis  of  a  study  of 
the  effect  of  pressure  on  the  XPS  signal,  presented  below.  Finally, 
the  effects  of  surface  roughness  may  be  quite  dramatic,  as  shown  in 
the  VAXPS  response  curve  for  the  gold  powder  sample. 
Results  of  Carbon  Studies:  Model  Comparison 

The  first  series  of  experiments  investigating  a  graphitic  carbon 
surface  were  designed  to  test  the  predictions  of  the  electron 
population  model  which  was  developed  as  a  routine  in  the  surface 
roughness  program.  As  discussed  above,  this  model  outputs  an 
approximation  to  the  normalized  VAXPS  curve  as  a  function  of  the 
topography  of  the  input  surface  height  vector  or  matrix.  The 
predictions  of  the  model  about  the  changes  in  the  VAXPS  curve  are 
compared  here  with  some  experimental  studies  of  carbon  surfaces  of 
various  roughness.        -«  ,  ..  ..„  ,, 

Experimentally,  the  samples  studied  included  a  carbon  foil 
(Goodfellow  Metals,  99.8+$  purity,  0.125  mm  thickness)  which,  although 
smooth,  is  not  as  smooth  as  the  gold  foil  sample.  Also,  studied  were 
two  carbon  powder  samples  (Alfa  Products,  99.5$  total  purity,  -325 
mesh  and  -20  +60  mesh).  For  the  -325  mesh  powder,  all  particles  are 
smaller  than  44  pm  in  diameter;  for  the  -20  +60  mesh  powder,  particle 
diameters  range  from  250-841  pm.  However,  for  the  latter  powder,  the 
particles  flake  together,  forming  smooth  sections  of  surface.  Both 
powders  were  sprinkled  onto  conductive  copper  tape,  then  affixed  to 
the  sample  holder  using  colloidal  silver.  The  carbon  powder 
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(-325  mesh)  sample  forms  a  surface  with  appreciable  topography;  it 
serves  to  act  as  a  very  rough  surface.  The  carbon  flake  (-20  +60 
mesh)  sample  consists  of  disjointed  sections  of  smooth  surface;  its 
surface  roughness  is  intermediate  between  the  powder  and  the  carbon 
foil.  The  instrumental  operating  conditions  for  data  collection  were 
the  same  as  in  the  case  of  the  gold  samples.  For  all  three  samples,  a 
series  of  VAXPS  experiments  was  carried  out  by  varying  the 
instrumental  setting  of  the  photoelectron  take-off  angle  9  in  the 
range  [5", 90°].  The  C  l3^/2  P®^'^  ai'eaa  were  normalized  in  each  case 
so  that  at  the  angle  where  the  peak  area  was  maximum,  the  area  had  a 
value  of  unity.  The  electron  population  model  results  over  the  same 
range  of  9  were  compared  for  a  series  of  wavelengths  using  the  one- 
dimensional  sinusoidal  function  and  matching  the  best  approximations 
to  the  experimental  results.  The  curves  are  compared  for  the  three 
samples  in  Figure  4-14  (carbon  foil),  Figure  4-15  (carbon  flakes),  and 
Figure  4-16  (carbon  powder).  The  average  percent  difference  between 
model  and  experimental  points  is  2.34$  for  the  carbon  foil,  4.57J  for 
the  carbon  flakes,  and  1.81$  for  the  carbon  powder.  The  carbon  foil 
shows  a  response  very  close  to  an  ideally  flat  surface.  In  fact,  the 
percent  difference  between  the  carbon  foil  and  a  model  flat  surface  is 
2.42$.  Thus,  the  carbon  foil  serves  to  act  as  an  experimental  flat 
surface  standard.  Note  the  good  agreement  between  model  and 
experimental  results.  As  is  expected,  the  wavelength  of  the 
sinusoidal  function  which  best  matches  the  experimental  curve 
increases  as  the  roughness  of  the  sample  decreases.  That  is,  a 
smoother  model  surface  vector  is  required  to  fit  a  smoother 


Figure  il-l4. 


Comparison  of  theoretical  electron  population  and 
experimental  VAXPS  response  curves  for  a  ruffle  function 
(A  =  20.0)  and  a  carbon  foil,  respectively. 
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Figure  4-15. 


Comparison  of  theoretical  electron  population  and 
experimental  VAXPS  response  curves  for  a  ruffle  function 
(A  =  10.0)  and  a  flaky  carbon  sample,  respectively. 
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Figure  4-16. 


Comparison  of  theoretical  electron  population  and 
experimental  VAXPS  response  curves  for  a  ruffle  function 
(A  =  0.2)  and  a  carbon  powder  sample,  respectively. 
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experimental  surface.  The  uneven  nature  of  the  theoretical  curve  in 
Figure  4-16  is  attributable  to  the  small  shading  fractions  for  this 
very  "rough"  surface.  The  results  are  obtained  by  averaging  over  a 
relatively  small  number  of  points  and  are,  thus,  more  subject  to 
random  error. 

Even  though  the  surface  roughness  varied  dramatically  among  these 
samples,  they  all  show  somewhat  similar  VAXPS  curves.  This  phenomenon 
is  due  to  the  offsetting  nature  of  the  two  effects  of  the  surface 
topography  on  the  magnitude  of  the  XPS  signal.  The  first  effect  is 
the  shading  fraction  Fg  variation.  For  a  rough  surface,  the  shading 
fraction  is  low,  and  a  smaller  signal  is  expected.  However,  for  this 
same  rough  surface,  the  electron  population  model  predicts  that  the 
regions  of  the  sample  which  are  capable  of  emitting  analyzable 
photoelectrons  emit  proportionately  more  photoelectrons.  An 
examination  of  Figure  4-17  helps  to  explain  this  result.  Here,  the 
relative  population  of  analyzable  photoelectrons  emitted  by  each 
surface  pixel  is  plotted  in  bar  chart  form  directly  below  that  pixel 
for  several  escape  angle  settings.  For  this  portion  of  a  hypothetical 
surface  roughness  vector,  the  electron  population  is  relatively  evenly 
distributed  for  9  -  89°.  As  e  is  decreased,  portions  of  the  surface 
are  shaded,  but  the  unshaded  portions  emit  a  greater  number  of 
photoelectrons.  At  higher  instrumental  take-off  angles,  the  decrease 
in  the  signal  due  to  a  lower  shading  fraction  is  more  important,  and 
an  overall  lowering  of  the  signal  is  observed.  At  lower  escape  angle 
settings,  the  increase  in  the  sampling  depth  is  more  important  than 
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Figure  4-17. 


Plots  of  the  relative  populations  of  emitted 
photoelectrons  from  each  pixel  of  a  small  portion  of  a 
surface  height  vector  for  several  different  values  of 
the  instrumental  take-off  angle  9.     The  vertical  scale 
is  in  arbitrary  units. 
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the  Shading  fraction  effect,  with  a  resultant  relative  increase  in  the 
signal.  .:;-?/:^'"  '      /:/  ■  ''■ 

The  agreement  between  the  electron  population"^model  and  the 
experimental  carbon  sample  VAXPS  curves  has  importance  for  several 
reasons.  The  surface  roughness  program  contains  the  first  surface 
roughness  model  to  allow  arbitrary  input  of  a  surface  height  vector  or 
matrix,  and  subsequently  output  a  series  of  parameters  defining  the 
resultant  effect  on  the  nature  and  magnitude  of  the  XPS  signal.  The 
various  factors  are  combined  to  give  rise  to  the  electron  population 
model,  which  successfully  mimics  the  behavior  of  experimental  VAXPS 
curves,  again  for  the  first  time.  That  the  agreement  is  so  close- 
even  though  a  Simple  one-dimensional  modeling  function  is  used- 
supports  the  argument  that  only  the  component  of  surface  roughness  in 
the  direction  of  the  analyzer  entrance  slits  significantly  affects  the 
XPS  Signal.  More  evidence  to  support  this  assumption  is  provided 
below.  Finally,  the  idea  of  the  two'major  effects  of  surface 
topography  on  the  XPS  signal  acting  to  offset  each  other  indicates 
that  samples  with  comparable  surface  topographies  should  possess  quite 
Similar  VAXPS  curves.  This  notion  supports  the  conclusions  of  the  SEM 
powder  Size  distribution  results  discussed  above.  In  that  experiment, 
even  though  the  two  powder  size  distributions  had  different  mean 
particle  diameters,  no  distinguishable  difference  in  surface  roughness 
effects  could  be  determined.  '^■: 

The  Effect  of  Pressure  on  the  VAXPS  Curve 

When  performing  solids  analysis  on  powder  samples  using  XPS.  a 
decision  must  be  made  regarding  the  method  of  sample  preparation.  The 
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two  major  options  available  are  analysis  of  the  free  powder  attached 
to  the  sample  holder  by  conductive  tape  and  compression  of  the  powder 
into  a  pellet,  A  study  of  the  effects  of  the  pressure  used  to  form 
the  pellet  has  been  carried  out.  and  the  results  are  presented  in  this 

section.  > 

■•>»•  ■■■.•■  \ . 

The  powder  used  was  the  -325  mesh  carbon  powder  mentioned 
above.  The  carbon  foil  sample  has  been  used  as  a  comparison  standard, 
representing  in  this  case  the  flat  surface  ideal.  Two  laboratory 
presses  were  used  in  this  experiment.  A  hand-held  press  (Wilmad 
Scientific  Co.)  with  a  stated  pressure  of  7  psi  was  used  to  make  one 
pellet  for  analysis.  Such  presses  are  commonly  used  for  XPS  sample 
preparation.  To  achieve  higher  pressures,  a  Carver  Model-M  laboratory 
press  With  a  25-ton  capacity  was  utilized,  along  with  a  Beckman  KBr 
die  (Wilmad  Scientific  Co.).  Laboratory  presses  such  as  this  model 
are  widely  available  to  the  researcher,  chiefly  for  preparing  solid 
samples  for  IR  analysis.  Two  pellets  were  prepared  using  this  press, 
one  at  2000  psi  and  the  other  at  10000  psi.  Using  the  same  operating 
conditions  as  in  the  gold  studies,  a  series  of  VAXPS  experiments  was 
carried  out  for  the  three  pellets,  the  carbon  foil,  and  the  carbon 
powder  sprinkled  onto  conductive  tape.  The  results  are  summarized  in 
Figure  4-18.  where  the  absolute  total  peak  areas  of  the  C  ^s,^^   peak 
are  plotted  as  a  function  of  the  escape  angle  setting.  Several 
Observations  are  apparent  from  an  examination  of  these  curves.  Note 
that  the  peak  area  is  appreciably  lower,  especially  at  higher  exit 
angles,  for  the  free  powder  and  the  lightly  pressed  pellet.  However. 
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use  of  a  sufficiently  high  pressure  in  pellet  manufacture  has  little 
effect  on  the  nature  or  magnitude  of  the  VAXPS  curve.  In  fact,  the 
pellets  pressed  at  higher  pressures  demonstrate  behavior  virtually 
indistinguishable  from  the  carbon  foil  results. 

Of  the  two  commonly  used  powder  sample  preparation  methods  used 
for  XPS,  use  of  a  pressed  pellet  clearly  produces  a  greater  signal, 
especially  at  the  higher  take-off  angles  preferred  for  analysis  with 
the  Kratos  XSAM  800.  However,  simply  pressing  a  powder  at  a  low 
pressure  does  not  result  in  an  optimum  signal.  Use  of  significant 
pressures  yields  samples  whose  VAXPS  curves  are  similar  in  shape  and 
magnitude  to  the  essentially  flat  surface  behavior  of  a  carbon  foil. 
As  long  as  a  sufficiently  large  pressure  (i.e.,  greater  than  about 
1000  psi)  is  used,  the  resulting  response  curves  differ  negligibly. 

It  is  thus  recommended  that  powder  samples  be  pelletized  at  high 
pressures  in  order  to  maximize  the  signal,  as  well  as  to  minimize  the 
resulting  surface  roughness  effects.  The  latter  conclusion  follows 
from  the  similarity  of  the  foil  and  high-pressure  pellet  VAXPS  curves. 
Peak  Structure  Studies  on  Abraded  Surfaces 

After  characterizing  the  effect  of  surface  roughness  on  the 
magnitude  of  the  XPS  signal,  the  question  arises  as  to  whether  the 
topography  of  the  surface  plays  any  part  in  determining  the  peak 
structure.  Given  the  capability  of  XPS  to  provide  chemical  state 
information  on  the  surface  species,  the  possibility  of  how  this 
information  might  be  altered  by  surface  topographic  effects  has  been 
investigated. 
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A3  a  model  experimental  system,   a  slightly  oxidized  carbon  foil 
has  been  used.     The  foil  consists  of  graphitic  carbon  in  the  bulk  of 
the  sample;   at  the  surface,  however,   the  carbon  is  oxidized  to  a 
higher  valence.     The  resulting  C   1s^/2  P^ak  envelope  thus  shows  the 
presence  of  bulk  and  oxidized  carbon  species.     At  high  photoelectron 
take-off  angles,  the  contribution  to  the  overall  signal  from  bulk 
graphitic  carbon  is  typically  high.     Figure  4-19  illustrates  the 
makeup  of  the  asymmetric  carbon  peak  for  a  smooth  carbon  foil  at 
9  =  90°.     The  relative  contributions  and  positions  of  each  peak  to  the 
total  peak  area  are  solved  for  through  the  use  of  a  peak-fitting 
algorithm  in  the  DS800  software  accompanying  the  Kratos  XSAM  800. 
Each  peak  is  initially  assigned  a  peak  shape,   full  width  at  half 
maximum  (FWHM),  and  position  on  the  binding  energy  scale.     The 
algorithm  varies  these  settings  iteratively  in  an  attempt  to  minimize 
the  fit  parameter  given  by 
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Where  n  is  the  number  of  channels  in  the  peak  envelope,  X    is  the 

1   "  <  -   •  -"  ff-,  max 

maximum  intensity  in  the  acquired  spectrum,  X,  is  the  acquired 
spectrum  intensity  for  channel  i,  and  Sj:  is  the  synthetic  envelope 
intensity  for  channel  i.  In  Figure  4-19,  three  peaks  are  apparent:  a 
predominant  bulk  graphitic  carbon  peak,  a  smaller  peak  due  to  carbon 
singly  bound  to  oxygen,  and  an  even  smaller  peak  attributable  to 
carbon  doubly  bound  to  oxygen.  At  this  angle,  the  bulk  carbon  is 


•t)***'*"^  "  #*'  ^'-^N.         tf***^;%  ."-■ 


y- : 


jQCB 


Figure  4-19. 


Results  of  a  peak  fitting  routine  for  a  carbon  foil  at 
9  =  90",  showing  the  main  graphitic  carbon  peak  at 
283.9  eV  and  the  smaller  oxidized  carbon  peaks  at 
285.0  eV  and  286.1   eV.     The  top  curve  shows  the 
experimental  spectrum  of  intensity  versus  binding 
energy,  while  the  bottom  curve  shows  the  matching 
simulated  spectrum. 
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responsible  for  approximately  75?  of  the  total  carbon  signal.  At 
lower  electron  escape  angles,  the  relative  contribution  of  the  surface 
carbon  species  to  the  carbon  peak  envelope  is  increased,  as 
Illustrated  in  Figure  4-20.  where  the  same  sample  is  analyzed  at  a 
take-off  angle  of  50.  m  fact,  an  additional  peak  due  to  the  presence 

Of  adsorbed  CO^  at  the  very  surface  may  be  resolved  as  well.  Here. 

the  bulk  carbon  contributes  only  about  50$  to  the  makeup  of  the 

asymmetric  peak,  due  to  the  reduction  in  the  effective  IMFP 

A^  ="  X  sin  e. 

To  study  the  effects  of  surface  roughness  on  the  contribution  of 
the  graphitic  carbon  peak  to  the  overall  peak  structure,   the  carbon 
foil  samples  were  subjected  to  abrasion  in  one  direction,  creating  a 
nonisotropic  surface  topography.     Three  abraded  samples  were  prepared. 
The  first  sample  was  polished  with  crocus  cloth  and  only  very  slightly 
roughened.     The  second  sample  was  scratched  with  i.00  grit  sandpaper, 
producing  noticeable  grooves  in  the  surface  of  the  foil.     The  final 
sample  was  abraded  with  220  grit  sandpaper;    the  surface  appeared  to  be 
Significantly  roughened.     Especially  in  the  latter  two  samples,   the 
anisotropy  of  the  abrasion  process  was  visually  apparent.     The 
instrumental  operating  conditions  were  essentially  the  same  as  those 
used  in  the  previously  discussed  experimental  studies.     However,   in 
this  case,   the  percentage  of  the  total  C  ^s,,,  peak  envelope 
attributable  to  bulk  graphitic  carbon-as  elucidated  by  peak-fitting- 
is  reported  as  a  function  of  escape  angle  9. 

By  abrading  the  samples  in  only  one  direction,   two  principal 
effects  of  abrasion  on  the   sample  can  be  separately  analyzed. 
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Figure  4-20, 


Results  of  a  peak  fitting  routine  for  a  carbon  foil  at 
0  =  50.  showing  the  increased  contribution  of  the 
oxidized  carbon  peaks  at  285.3  eV,  286.4  eV.  and 
?88  3  eV  relative  to  the  main  graphitic  carbon  peak  at 
284  3  eV   The  top  curve  shows  the  experimental  spectrum 
of  intensity  versus  binding  energy,  while  the  bottom 
curve  shows  the  matching  simulated  spectrum. 
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Abrasion  serves  to  remove  the  top  layers  of  the  surface,  including  the 
oxidized  carbon  layers.  However,  the  fresh  carbon  surface  is  again 
subject  to  oxidation.  Each  sample  is  expected  to  exhibit  slightly 
different  rates  of  oxidation.  This  phenomenon  is  illustrated  in 
Figure  4-21 .  There  is  a  general  increase  in  the  percentage  of  signal 
due  to  bulk  carbon  as  the  level  of  abrasion  is  increased.  This 
observation  follows  from  the  increasingly  clean  surface  obtained  by 
greater  roughening  of  the  sample.  As  a  result,  the  bulk  carbon  signal 
contributes  to  the  overall  peak  area  to  a  greater  degree.  By 
analyzing  the  sample  so  that  photoelectrons  are  emitted  parallel  to 
the  polishing  direction,  any  effects  of  the  surface  topography  induced 
by  the  abrasion  process  are  bypassed,  and  any  change  in  the  signal  may 
be  attributed  to  surface  oxidation  effects.  When  the  sample  is 
analyzed  perpendicular  to  the  polishing  direction,  the  surface 
roughness  effects  on  the  VAXPS  signal  are  present  as  well.  The  VAXPS 
curves  for  analyses  in  the  parallel  and  perpendicular  directions  are 
shown  for  the  three  levels  of  abrasion  in  Figures  4-22,  4-23,  and 
4-24.  Note  that  the  change  in  the  shape  of  the  response  curves  is 
much  more  significant  for  the  results  of  data  collected  perpendicular 
to  the  direction  of  surface  roughness.  Most  significantly,  the  area 
percentage  of  the  bulk  carbon  signal  drops  off  at  high  take-off 
angles,  mainly  due  to  the  surface  roughness  effect.  For  the  highly 
abraded  sample  (Figure  4-24),  this  decrease  is  substantial.  The 
effects  of  surface  roughness  on  the  variation  in  the  bulk  carbon  area 
percentage  can  be  better  seen  in  Figure  4-25,  where  the  ratio  of 
perpendicular-to-parallel  area  percentage  is  shown  for  the  three 
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samples  as  a  function  of  the  instrumental  escape  angle.  In  taking 
this  ratio,  the  change  in  area  percentage  due  to  different  rates  of 
surface  oxidation  is  effectively  cancelled  out.  leaving  the  surface 
topography  effect  as  a  function  of  abrasion.  The  decrease  in  the  area 
percentage  at  high  take-off  angle  settings  is  evident.  This 
observation  provides  supporting  evidence  for  the  results  of  the 
surface  roughness  program,  it  was  stated  in  an  earlier  section  that 
the  model  predicts  a  signal  intensity  drop-off  at  high  angles  due  to 
the  smaller  effective  sampling  depth.  If  this  assertion  is  true,  then 
the  area  percentage  should  drop  similarly,  since  in  this  case  more  of 
the  total  carbon  signal  will  be  due  to  the  oxidized  carbon  species 
present  near  the  surface  and  less  will  arise  from  the  graphitic  carbon 
in  the  bulk.  An  explanation  of  the  VAXPS  curve  shape  at  lower  take- 
off angles  illustrates  the  changing  relative  contributions  of  the  two 
primary  surface  roughness  effects.  For  the  two  surfaces  with  the 
least  abrasion,  the  area  percentage  increases  with  decreasing  e  due  to 
the  greater  effective  sampling  depths  at  low  escape  angles  for  rough 
surfaces  predicted  by  the  roughness  model.  However,  if  the  surface 
roughness  is  great  enough,  the  electron  shading  effect  reduces  the 
Signal  from  all  but  the  tops  of  ridges  on  the  sample,  for  which  the 
sampling  depth  is  reduced,  for  the  same  reasons  as  in  the  high  take- 
off angle  case.  i.e..  glancing  photoelectron  exit  angles  from  the 
surface  resulting  in  lower  effective  sampling  depths  and  greater 
surface  species  contribution  to  the  signal 

Several  oo;,=iualon3  n,ay  be  drawn  free  thl3  study.  An  examination 
Of  Figure  *-2,  indicates  little  actual  ohanje  In  the  shape  of  the 


i  .  « 


.^3  *^ 


^if^i-'^lrxJini^:-7y:V\ 


A, 


118 


curves  with  increasing  abrasion  which  could  be  due  to  a  surface 
roughness  effect,  with  the  exception  of  a  minor  lowering  of  the  area 
percentage  at  higher  escape  angles  for  the  more  harshly  treated    :  . 
samples.  This  experiment  thus  supports  the  claim  made  previously  that 
only  the  surface  roughness  in  the  direction  of  detectable 
Photoelectron  exit  has  a  significant  effect  on  the  resulting  XPS 
Signal.  Characterization  of  the  surface  roughness  using  a  one- 
dimensional  model  is  justified,  provided  that  the  topography  is  shown 
to  be  directionally  isotropic.  Also,  the  effects  of  surface 
topography-Shading  and  sampling  depth  variation-on  the  VAXPS  signal. 
Which  predicted  the  shapes  of  the  response  curves,  are  here  supplied 
with  further  evidence  to  support  their  existence.  Finally, 
examination  of  Figure  4-25  indicates  that  if  a  sample  with  a  roughened 
surface  is  to  be  analyzed,  the  results  should  be  most  reliable  for 
escape  angles  between  30  -  and  60 ».  where  the  ratio  between  the  "rough" 
and  "flat"  surface  is  relatively  close  to  unity. 
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CHAPTER  5 

THEORY  AND  TESTING  OF  A  METHOD  FOR 

NONDESTRUCTIVE  DEPTH  PROFILING 


Introduction 
The  development  of  methods  for  nondestructive  depth  profiling  has 
been  reviewed  in  detail  previously  (53).  Thus,  only  a  brief  overview 
Of  this  field  is  presented.  Initial  intimations  of  the  possible  usage 
Of  XPS  data  to  obtain  information  about  the  stratification  of  the 
near-surface  region  were  published  by  Fadley  and  BergstrSm  (6i^).  The 
first  series  of  variable  angle  XPS  (VAXPS)  experiments  was  reported  by 
Fraser  et  al.  (65),  confirming  the  predictions  of  Fadley  and 
Bergstrsm.  The  methodology  for  recognizing  the  effects  of  a  thin 
surface  film  on  the  XPS  signal  was  first  presented  by  Fadley  et  al. 
(5).  Thin  film  characterization  using  VAXPS  has  since  been  used  to 
study  oxidized  silicon  (56-68).  Nb205  on  Nb  (69).  PTFE  layers  on  Au 
(70).  and  hydrocarbon  contamination  overlayers  (71).  Several 
researchers  have  developed  methodologies  specifically  for  thin  film 
analysis  (70.72.73).  These  methods,  although  different  in  their 
approach,  share  the  common  assumptions  of  a  uniform  film  thickness 
With  a  Sharp  transition  gradient  between  film  and  substrate,  m 
addition,  they  assume  that  only  two  basic  phases- the  film  phase  and 
the  bulk  phase— are  present. 

Several  researchers  have  developed  more  general  models  for 
nondestructive  depth  profiling.  Iwasaki  et  al.  (7^.)  realized  the 
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olose  relationship  between  the  Intenalty  expression  for  VAXPS  and  the 
formula  for  the  Laplace  transform  of  a  function.  They  developed  a 
method  based  on  this  realization  and  ;pplied  It  to  the  analysis  of  the 
data  Of  Praaer  et  al.  (65).  Nefedoy  (75)  developed  a  Method  of  thin 
layer  analysis  wherein  the  layer  could  follow  one  of  four  analytic 
concentration  gradients  as  a  function  of  depth,  no  experimental 
supporting  evidence  was  given.  Pijolat  and  Hollinger  (76)  assumed 
that  the  surface  could  be  broken  up  into  a  discrete  number  of  layers, 
then  calculated  and  summed  the  attenuated  intensities,  comparing  the' 
results  to  experimental  VAXPS  data.  However,  the  model  was  developed 
for  two-component  systems  only.  An  interesting  qualitative  depth 
profiling  approach,  utilizing  the  strong  depth  dependence  of  a  core 
level  peak  energy  distribution  due  to  elastic  and  inelastic 
scattering,  was  proposed  by  lougaard  and  Ignatiev  (77).  A  thin  layer 
method  based  on  Fadley's  fundamental  work  (5),  yet  allowing  for  the 
calculation  of  fractional  monolayer  coverages,  has  also  been  proposed 
(78).  Bussing  and  Holloway  (79)  have  recently  Investigated  obtaining 
depth  profiles  for  two-component  systems,  using  a  Laplace  transform 
formulation.  They  also  present  a  useful  discussion  of  the  hazards  of 
Obtaining  depth  profiles  by  an  inverse  Laplace  transform  of  VAXPS 
data.  Hazell  et  al.  (80)  have  expanded  upon  the  approach  of  Pijolat 
and  Hoiunger  in  their  development  of  a  model  assuming  a  multilayer 
surface  structure.  Their  approach  also  allows  for  an  Interactive 
determination  of  the  Inelastic  mean  free  path  (IMFP,.  Hecently,  Yih 
and  Hatner  (8,,  have  taken  a  different  approach  to  obtaining  the  depth 
profile  from  VAXPS  data.  Their  method,  based  on  inverse  theory. 
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involves  the  questionable  assumption  that  a  delta  function  can  be 
approximated  by  a  linear  combination  of  a  small  number  of  exponential 
decay  functions.  One  result  of  this  assumption  is  a  smoothing  of  the 
depth  profile,  making  the  method  unsuitable  for  the  determination  of 
depth  profiles  «ith  sharp  transitions.  All  these  models  suffer  from 
instabilities  of  one  kind  or  another  in  the  calculation  of  the  depth 
profile.  In  the  depth  profiling  method  presented  in  this  chapter, 
several  steps  are  taken  to  minimize  the  potential  instability. 

Several  other  methods  exist  for  obtaining  depth  profiles  of  the 
near-surface  region  of  a  sample.  Many  such  methods  make  use  of  Ar^ 
ion  etching,  in  conjunction  with  either  XPS.  AES.  secondary  ion  mass 
spectroscopy  (SIMS),  or  ion  scattering  spectroscopy  (ISS).  However, 
typical  sputtering  rates  of  approximately  10  nm/min  limit  the  depth 
resolution.  Preferential  sputtering  of  species  in  the  sample  matrix 
can  be  difficult  to  take  into  account  iS2).     Sputtering  also  has  an 
irreversible  effect  on  the  surface  topography  (83.8^).  Electron 
energy  loss  spectroscopy  has  been  used  to  obtain  qualitative  profile 

information  (85):  howevpr  i-hia  i-^^i, i 

oj.   however, ^thi3  technology  is  still  developing,  and  it 

involves  the  use  of  a  surface-disruptive  electron  beam. 
Nondestructive  methods  thus  possess  the  advantages  of  higher  depth 
resolution,  preservation  of  the  integrity  of  the  sample  surface,  and 
avoidance  of  the  problem  of  preferential  sputtering. 

Numerical  Inverse  r.aplace  Transform  Method  DevelonmPnf. 
Theoretical  Basis 

Experimentally,  the  VAXPS  method  involves  a  variation  in  the 
angle  of  photoelectron  escape  9  and  the  collection  of  peak  intensities 
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122 


(or  areas)  as  a  function  of  this  angle.  By  varying  0,  one  may  vary 
the  effective  electron  mean  free  path  and  thus  the  sampling  depth.  A 
greater  value  of  9  results  in  a  larger  sampling  depth.  The  depth 
profile  N^(z)  is  contained  in  this  VAXPS  data;  however,  the  data  must 
be  manipulated  before  the  depth  profile  becomes  apparent.  In  Chapter 
2,  the  equation  of  importance  for  analyzing  VAXPS  data  to  obtain  depth 
profile  information  was  shown  to  be    ■ 

1^(9)  =  K^(e)  /*  N^(z)  exp  C-z/X(E^)-3in  9]  dz. 

The  intensity   (or  peak  area)  is  given  by  an  integral  equation 
involving  the  depth  profile  N^(z)  multiplied  by  an  exponential   term 
reflecting  the  signal  attenuation  due  to  inelastic  scattering.     If  the 
assumption  of  a  constant  value  for  the  IMF?  A(E^)  is  made,   this 
equation  can  be  rewritten  as 

■^A^^^  '  ^•Kj^(3)   Iq  N^(z/X)  exp   C-zs/X]  d(z/\), 

where  s  =  1/sin  9.     The  Laplace  transform  F(s)  of  a  function  f(t)  is 
defined  for  t  >  0  by  the  integral  equation  ;'""v;\ 

■'^•C--.\:  F(3)  -  £[f(t)]  =  /^  f(t)  exp   (-St)  dt.     ■    ^; 

Recognizing  the  analogy  between  (z/X)  and  t,  the  intensity  equation 
can  be  written  as  i   p  f  \  }  r\  P:J^  ».  *  .  ff  ;*1  C^ 


c  ,     ■:•  y^         ^:'rf*'\t 


1,(3)  =  X-K  fs)-£[N  fz/A)]."^   ^' 
A         A       A 


Mathematically,  the  depth  profile  N^(z/X)  as  a  function  of  the  reduced 
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depth  z/X  should  then  be  obtained  by  the  relation 
-:^  /         N^(z/X)  =  £"''[I^(3)/(A.K^(3))]  =  £"^[1^(3)]. 

where  I.Cs)  is  the  normalized  peak  area. 

However,  the  inversion  process  is  inherently  instable,  and 
special  precautions  must  be  taken  to  ensure  that  the  resulting  depth 
profile  is  a  realistic  approximation  to  the  true  species  depth 
profile.  The  basic  reason  for  the  instability  of  the  inverse  Laplace 
transform  is  illustrated  in  Figure  5-1.  Shown  are  the  depth  profiles 
for  three  experimental  power  law  functions  and  their  corresponding 
intensity  curves.  Note  the  similarity  between  the  intensity  curves, 
even  though  the  depth  profiles  are  significantly  different.  A  small 
change  in  the  intensity  function  can  cause  a  resultant  large  change  in 
the  depth  profile. 

The  experimental  data  consist  of  a  series  of  peak  area  values  as 
a  function  of  the  instrumental  escape  angle  setting.  These  data  are 
input  into  the  model;  the  depth  profile  for  the  selected  species  is 
output  as  a  function  of  the  reduced  depth  (z/X).  The  main  reason  for 
reporting  depth  profiles  as  a  function  of  (z/X)  lies  in  the  relative 
immaturity  of  the  field  of  IMFP  determination.  For  all  but  the 
simplest  sample  matrices,  the  IMFP  cannot  be  calculated  with  great 
accuracy.  Thus,  to  avoid  compromising  the  utility  of  the  resultant 
depth  profiles,  they  are  reported  as  a  function  of  the  reduced  depth. 
Mathematical  Model 

The  numerical  inverse  Laplace  transform  (NILT)  model  developed 
here  consists  of  two  levels.  In  order  to  minimize  the  appearance  of 


1         '  !      T 


A 


Figure  5-1.   Comparison  of  three  simulated  intensity  functions  as  a 
function  of  s  »  1/sin  9,  along  with  their  corresponding 
inverse  Laplace  transforms.  The  vertical  axes  are  in 
arbitrary  units. 
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instabilities  in  the  output  depth  profile,  an  initial  approximate  fit 
to  the  data  is  necessary.  This  first  level  consists  of  a  number  of 
potential  analytic  functions  (Table  5-1  )  whose  Laplace  transforms  are 
known.  Thus,  the  form  of  the  intensity  as  a  function  of  s  =  I/sin  9 
can  be  matched  to  the  analytic  function  which  gives  the  best  least 
squares  fit.  For  example,  if  a  particular  species  is  thought  to 
possess  a  step  function  depth  profile  such  that 


NAz/X) 
A 


i  n 


if  0  <  z/X  <  k 


if  z/X  >   k 


then  the  intensity  as  a  function  of  s  should  follow  the  relationship 


I^(s)  =  C  exp  (-ks)  s 


■1 


^■-^i 


By  matching  the  actual  intensity  curve  to  its  expected  form,  best-fit 
values  for  C  and  k  can  be  obtained.  In  this  manner,  a  suitable  first 
approximation  to  the  depth  profile  can  be  obtained. 

At  this  point,  several  options  are  available.  If  a  rough, 
qualitative  idea  of  the  depth  profile  is  all  that  is  desired,  then  the 
program  output  at  this  point  typically  suffices.  However,  to  gain  a 
true  understanding  of  the  species  depth  profile  behavior,  a  more 
detailed  analysis  is  necessary.  This  second  level  makes  use  of  a 
numerical  inverse  Laplace  transform  technique,  where  fewer  constraints 
as  to  the  actual  output  depth  profile  exist.  The  variables  are  first 
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Analytic  Depth  Profile  Functions  Available  in 


the  NILT  Program 


'V: 


step  Function 
Linear  Increase  with  Cutoff  Depth 
^*'         Power  Law  Function      "' 
Power  Law  Increase  with  Cutoff  Depth 
Overlayer  Distribution  .jt''..  ; 
•  '^-^■'■'~   Exponential  Function 
Linear  Function  with  Initial  Concentration 
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transformed  to  give  the  definite  integral         ,  v 

"-  '■  -i-  '  ;"  • 

■  /•  ^A^^^  °  -fo  N^^-log  z/O  {z/\f~^   d(z/X), 

which  is  subjected  to  Gaussian  quadrature.  If  n  experimental  points 

are  taken,  then  a  system  of  linear  equations  can  be  set  up,  such  that 

::  ^  ■   •    'a^  V-'  /  V?  »v*  ■>-.'  »  <  • 

n        s  -1 
I^(Sj)  =  I  w^(z^/A)  J   N^(-log  z^/X),     j=1,2,...,n. 

Here,  the  w.   and  z./X  are  tabulated  values  of  the  weights  and  roots, 
respectively,  of  the  shifted  Legendre  polynomial  of  degree  n.     Thus, 
the  depth  profile  is  available  as  a  vector  of  values  of  N  (-log  z  /\) 
by  solving  this  system  of  n  linear  equations.     However,  a  straight- 
forward matrix  solution  approach  is  subject  to  instabilities,  as 
illustrated  in  Figure  5-2.     As  the  number  of  points  n  is  increased—a 
procedure  necessary  for  achieving  higher  depth  resolution— the 
resultant  depth  profile  deviates  more  and  more  from  the  true  output 
curve.     Thus,  some  means  of  damping  the  instabilities  must  be 
implemented  in  order  to  obtain  realistic  depth  profiles.     A  dynamic 
programming  approach   (86)  has  been  used  to  obtain  high  depth 
resolution  and  a  reasonably  stable  output  vector.     A  flow  chart 
illustrating  the  basic  steps  in  determining  species  depth  profiles  is 
given  in  Figure  5-3.     By  initially  obtaining  stable,  low  resolution 
approximations  to  the  profile  and  gradually  approaching  a  higher 
resolution  solution  iteratively,   the  problems  of  instability  are 
minimized.     The   basic  premise  of  the  dynamic  programming  algorithm, 
which  was  developed  to  solve  ill-conditioned  matrix  systems,  involves 
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Figure  5-3- 


Flow  chart  illustrating  the  basic  approach  of  the 
dynamic  programming  routine. 
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the  minimization  of  the  expressi 


on 


(Ax  -  b,  Ax  -  b)  *  r  (X  -  c.  X  -  c), 


where 


r'' 


'  f  '  i.         J 


(x.y) 


•  '   i   t   ~  1   "t"  I  ■» "  \  I  •  ^ ■ 


^    Vi* 


,^  J.  /-  ■ .'  /%-r-t 


i     f  ,i'^ 


H«re,  A  represents  the  ooemolent  matrix  In'the  lysVe.  of  ainelr 
eguatlons.  x  represents  the  depth  profile  vector,  b  Is  the  vector  of 
experimental  Intensity  values,  c  Is  an  Initial  guess  to  the  depth 
profile  vector,  and  r  Is  a  mixing  parameter  In  the  range  [o.n.  The 
-aller  r  Is,  the  faster  Is  the  approach  to  a  convergent  solution, 
however,  the  greater  are  the  resulting  solution  vector 
instabilities.  An  Iterative  approach  to  finding  the  depth  profile  Is 
used,  the  vector  c  represents  the  approximation  to  the  profile  vector 
fro™  the  previous  Iteration.  By  mixing  in  a  fraction  of  this  second 
term,  provision  is  made  for  allowing  the  depth  profile  vector  to  be 
improved  while  keeping  instabilities  to  a  minimum.  Although  the 
method  was  developed  for  use  with  absolute  normalized  pea.  areas, 
.epth  profile  ratios  may  be  used  to  eliminate  the  need  for  calculation 
^r  the  K,(s)  term.  Thus,  a  potential  source  of  error  may  be  removed, 
improving  the  reliability  of  the  output  data  vector. 

Testing  with  Theoretical  Functions  ■  ii:  ; 
in  order  to  determine  the  reliability  of  the  NILT  method,  several 
known  theoretical  functions  were  used  to  simulate  VAXP3  data.  The 
program  output,  after  ,00  Iterations  of  the  dynamic  programming 
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routine,  is  shown  for  these  functions  in  Figure  5-^  (Linear  Increase 
With  cutoff).  Figure  5-5  (Overlayer  Distribution).  Figure  5-6  (Power 
Law  Function),  and  Figure  5-7  (Exponential  Increase).  The  intensity 
distributions  are  shown  in  (a)  and  the  depth  profiles  in  (b).  The 
curves  in  these  figures  represent  expected^ depth^ prof iles.  while  the 
symbols  correspond  to  program  output.  ■^■^'*^^^^ 

Several  conclusions  about  the  nature  of  the  dynamic  programming 
solution  vector  may  be  drawn  from  an  ^alysis  of  these  simulated  data 
results,  in  general,  the  agreement  is  quite  good,  although  the 
program  does  not  account  well  for  sharp  transition  regions.  This 
apparent  weakness  in  the  model  results  from  its  original  design 
criteria.  The  output  vector  was  assumed  to  be  relatively  smooth; 
thus,  any  possible  large  changes  in  the  vector  could  be  attributed  to 
the  inherent  instability  of  the  inverse  Laplace  transform  process  and 
thus  could  be  ruled  out.  If  enough  outside  information  about  a 
particular  surface  system  is  available  to  conclude  that  a  sharp 
transition  in  species  concentration  occurs  in  the  near-surface  region. 
then  one  of  the  several  analytic  functions  in  the  program  which  assume 
a  Sharp  cutoff  region  may  be  used  to  achieve  the  best  fit  to  the  VAXPS 
data.  In  addition,  the  dynamic  programming  routine  can  be  made  to 
approach  the  solution  vector  more  slowly  by  increasing  the  value  of 
the  mixing  parameter  r.  at  the  expense  of  longer  computation  time. 
This  Choice  allows  smaller  changes  to  be  made  in  the  solution  vector; 
however,  the  same  instabilities  will  manifest  themselves,  given  a 
sufficient  number  of  iterations.         ' 


Figure  5-^.       Test  results  for  the  NILT  program.  Linear  increasing 
function  with  cutoff  value  for  z/X:  a)  the  simulated 
normalized  VAXPS  intensity  curve,  and  b)  the 
corresponding  normalized  depth  profile.  The  curve  shows 
the  predicted  response,  while  the  program  output  (+)  is 
shown  in  comparison. 
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Figure  5-5.   Test  results  for  the  NILT  program.  Over layer 

function:  a)  the  simulated  normalized  VAXPS  intensity 
curve,  and  b)  the  corresponding  normalized  depth 
profile.  The  curve  shows  the  predicted  response,  while 
the  program  output  (+)  is  shown  in  comparison. 
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Figure  5-6.   Test  results  for  the  NILT  program.  Po^^^^  -^^^ 
**  distribution:   a)  the  simulated  normalized  VAXPS 

■'-■     intensity  curve,  and  b)  the  corresponding  normalized 

depth  profile.  The  curve  shows  the  predicted  response, 
while  the  program  output  (+)  is  shown  in  comparison. 
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Figure  5-7.   Test  results  for  the  NILT  program.  Exponential 

increase:  a)  the  simulated  normalized  VAXPS  intensity 
curve,  and  b)  the  corresponding  normalized  depth 
profile.  The  curve  shows  the  predicted  response,  while 
the  program  output  (+)  is  shown  in  comparison. 
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Experimental  Testing  of  the  NILT  Method 
The  results  of  the  NILT  method  have  been  reported  for  an  oxidized 
aluminum  foil  sample  (63).  Two  other  examples  of  its  application  to 
the  determination  of  uncomplicated  vertical  inhomogeneity  in  the  near- 
surface  region  are  presented  here.  The  samples  discussed  are  the  gold 
foil  and  oxidized  carbon  foil  samples  whose  data  collection  was 
outlined  in  Chapter  4. 

The  gold  foil  VAXPS  response  curves  for  the  Au  If  peaks  using 
both  Al  Ka  and  Mg  Ka  radiation  were  illustrated  in  Figure  4-12.  These 
data  were  subjected  to  the  NILT  method.  Initially,  a  power  law 
functional  form  for  the  depth  profile  was  assumed,  to  analytically 
obtain  a  reasonable  initial  guess  to  the  actual  depth  profile.  These 
initial  approximations  were  next  input  into  the  dynamic  programming 
routine,  generating  the  output  depth  profiles  shown  in  Figure  5-8. 
Several  conclusions  can  be  drawn  from  an  examination  of  this  plot. 
First,  an  increase  in  the  gold  concentration  of  approximately  an  order 
of  magnitude  over  the  range  shown  indicates  that  the  ubiquitous 
hydrocarbon  contamination  is  localized  within  the  near-surface 
region.  Of  more  general  interest  is  the  observation  that  both  choices 
of  exciting  radiation  give  similar  resulting  depth  profiles.  Thus, 
either  type  of  exciting  radiation  may  be  used  to  carry  out  VAXPS 
experiments  for  nondestructive  depth  profiling  applications.  Finally, 
this  analysis  was  carried  out  without  the  use  of  ratios. 

An  oxidized  carbon  foil  sample  was  studied  to  determine  the 
nature  of  the  oxidation  layer  in  the  near-surface  region.  The  signal 
was  corrected  for  instrumental  effects,  and  the  graphitic  carbon  peak 
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area  was  separated  from  the  oxidized  carbon  peak  area  using  a  peak 
fitting  routine.     The  resulting  VAXPS  response  curves  are  shown  in 
Figure  5-9  for  the  overall  peak  area,   graphitic  carbon  peak  area,  and 
oxidized  carbon  peak  area.     The  latter  two  curves  were  then  input  into 
the  NILT  program  in  the  manner  described  above.     A  ratio  of  the 
graphitic-to-oxidized  carbon  depth  profile  is  shown  in  Figure  5-10. 
As  expected,  the  curve  indicates  that  oxidation  is  confined  to  the 
near-surface  region.     The  values  for  the  ratio  are  realistic,  due  to 
the  cancellation  of  the  K^(s)  term. 
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CHAPTER  6 

QUANTITATIVE   XPS  OF 

LEAD  ION-SELECTIVE  ELECTRODE  MEMBRANES 


Introduction 
The  reaearch  described  up  to  this  point  In  this  dissertation  has 
had  a  two-fold  purpose.     First,   It  has  been  carried  out  to 
charaoterue  several  general  problems-the  effects  of  surface 
topography  and  a  „eans  of  nondestructive  depth  proflung-which  have 
not  received  .uch  Interest,   but  which  are  of  obvious  Importance  to  the 
development  of  quantitative  XPS.     However,   the  driving  force  of  this 
wor.  since  Its  Inception  has  been  related  to  how  these  areas  can  be 
applied  to  the  study  of  lead  ISE  membrane  surfaces.     Even  though  the 
methods  developed  to  tac.le  these  problems  have  been  general  In  scope, 
their  effectiveness  in  the  study  of  lead  ISE  membranes  has  been  of 
prime  concern.  -  .-.  ■•.. 

Specifically,  the  experimental  surface  roughness  work  described 
in  a  previous  chapter  demonstrated  that  a  smooth  surface  could  be 
formed  from  a  powder  by  pressing  it  in  a  die.  provided  that  the 
pressure  was  sufficiently  high.     In  this  fashion,   the  surface 
topographical  effects  demonstrated  could  be  neglected.     Thus,   by 
pressing  the  lead  ISE  powders  at  sufficiently  high  pressures,   any 
surface  topographical  effects  are  expected  to  be  negligible. 

The  application  of  the  numerical   inverse  Laplace  transform   (NILT) 
method  to  VAXPS  data  collected  from  various  Pb  ISE  membranes  is  the 
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focus  of  the  present  chapter.  First,  the  data  treatment  necessary  to 
prepare  the  VAXPS  data  for  NILT  analysis  is  discussed.  Then  three 
separate  experimental  studies  of  Pb  ISE  membranes  are  presented.  The 
effects  of  membrane  composition  and  preparation  on  the  distribution  of 
chemical  species  in  the  near-surface  region  are  investigated.  Also,  a 
chemical  means  of  electrode  surface  cleaning  is  examined.  Finally, 
the  nature  of  the  poisoning  effect  of  Fe^"^  ions  on  the  Pb  ISE  signal 
is  elucidated.  '"■■ 

Experimental  Analysis  of  Pb  ISE  Membranes 
Operating  Conditions  and  Preliminary  Data  Treatment 

The  study  of  Pb  ISE  membranes  outlined  here  involves  an 
examination  of  the  effects  of  composition,  preparation,  cleaning,  and 
poisoning.  The  membranes  consist  of  a  co-precipitate  of  PbS  and 
AggS.  The  PbS  serves  as  a  site  at  the  surface  where  Pb^"^  ions  in 
solution  can  participate  in  the  solubility  equilibrium,  which  is 
related  to  the  electrode  potential.  However,  PbS  is  a  narrow  band-gap 
semiconductor,  an  ill-suited  choice  for  an  electrode  material.  The 
Ag2S,  which  is  an  excellent  ionic  conductor,  thus  serves  to  allow  the 
[Pb  ]  information  to  be  couched  in  terms  of  a  potential— which  can  be 
measured  and  related  back  to  [Pb^^D—across  the  ISE  membrane  (16). 
This  discussion  suggests  several  crucial  requirements  which  must  be 

met  in  order  to  have  a  stable,  responsive  Pb  ISE  of  this  type.  First, 

2+ 
Pb   sites  must  be  present  at  the  surface  in  order  to  allow  an 

exchange  equilibrium  to  be  set  up.  Secondly,  these  sites  must  not  be 

oxidized.  Finally,  there  must  be  sufficient  Ag  S  throughout  the 

membrane  to  provide  a  conductive  pathway. 
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The  operating  conditions  and  data  pretreatment  for  the  Pb  ISE 
membrane  samples  studied  were  uniform.  The  XSAM  800  was  run  in  fixed 
retarding  ratio  (FRR)  mode,  using  the  Al  Ka  excitation  source 
operating  at  15  kV  accelerating  potential  and  15  mA  filament  emission 
current.  The  sample  powder  was  pressed  for  1  minute  at  25000  psi  in 
an  evacuated  die  yielding  pellets  16  mm  in  diameter;  the  mass  of 
powder  used  in  each  case  was  0.i45  g.  The  sample  was  cut  out  of  the 
pellet  to  fit  the  sample  holder,  to  which  it  was  affixed  using 
conductive  colloidal  silver.  The  spectra  were  collected  at  high 
resolution  and  high  magnification  for  six  values  of  the  instrumental 
take-off  angle:  70o.  25°.  20o.  150.  iQo.  and  50.   It  was  assumed  that 
the  only  elements  present  in  the  near-surface  region  were  oxygen. 
Silver,  carbon,  sulfur,  and  lead;  this  assumption  was  verified  using 
survey  scans  of  the  samples,  an  example  of  which  is  shown  in  Figure 
6-1.  Note  that  several  of  the  peaks  are  difficult  to  discern.  Their 
peak  Shape  and  area  can  be  determined,  however,  given  sufficient  data 
collection  time.  For  each  sample  at  each  take-off  angle,  several  high 
resolution  energy  windows  were  collected,  for  the  0  Is.  Ag  3d,  C  Is, 
S  2s,  and  Pb  Hf  regions.  The  peak  areas  were  corrected  using  sample- 
specific  instrumental  response  factors,  peak  photoionization  cross- 
sections  (31),  angular  distribution  asymmetry  parameters  (35),  and  the 
assumption  of  an  IMFP  energy  dependence  of  the  form  A  =  K  •  E^'^  (39). 
Application  of  the  NILT  depth  profiling  method  was  carried  out  in 
a  systematic  fashion  for  each  sample.  Of  chief  interest  in  these 
studies  were  the  variations  in  the  ratios  of  Ag^  to  ?t'\   PbSO,  to 
PbS,  and  total  sulfate  to  total  sulfide  as  a  function  of  depth.  In 
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addition,  the  total  carbon  depth  profile  was  investigated.  These 
depth  profiles  show  the  degree  of  oxidation,  as  well  as  the 
availability  of  Pb^*  for  exchange  at  the  surface.  All  the  requisite 
peaks  were  normalized  ir.  the  above  fashion,  then  input  into  the  depth 
profiling  program.  A  power  law  function  was  uniformly  used  to 
generate  the  initial  approximation  to  the  solution  vector,  which  was 
subsequently  fed  into  the  dynamic  programming  routine. 
Composition  and  Preparation  Effects 

The  NILT  method  was  first  used  to  test  whether  it  could  discern 
the  ratios  of  Ag*  and  Pb^*  as  a  function  of  depth  for  three  samples 
differing  in  composition.  The  three  samples  contained  a  mole  ratio 
for  PbS:Ag2S  of  3:7.  1:1.  and  7:3.  The  depth  profile  ratios  are 
compared  in  Figure  6-2.   They  show  the  expected  order  for  the  AgVPb2- 
ratio,  in  addition,  all  three  curves  share  a  similar  shape,  with  no 
great  changes  in  the  concentration  ratio  over  the  range  of  reduced 
depth  values.  These  powders  were  all  prepared  in  the  same  fashion,  by 
adding  the  salts  in  a  dropwise  fashion  to  the  sulfide  solution.  Thus, 
the  continual  presence  of  excess  sulfide  prevented  any  competition 
between  formation  of  PbS  (K^^   =  3.^*  x  lO'^S  ^t  l8oc)  and  Ag^S 
%  =  1.6  X  10-^9  ^,  ,g„^^  ^^^^  occurring.  A  comparatively  flat 
profile  is  expected,  with  no  general  increasing  or  decreasing  trend. 
The  curves  reflect  the  expected  behavior.  An  interesting  contrast  may 
be  Observed  by  comparing  the  AgVPb^*  depth  profiles  for  1:1  mole 
ratio  powders  prepared  either  in  the  above  manner  (procedure  2. 
discussed  m  Chapter  .).  or  by  adding  sulfide  to  the  salt  solution 
(procedure  1).  since  the  latter  method  ensures  that  the  sulfide  will 


-,:'v^"--".V  -  .  ■"  ■■■■'  "■■ 


m 
o 

-.J: 

•iH 

4J 

OS 

L 

<0 

iH 

i 

- .  _.  .1      ■ 

5 

'.;  ;  .■;. 

•^ 

» 

a 

t. 

■  ^-  '■  • ,  ■     , 

o 

•o 

-            ',     '       ."              '■'':■.       ■' 

•- 

0 

'■       '\'       ■( 

o. 

m 

<NJ 

bO 

< 

-                                        "             i'''- 

•• 

CO 

o 

ou 

£. 

>:.• 

o 

4^ 

VO  f*f'^ 

n 
o 

'i\':V-V. 

•H 

^ 

2 

+ 

CM 

£ 

*.  ■■   ■f* 

f«fc,-*-W  V. 


5 


o  •• 
s-  c- 
a 

£  c 

o. 

<D      « 

+ 

<i-.  .- 
O 

C   '^ 

O    X 

n  ^ 

nj    .. 
a.  t^ 

E 
O  *-i 

o  o 


I 
SO 


158 


WC  -v^ 


'■.;'  «"i  f^   ?  -.  ?  ^  a  /^ .'  *■  I   i 


^r"*. 


j'lt  .», . .  ^ 


;^;v 


-■  U  iJ  W  □ 


X 


X 


X       + 


'& 


1^ 


n 


(M 


■■CO 


t    N 


Q 


oiiyy  ad/9y 


be  the  limiting  reagent  in  the  reaction,  the  leas  soluble  AgoS  would 
be  expected  to  form  preferentially.  As  indicated  in  Figure  6-3,  this 
is  indeed  the  case.  The  Ag  /Pb   ratio  clearly  increases  with 
increasing  depth,  lending  credibility  to  this  line  of  reasoning. 
Indeed,  this  method  of  preparation  can  lead  to  powders  with  a  very 
small  amount  of  Ag2S  at  the  surface.  Electrodes  manufactured  from 
such  powders  cannot  function  properly,  due  to  the  lack  of  the  normal 
ionic  conductive  pathway  across  the  ISE  membrane.  Also  of  importance 
to  the  manufacture  of  suitable  powders  for  constructing  lead  ISEs, 
procedure  2  seems  to  yield  surfaces  with  a  significantly  lower  degree 
of  oxidation,  which  can  be  measured  by  the  relative  peak  areas  of 
sulfate  and  sulfide  species.  These  species  give  clearly  separated 
peaks  in  the  S  2s  spectrum  due  to  the  large  difference  in  the  valence 

of  sulfur,  as  illustrated  in  Figure  6-4.  The  depth  profile  plot  of 

2-  2- 
the  SO^  /S   ratio  for  these  two  species  is  shown  in  Figure  6-5.  A 

threefold  difference  in  the  degree  of  oxidation  is  apparent.  The 

difference  is  even  more  dramatic  if  one  considers  the  ratio  of  PbSOi, 

to  PbS.  The  use  of  XPS  allows  these  species  to  be  resolved,  as  shown 

in  Figure  6-6.  The  Pb  ^f^^g  °an  be  resolved  into  three  components, 

present  at  137.6  eV,  139.0  eV,  and  139.6  eV.  The  first  peak  can  be 

assigned  to  PbS  and  the  third  peak  to  PbSOj^  (27,87).  There  is  still 

some  uncertainty  as  to  the  origin  of  the  second  peak.  It  is  in  the 

general  region  where  a  basic  lead  carbonate  peak  would  be  located 

(88),  and  the  presence  of  a  carbonate  peak  at  289.0  eV  in  the  C  Is 

spectrum,  shown  in  Figure  6-7,  supports  this  hypothesis  (27).  The 
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Figure  6-6. 


Illustration  of  the  peak  fit  procedure  yielding  three 
peak  doublets  in  the  Pb  4f  spectrum  of  an  oxidized  Pb 
ISE.  The  ^fy/p  peaks  correspond  to  lead  sulfide 
(137.6  eV),  a  lead  carbonate  species  (139.1  eV),  and 
lead  sulfate  (139.5  eV) .  The  top  curve  shows  the 
experimental  spectrum  of  intensity  versus  binding 
energy,  while  the  bottom  curve  shows  the  matching 
simulated  spectrum. 
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Fteure  6-7.   Illustration  of  the  peak  fit  procedure  on  the  C  is 
rlgure   7.   ^^^^^   indicating  the  presence  of  a  carbonate  species 
at  ?8q  0  eV.  The  top  curve  shows  the  experimental 
spectrum  of  intensity  versus  binding  energy,  while  the 
bottom  curve  shows  the  matching  simulated  spectrum. 
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effect  of  sample  preparation  on  the  relative  amounts  of  PbS  and  PbSOn 
can  be  seen  in  Figure  6-8.  The  sample  prepared  using  procedure  2 
shows  a  sixfold  smaller  degree  of  lead  species  oxidation  at  the 
surface,  with  an  even  greater  difference  in  the  bulk. 

Thus,  addition  of  the  salt  solution  to  excess  sulfide  solution  is 
to  be  preferred  in  the  preparation  of  Pb  ISEs,  since  this  method 
provides  a  greater  proportion  of  unoxidized  species  in  general,  and 
the  lead  ion-sensing  material  PbS  in  particular,  in  the  near-surface 
region.  At  the  same  time,  the  order  of  addition  ensures  that  no 
significant  segregation  of  species  should  occur,  thus  allowing 
suitable  membrane  ionic  conductivity. 
Effects  of  HClOi, 

The  use  of  HClOn  solutions  to  act  as  cleaning  agents  for  lead 
ISEs  has  been  reported  previously  (24,27,89).  However,  the  effects  on 
the  near-surface  region  were  not  clearly  defined.  A  VAXPS  study  was 
carried  out  to  determine  the  action  of  such  a  solution  on  the  surface 
composition.  The  powder  was  prepared  using  procedure  2  discussed 
above.  It  was  then  allowed  to  sit  in  the  open  atmosphere  for  several 
months,  to  ensure  its  need  for  cleaning.  It  was  analyzed  using  VAXPS, 
then  exposed  to  a  0.1  M  HCIO^  solution  for  15  minutes,  rinsed  in  a 
stream  of  deionized  water,  and  blotted  dry.  The  sample  was  analyzed 
using  VAXPS  after  this  treatment  as  well.  The  two  experimental 
results  are  compared  in  Figure  6-9  for  the  three  depth  profile  plots 
most  likely  to  show  aging  affects:  SOJj^/s"^  ratio,  PbSO^/PbS  ratio, 
and  C  Is  total  carbon.  The  plots  uniformly  show  significant 
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Figure  6-9.       Comparison  of  the  depth  profiles  of  the  SO^  /S       and 

PbS02|/PbS  ratios,  and  normalized  C  Is  total  carbon  for  a 

Pb  ISE  membrane  before   (+)  and  after    (x)  exposure  to 

f|       HClOjj  solution. 
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differences  in  the  profiles  for  z/\  <   -1 ,  with  the  cleaning  action 
being  predominantly  localized  for  z/A  <  -0.5.  The  near-surface 
oxidized  species  have  been  cleaned  off,  as  has  the  hydrocarbon 
contamination  in  this  region.  Thus,  HClOi,  serves  the  purpose  of 
reducing  the  degree  of  surface  oxidation,  while  decreasing  the 
hydrocarbon  contamination  as  well.  However,  the  need  for  extended 
exposure  times  is  indicated  by  the  limited  amount  of  penetration 
achieved  by  the  acid,  as  illustrated  by  the  VAXPS  results. 
Effects  of  Fe^"^ 


•  1 


Several  authors  have  previously  reported  the  irreversible  changes 
in  the  sensitivity  of  lead  ISEs  caused  by  exposure  to  Fe^"^ 
(15,16,21,22).  However,  no  research  into  the  changes  in  surface 
composition  caused  by  this  ion  has  been  reported.  Thus,  a  study  was 
carried  out  to  define  the  nature  of  the  Fe^*  effect  on  the  lead  ISE 
response.  A  membrane  with  a  1:1  mole  ratio  of  PbStAgpS  was  prepared 
using  procedure  2  and  analyzed  using  VAXPS.   It  was  then  exposed  for 
15  minutes  to  a  0.1  M  Fe(N02)3  solution,  rinsed  with  deionized  water, 
and  blotted  dry.  The  pellet  was  again  analyzed  by  VAXPS.  The  chief 
difference  in  the  sample  composition  is  shown  in  Figure  6-10.  The 
Fe   acts  to  reduce  the  lead  concentration  throughout  the  near-surface 
region  by  approximately  an  order  of  magnitude,  leaving  a  large  excess 

of  silver  species  at  the  surface.  Essentially  all  the  exchange  sites 

2+ 
for  Pb   are  removed,  thus  poisoning  the  electrode  response.  Note 

also  that  the  effect  is  manifested  at  much  greater  depths  than  the 

cleaning  action  of  HCIO^.  No  evidence  of  the  presence  of  iron  species 
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at  the  surface  of  the  treated  membrane  could  be  found;  thus,  the  Fe^* 
does  not  exchange  with  the  Pb^*,  as  has  been  observed  with  Cu^*  ions, 
for  example  (25).  The  depth  profile  for  the  SoJ~/S^"  ratio  (Figure 
6-11)  shows  a  major  change  as  well.  This  graph  provides  convincing 
evidence  that  the  poisoning  action  of  Fe^"*"  ions  is  not  due  solely  to 
surface  oxidation.  This  ratio  decreases  by  at  least  a  factor  of  five, 
due  chiefly  to  the  reduction  in  the  lead  content,  since  PbSOj^  is  the 
predominant  sulfate  species.  The  lead  species  which  remain  show  a 
similar  degree  of  oxidation  near  the  surface,  as  shown  in  Figure 
6-12.  However,  there  is  an  indication  that  the  degree  of  oxidation 
increases  with  increasing  depth.  It  is  possible  that  the  ions 
penetrate  into  the  sample  and  act  as  oxidizing  agents  below  the 
surface,  although  the  evidence  indicated  is  not  conclusive  enough  to 
make  this  statement  with  certainty. 

To  summarize  the  results  of  the  study  of  PblSEs,  it  was  found 
that  a  variation  in  composition  could  be  detected  by  VAXPS,  but  the 
ratios  studied  did  not  have  a  substantial  effect  on  the  distribution 
of  species  near  the  surface.  In  contrast,  the  method  of  powder 
preparation  was  found  to  yield  substantial  differences  in  the  degree 

of  surface  oxidation,  and  a  rather  less  dramatic  difference  in  the 

+   2+ 
variation  of  Ag  /Pb   with  depth,  the  latter  of  which  may  be  explained 

by  the  competition  for  sulfide  ion  when  it  is  the  limiting  reagent. 

It  is  recommended  that  of  the  two  procedures  studied,  addition  of  the 

salt  solution  to  the  sulfide  solution  is  preferable  because  of  the 

resultant  lower  degree  of  oxidation,  as  well  as  the  homogeneity  of  the 
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+   2+ 
Ag  /Pb   ratio.  The  cleaning  effect  of  HCIO^  was  proven  to  be 

localized  in  the  near-surface  region;  the  action  of  this  acid  served 

to  reduce  the  proportion  of  oxidized  species,  as  well  as  the 

hydrocarbon  overlayer  content.  Finally,  the  irreversible  degradation 

of  Pb  ISE  response  caused  by  exposure  to  Fe^*  was  shown  to  result  from 

a  depletion  through  a  large  region  of  the  near  surface  in  the  amount 

2  + 
of  Pb  ,  with  the  possibility  of  sub-surface  oxidation  processes  being 

indicated.        ,  -  ,  , 
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CHAPTER  7 
CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 


Summary  and  Conclusions 
Three  major  research  projects  have  been  described  in  previous 
chapters.  First,  a  study  of  the  potential  effects  of  surface 
roughness  on  the  nature  and  magnitude  of  the  signal  in  variable  angle 
X-ray  photoelectron  spectroscopy  (VAXPS)  has  been  carried  out  from 
both  a  theoretical  and  an  experimental  standpoint.  The  theoretical 
aspect  involved  the  development  and  testing  of  a  program  designed  to 
output  several  factors  characterizing  the  surface  topographical 
effect,  as  well  as  a  theoretical  VAXPS  curve  modelling  the  electron 
population  as  affected  by  the  surface  roughness.  Testing  was  carried 
out  using  one-dimensional  sinusoidal  and  grating  functions,  each  of 
whose  results  could  be  compared  to  values  arrived  at  analytically. 
The  program  was  initially  applied  to  determine  whether  lead  ion- 
selective  electrode  (ISE)  membrane  powders  prepared  in  different 
fashions  would  show  different  surface  roughness  effects.  A 
particulate  surface  model  (PSM)  was  developed  to  transform  particle 
size  data  obtained  from  scanning  electron  microscopy  (SEM)  into  a 
suitable  input  for  the  program.  The  results  demonstrate  that, 
although  the  SEM  photomicrographs  indicate  noticeably  different 
particle  size  distributions,  the  differences  are  insufficient  to 
produce  a  change  in  the  surface  roughness  effect.  Experimentally,  a 
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number  of  studies  have  been  carried  out.  The  roughness  effect  was 
demonstrated  using  a  series  of  gold  samples  of  varying  topography. 
Next,  the  predictions  of  the  model  as  to  the  effect  of  surface 
roughness  on  the  VAXPS  response  curve  were  lent  strong  credence  by  a 
comparison  of  theoretical  results  and  experimental  VAXPS  data  from  a 
series  of  carbon  samples  differing  in  topography.  It  was  then  shown 
that  the  use  of  a  sufficiently  high  pressure,  greater  than  1000  psi, 
on  a  carbon  powder  specimen  could  yield  a  pellet  with  a  VAXPS  response 
curve  virtually  indistinguishable  from  that  of  a  smooth  carbon  foil. 
Finally,  several  peak  structure  studies  were  carried  out  for  a 
unidirectionally  abraded,  oxidized  carbon  foil.  The  variation  in  the 
relative  contributions  of  surface  (oxidized)  and  bulk  (graphitic) 
carbon  to  the  overall  signal  was  studied  as  a  function  of  surface 
roughness  and  photoelectron  take-off  angle.  Again,  the  results 
provide  supporting  evidence  for  the  validity  of  the  surface  roughness 
model.        v->^'  i^  ii>....v*  ^^.'  ,e-  '■  t  i^  ^j.  .,  ,^  ^^^,-  .  ;v^__.^.r  ■  %,^,^-  '. 

Secondly,  a  means  of  nondestructively  examining  the  near-surface 
region  for  compositional  changes  as  a  function  of  depth  has  been 
developed.  This  method  is  based  on  the  inverse  Laplace  transform 
relationship  between  the  species  depth  profile  and  the  peak  area  for 
that  particular  species.  The  method  attempts  to  minimize  the 
instabilities  inherent  in  performing  a  numerical  inverse  Laplace 
transform.  It  has  been  tested  on  several  simulated  data  vectors  and 
shown  to  yield  satisfactory  results.  It  has  also  been  shown  to  yield 
expected  results  for  the  analysis  of  a  gold  foil  sample  with  a 
hydrocarbon  contamination  layer,  as  well  as  the  analysis  of  an 
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oxidized  carbon  foil  sample.  The  studies  indicate  that  depth  profile 
ratios  should  be  reported  in  all  but  the  simplest  systems.  In 
addition,  the  method  clearly  reflects  any  near-surface  segregation  of 
species  which  might  be  present. 

Finally,  these  methods  were  applied  to  the  study  of  the  surface 
composition  of  lead  ISEs.  The  effects  of  surface  roughness  were 
demonstrated  to  be  minimal,  provided  that  sufficient  pressure  was  used 
m  preparing  the  sample  pellets  for  analysis.  Various  experiments 
were  then  carried  out  using  VAXPS  to  probe  the  near-surface  speciation 
of  these  samples.  It  was  found  that  changes  in  the  relative  amounts 
of  PbS  and  Ag2S  making  up  the  pellet,  although  noticeable,  did  not 
significantly  alter  the  surface  composition.  Method  of  preparation, 
on  the  other  hand,  was  found  to  substantially  affect  the  degree  of 
surface  oxidation,  which  was  lessened  by  using  a  surface  cleaning 
agent  such  as  HC10|,.  The  Interference  of  Fe^*  ions  was  shown  to 
result  from  the  severe  depletion  of  lead  species  throughout  the 
analyzed  depth.  A  suggestion  of  oxidation  processes  taking  place 
below  the  surface  of  iron-treated  Pb  ISEs  was  indicated,  but  the 
experimental  evidence  was  inconclusive.        ■  -. 

-,       Future  Work 
The  surface  roughness  studies  present  a  number  of  possibilities 
for  further  investigation.  Chief  among  these  ideas  is  a  coupling 
between  some  experimental  means  of  surface  topography  determination 
and  the  surface  roughness  program  developed.  Among  current  methods, 
scanning  tunneling  microscopy  offers  the  highest  potential  spatial 
resolution  and  the  greatest  promise.  However,  this  technique  is  still 
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m  its  infancy.  More  mature  methods  such  as  scanning  electron 
microscopy  (90)  and  profilometry  (60)  offer  more  immediate  solutions 
to  the  provision  of  surface  height  matrices.  The  program  itself  may 
serve  as  a  framework  upon  which  more  elaborate  models  may  be  built. 
For  example,  the  effects  of  elastic  scattering-which  have  been 
neglected  in  the  surface  roughness  model  adopted-could  be 
incorporated  into  the  program.  In  addition,  the  ideas  of  fractal 
geometry  as  applied  to  heterogeneous  surfaces  (91)  could  be  used  to 
increase  the  sophistication  of  the  surface  roughness  approximation. 

Future  work  on  the  NILT  method  of  nondestructive  depth  profiling 
Should  be  geared  toward  improving  the  stability  of  the  output  data 
vector.  At  present,  the  method  provides  a  reliable  approximation  to 
the  true  surface  composition  profile.  Research  efforts  to  more 
Clearly  define  the  magnitude  of  the  inelastic  mean  free  path  should 
Obviate  the  need  for  reporting  profiles  as  a  function  of  the  reduced 
depth.  The  potential  applications  of  the  NILT  method  are  limited  only 
by  the  sampling  depth  of  the  photoelectrons.  which  can  be  increased  by 
increasing  their  kinetic  energy,  i.e..  using  more  energetic  sources. 

The  study  of  Pb  ISEs  using  VAXPS  offers  a  number  of  opportunities 
for  research  as  well.  Chief  among  these  potential  studies  is  a   ' 
further  probe  into  the  surface  composition  of  Fe3*-treated  ISEs. 
perhaps  by  a  combination  of  VAXPS  and  Ar*  etching,  in  order  to 
determine  whether  sub-surface  oxidation  is  indeed  occurring.  There 
are  a  large  number. of  preparation  methods  described  in  the  literature 
Which  have  produced  workable  electrodes.  It  would  be  interesting  to 
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investigate  whether  these  preparation  methods  share  any  common 
features  near  the  surface. 

Finally,  a  long-term  research  goal  related  to  the  development  of 
quantitative  XPS  may  be  defined.  In  this  scenario  are  included  both 
the  surface  roughness  and  NILT  models,  along  with  suitable 
instrumental  response  functions,  parameter  definitions,  and  IMFP 
corrections.  Such  a  method,  when  implemented  on  a  state-of-the-art 
XPS  instrument,  would  be  able  to  provide  a  total  quantitative  picture 
of  the  near-surface  region  for  a  sample,  including  segregation 
effects,  correlation  of  peaks,  and  even  lateral  resolution  capability 
using  small-spot  XPS.  Many  of  the  components  of  such  a  system  have 
been  developed,  and  their  association  provides  an  interesting 
challenge  for  the  future. 
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APPENDIX  A 
SURFACE  ROUGHNESS  PROGRAM  SOURCE  CODE 
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c„Io«L^'^°°^^^  RECEIVES  AS  INPUT  A  MATRIX  REPRESENTING  DIGITIZED 

BOTH  THE  X  'amh' ?  nrp.".'"'  '"''  '^''^  ™^  '^^^"^^  RESOLUTiS  IN 
BOTH  THE  X-  AND  Y-DIRECTIONS,  OUTPUTS  THE  SHADING  FRACTION 

SURFACE  AREA.  EFFECTIVE  AVERAGE  THETA'  AND  SIN  ^HE^AM  JaLUES 
AND  RANGE  INFORMATION,  AND  INFORMATION  NECESSARY  TO  CONSTRUCT  AN 
ELECTRON  POPULATION  DISTRIBUTION  FOR  A  GIVEN  SERIES  qf     " 
INSTRUMENTAL  SETTINGS  OF  THE  ANGLE  THETA  OF  PHOTOELECTRON 
tibOAPE, 

o,mS2yP^°^  ^^  *^^°  ^^^^  ^°R  PROGRAM  INPUT  OF  VARIOUS  MODELING 
behavior!  '°  ""'""  ""  TEST  FUNCTIONS  FOR  SURFACE  ROUGHNESS 


DIMENSION  PARDK 70) 

C0MM0N/HARY/HEIGHr(5000,2)/IARY/BL0CK(5000  2)  FRACT 

COMMON/ JARY/PIXAR ( 5000, 1 ) /KARY /DEPAL ( 5000  2 ) 

DATA  PI  /3.m59265il/ 

OPEN  ( UNIT=8, FILE- 'TEST. DAT •,STATUS=' OLD') 

OPEN  (UNIT=12.FILE='DA0UT.DAT',STATUS='NEW') 

FIRST,  WE  READ  IN  THE  INPUT  CONTROL  PARAMETERS. 

READ (8, 100)  ICON,ISAVE  '  - 

F0RMAT(2(5X,Iil))  ,./;':.•; 


IF  ICON 
ICON 
ICON 
ICON 
ICON 
ICON 


1 ,  THEN  INPUT  IS  WITHIN  INPUT  DATA  FILE 

2,  THEN  INPUT  IS  VIA  A  BROWNIAN  FUNCTION 

3,  THEN  INPUT  IS  VIA  A  SINUSOIDAL  FUNCTION 
it,  THEN  INPUT  IS  VIA  A  "RUFFLE"  FUNCTION 

5,  THEN  INPUT  IS  VIA  A  GRATING  FUNCTION 

6,  THEN  INPUT  IS  VIA  A  CLOSEST  PACKING  FUNCTION 


IF  ISAVE  =  0,  THEN  THE  INPUT  DATA  FILE  IS  NOT  SAVED 

ISAVE  =  1,  THEN  THE  INPUT  DATA  FILE  IS  SAVED  TO  UNIT  10. 
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C     NX  AND  NY  ARE  THE  NUMBER  OF  POINTS  IN  X  AND  Y  DIRECTIONS. 
C 

READ(8,100)  NX, NY 
C 

C     PROVISION  IS  MADE  FOR  AUTOCORRELATION  OF  INPUT  VECTOR. 
C 

READ (8, 100)  ICOR.LAG 
C 

C     XTAU  AND  YTAU  ARE  X  AND  Y  RESOLUTIONS,  RESPECTIVELY. 
C 

READ(8, 110)  XTAU, YTAU  .""■-:.■'- 

110  F0RMAT(2(5X,F11.8)) 

WRITE(12,120)  • 

120  FORMAT ('  ',60 (•*')/)   r 

'rfRITE(12,130)  NX, NY   ' '" 
130  FORMATC  ',5X, 'THERE  ARE  ' ,  I4  • ,  AND  \1H, 
&'  POINTS  IN  THE  X  AND  Y  DIRECTIONS.  '/) 
MRITE(12,1ilO)  XTAU, YTAU 
140  FORMAT  ('  ',5X,'X  AND  Y  RESOLUTION  ARE  '.FII.S, 
&'  AND  ',F11.8/) 
C 

C     THE  MATRIX  OF  ROUGHNESS  VALUES  IS  STORED  IN  HEIGHT (I. J) 
C 

IF  (IC0N.EQ.1 )  THEN 

READ(8,150)  ((HEIGHT(I,J),I=1,NX),J=1,NY) 
150    FORMAT (5X,F8. 5) 

WRITE (12, 160) 
160    FORMATC  ',5X, 'INPUT  DATA  IS  WITHIN  THE  DATA  SET.'/) 
ELSE  IF  (IC0N.EQ.2)  THEN   •      - 
C 

C     INPUT  RANDOM  SEED  AND  RANGE  OF  POTENTIAL  VALUES. 
C 

READ(8,100)  IRN,NRANGE 

CALL  BROWN(NX,NY,IRN,NRANGE,ISAVE) 
ELSE  IF  (IC0N.EQ.3)  THEN  ' 

CALL  SINUS(NX,NY,ISAVE)  '"^O'y  .   -  . 

ELSE  IF  (IC0N.EQ.4)  THEN     "*  "' ■  ^     .'  ^V' 

CALL  RUFFLE  (NX ,  NY ,  ISAVE )  :  '.  -  -r:   . 

ELSE  IF  (IC0N.EQ.5)  THEN    I    {   p.  J 1  H      Sv   f 
Q  '  .it  t  '■  *-»■  ■"•««■ 

C     INPUT  PRIMARY  BLAZE  ANGLE. 

READ (8, 180)  RHO  \  "'l'  !  '  "' [  \   .....  . 

1 80    FORMAT  ( 5X ,  F7 . 4  )  '  ■  ■*'  "'-'■  > .'  "  '—  ''  ^ -^    '^   '  — -' 
CALL  GRATE (NX, NY, IS AVE, RHO, XTAU) 
ELSE  IF  (IC0N.EQ.6)  THEN 

OPEN  (UNIT-1 8, FILE- 'PART. DAT' , STATUS- 'OLD' ) 
C 

C     INPUT  RANDOM  SEED  AND  THE  NUMBER  OF  PARTICLES. 
C 


'«» 
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i  i   fV 
READ    (18,100)    IRN.NP       '      '   •     ''*►♦■•-    ■ 
DO  1-1, NP 
C 

C     READ  THE  PARTICLE  DIAMETERS  INTO  PARDI(I) 
C 

READ(18,190)  PARDI(I)    '-:■■ 
190       FORMAT (F4. 2) 
C 

C     WE  CORRECT  FOR  THE  RELATIONSHIP  BETWEEN  DISTANCE  AND  LENGTH 

C     IN  THE  SEM  PHOTOGRAPH 

C 

PARDI(I)=PARDI(I)/2.0 
END  DO 

CALL  PARSUR ( NP , PARDI , NX, NY , IRN , XTAU , YTAU ) 
END  IF 
C 

C     NOW,  WE  READ  IN  THE  NUMBER  OF  THETA  VALUES  TO  BE  ANALYZED, 

C     AND  THE  NUMBER  OF  SCANNING  LIMIT  VALUES. 

C 

READ (8, 100)  NTHETA.NLIM 
IF  (NTHETA.EQ.O)  GOTO  310 
C 

C     WE  NOW  EVALUATE  THE  PARAMETERS  FOR  EACH  THETA  VALUE,  AND  EACH 
C     SHADING  FRACTION  LIMIT  PARAMETER  VALUE, 
C 

DO  I-1,NTHETA  re 

READ(8,180)  THETA 
WRITE (12, 200) 
200    FORMATC  '.SOC-')) 
WRITE (12, 210)  THETA 

210    FORMATC  ',5X,'THE  INSTRUMENTAL  THETA  VALUE  IS  '.FT.M/) 
C 

C     THETA  IS  READ  IN  HAVING  UNITS  OF  DEGREES,  SO  A  TRANSFORM- 
C     ATION  TO  RADIANS  IS  NECESSARY. 

c  ,-.•.- 

THETA- THETA* PI/ 180.0 
DO  J-1,NLIM  ;  "  -  ■■■  C  ,'-'  s 

READ(8,220)  LIM  -  -   ' 

220       F0RMAT(5X,I4)  „  .  •  .. 

WRITE (12, 230)  ''  /••^■■^  ■ 

230       FORMATC  MOX,36C.')) 

WRITE(12,240)  LIM 
240       FORMATC  M  OX, 'SCAN  LIMIT  IS  ',14, 
&•  BEYOND  EACH  POINT.'/) 

CALL  SHADE  (NX, NY, LIM, XTAU, THETA) 
END  DO 

WRITE( 12,230)  '  '  "  v  ,  ' 

CALL  ACTHET  ( NX ,  NY ,  XTAU ,  THETA )  ' '  '  ~  ' '  •   .  •.■="v-  ■. 
IF  (I.EQ.1  )  THEN  .  '  '  ".  . 

PNORM-0.0 
TNORM=THETA   :    .:}. 
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END  IF 

CALL  AREA  (NX,NY,XTAU,YTAU,SURFAR ) 
CALL  EPOP(NX, NY , XTAU , YTAU , THETA, PNORM, TNORM ) 
WRITE (12, 200) 
END  DO 

WRITE (12, 300)  SURFAR 
300  FORMAT (•  SSX.'THE  APPROXIMATE  SURFACE  AREA  IS  '.FU.S//) 
310  WRITE(12,120) 

IF  (IC0R.EQ.1 )  THEN 

CALL  INOU(NX,NY,LAG) 
END  IF 
STOP 
END 
C 
C 
C 
C 
C 

SUBROUTINE  BROWN  (IX.IY.IRN, N, ISAVE) 
C 

C  THIS  SUBROUTINE  IS  USED  TO  INPUT  DATA  INTO  THE  SURFACE  CONTOUR 

C  MATRIX  SUCH  THAT  ALONG  THE  AXIS  OF  OBSERVATION  EACH  VECTOR 

C  REPRESENTS  A  BROWNIAN  LINE-TO-LINE  FUNCTION.   EACH  VALUE  DIFFERS 

C  RANDOMLY  FROM  THE  PREVIOUS  VALUE  BY  +1  OR  -1 .   THE  RANGE  OF 

C  POTENTIAL  VALUES  IS  GIVEN  BY  [-N,N]. 

C 

COMMON/HARY/H ( 5000 , 2 ) 

ITMP=IRN 

BLOW=0.0 

N1=2*N+1 

DO  1=1, IX 

KR-JINT(RAN(IRN)»FL0ATJ(N1)) 

BL0W=BL0W+(FL0ATJ(KR-N) ) 

DO  J=.1,IY 

H(I,J)=BLOW 

END  DO  ,-M  /^> 

END  DO 

WRITE(12,100)  *  *  J  s  , 

100    FORMATC  SSX.'A  BROWNIAN  LINE-TO-LINE  FUNCTION  IS  THE  DATA 
&  GENERATOR.'/) 

WRITE(12,110)  ITMP       -•■   •-  -- 
110  FORMATC  SSX, 'BROWNIAN  RANDOM  SEED  NUMBER  IS  ',13/) 

WRITE(12,120)  N,N 
120  FORMATC  ',5X,'THE  RANGE  IS  [- '  ,  12,  • ,  ' ,  12,  •  ] '/ ) 
IF  ( ISAVE, EQ.1)  THEN 

OPEN  (UNIT-1  0 , FILE- • BROWN . DAT • , STATUS- ' NEW' ) 
WRITE(10,110)  ITMP 
WRITE(10,120)  N,N 
DO  1=1,  IX  ■;:; 

WRITE(10,130)  I,H(I,1)  '''■' 
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130  FORMATC    M4,5X,F10.5) 

END  DO 
CLOSE    (UNIT=10)  .  ,:'    r-  ■.:  - 

END  IF  -^   -         -S 

RETURN  ■'    -  .     '^v- 

END 
C  -• 

c 
c 

C  '-'"-":  •■  ■ 

SUBROUTINE  SINUS  (IX.IY, ISAVE) 
C 

C        THIS  SUBROUTINE  CALCULATES  A  SURFACE  CONTOUR  WHICH  IS 

C  SINUSOIDAL  IN  BOTH  THE  X-  AND  Y-DIRECTIONS.   THESE  TWO 

C  SINE  FUNCTIONS  ARE  IN  PHASE.  — ' 
C 

COMMON/HARY/H(5000,2)       n  /-  ^  ■'• 


DATA  PI  /3. 1  41592654/  V  {  !';'?*■  jr- 

DO  j-i.iY  ■■  V  V  u '■>:::■■' V 


DO  1-1, IX 

FI=FLOAT(I)/FLOAT(IX)       -   ' 
FJ=FLOAT(J)/FLOAT(IY) 

H(I,J)=SIN(FI»2.*PI)»SIN(FJ*2.*PI) 
END  DO  v: 

END  DO 
WRITE(12,100) 

100  FORMATC  '.SX, 'INPUT  FILE  IS  A  TWO-DIMENSIONAL  SINUSOIDAL 
&  FUNCTION.'/) 
IF  (ISAVE.EQ.1)  THEN 

OPEN  (UNIT=1 4, FILE=' SINUS. DAT', STATUS- 'NEW') 
DO  J-1,IY 
DO  1=1, IX 

WRITE(U,110)  I,J,H(I,J) 
110  FORMATC  '.I^.SX.Iil.SX.FIO.S) 

END  DO 
■  END  DO  ■■■:^'-  '    .' 

CLOSE  (UNIT-14)  ..  V       . - 

END  IF  ■  :  '     ,..;   .  -  .' 

RETURN  "■ 

END 


C 
C 

c 
c 
c 


SUBROUTINE  RUFFLE  (IX, lY, ISAVE) 
C 

C        THIS  SUBROUTINE  CONSTRUCTS  A  SURFACE  WHICH  RESEMBLES  A 
C     RIDGED  POTATO  CHIP.   IT  IS  SINUSOIDAL  IN  ONLY  ONE  DIRECTION. 
C     THE  DIRECTION  OF  OBSERVATION. 
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COMMON/HARY/H ( 5000 , 2 ) 
DATA  PI  /3. 141 592654/ 
C 

C     NWAVE  IS  THE  NUMBER  OF  WAVELENGTHS  OF  THE  FUNCTION  ALONG  THE 
C     VECTOR.   THE  NUMBER  OF  POINTS  IN  EACH  WAVE  IS  THEN  GIVEN  BY 
C     IX/NWAVE. 
C 

NWAVE='50 

DO  1-1,  IX  ,        >?': 

FLI=FLOAT(I)/FLOAT(IX)  "   -■.  .' 

SVAL=SIN(2.0»FL0AT(NWAVE)*PI»FLI)  !    ■ '.  ^   '  :  •. 
DO  J»1,IY 

H(I,J)=SVAL  '■ 

END  DO 
END  DO 
WRITE(12,100) 

100  FORMATC  •,5X, 'INPUT  FILE  IS  A  ONE-DIMENSIONAL  SINUSOIDAL 
&  FUNCTION.'/) 
IF  (ISAVE.EQ.1)  THEN 

OPEN  (UNIT=16,FILE='RUFF.DAT',STATUS='NEW') 

DO  1=1, IX 

WRITE(16,110)  I,H(I,1  )   , 

110       FORMATC  •,I4,5X,F10.5)  ^niU        •>/.-:■ 

END  DO  ^  i<-*v^.i  ;-..;.:  ;:,  , 

CLOSE  (UNIT=16),  ,  *  j .  ^  ...  ^  .  >    .■•>nA  '   '" 

END  IF         '  HIj/'  ^'-  i  WO);^ 

RETURN 

END  .    ^        , 

c       f  -  ^     "•^flVV.Pv^ 

C  I  ,  :  i  '  -''■■       ■     ••  K^V-    .    . 

Q        -  --.  .,  '....V-. -^  .  J  1/  .  :,  \j  v":  .  ..^,, 

c     ■;■-  '''■''  -^  "  ■''. . 

SUBROUTINE  GRATE(IX,IY,ISAVE,R,XT)  "   ■.' 

C 

C        THIS  SUBROUTINE  CONSTRUCTS  A  ONE-DIMENSIONALLY  VARIANT 

C  SURFACE  ANALOGOUS  TO  A  DIFFRACTION  GRATING  OF  BLAZE  ANGLE 

C  R.   THE  VECTOR  CONSISTS  OF  A  NUMBER  PEAKS,  THE  PEAK-TO-PEAK 

C  DISTANCE  BEING  DETERMINED  BY  THE  VALUE  OF  XT. 


C 


COMMON/HARY/H ( 5000 , 2 ) 
DATA  PI   /3. 141592654/ 


C 

C     NPEAK  IS  THE  NUMBER  OF  PEAKS. 
C 

NPEAK-20 

ND-IX/NPEAK 

DIST=XT» FLOAT (ND) 

D1=2.0/TAN(R*PI/180.0) 

D2=DIST-D1 
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c  ■    '■"'  "'  ■    .  '  . 

C     THE  SECONDARY  BLAZE  ANGLE  IS  GIVEN  BY  GAMMA. 
C  -■' 

GAMMA-ATAN(2.0/D2)*180.0/PI 

DO  I=.1,ND  .  ^ 

FLI=DIST*FLOAT(I)/FLOAT(ND) 

IF  (FLI.LE.D1)  THEN  :'" 

DO  J=i1,IY  :  '  ■  :  •  •:  ^  . 

H(I,J)=FLI*TAN(R»PI/180.0)   ■:  '^^.;   .'' 
END  DO 
ELSE 

DO  J-1,IY 

H(I,J)=2.0+(2.0/TAN(R*PI/l80.0)-FLI)» 
4TAN(GAMMA*PI/180.0) 
END  DO 
END  IF 
END  DO 
NM-NPEAK-1 
DO  N-1,NM 
DO  I-I.ND 

DO  j=i,iY         ■■■' •    -  ,  >:^&     "^•:  ■■ 

IPN=I+ND»N  -  '"V  ^:^■  ■^v. 

H(IPN,J)-H(I,J)  .  .  ,  '   ^  ■ 

END  DO 
END  DO  .  ^  ■:''-■' ^'' 

END  DO  >;   ■;"■-  ' 

WRITE (12, 100) 
100  FORMATC  ',5X, 'INPUT  FILE  IS  A  GRATING  FUNCTION.'/) 

WRITE(12,110)  R 
110  FORMATC  ',5X,'THE  BLAZE  ANGLE  IS  '.FT.^/) 

WRITE  (12,120)  GAMMA  ^ 

120  FORMATC  •,5X,'THE  SECONDARY  BLAZE  ANGLE  IS  ',F7.4/) 
IF  (ISAVE.EQ.l)  THEN 

OPEN  (UNIT-1 8, FILE= 'GRATE.DAT' ,STATUS= 'NEW' ) 
WRITE(18,110)  R       -  'T  .—  .w  -— .,    ^ 
WRITE(18,120)  GAMMA  '   ' 

WRITE(18,130)  ND      V  /-:  V  -  -.  ^«'^  /^ 

130    FORMATC  SSX.'THE  VECTOR  REPEATS  AFTER  ',  I4, 'POINTS. '/ ) 
DO  1=1, ND 

WRITE (18, 140)  I,H(I,1) 
140       FORMATC  •,I4,5X,F10.5) 
END  DO 

CLOSE  (UNIT=18) 
END  IF      ■^:  wt'v 

RETURN  .    -    ::v:.  ■-  y      '": 

END 


C 

c 
c 
c 
c 
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C 
C 

c 
c 


SUBROUTINE  PARSUR(N,XI, IX,IY,IRN,DINC ,YT) 

THIS  SUBROUTINE  GENERATES  A  CLOSEST  PACKING  OF  PARTICLES  AS 
A  SURFACE  MODEL,  GIVEN  A  DISTRIBUTION  OF  PARTICLE  SIZES. 

DIMENSION  XI(70),RAD(70),C(70) 

COMMON/HARY/H( 5000 , 2 ) 

ITMP=IRN 

SLEN=0.0  ■'.  ^'C'  -;■ 

DO  I=-1,N  ^  <  '  ' 

RAD(I)=XI(I)»RAN(IRN)/2.0 

C(I)=SQRT(4.0»RAD(I)»(XI(I)-RAD(I))) 

SLEN=SLEN+C(I) 
END   DO  .       .  ,    .,., 

ICUR-1  ,;   '-    i;^■- 

DINC=>SLEN/FLOAT(IX)  '    .  '    •■ 

P0=0.0 
PI =0.0 
P2=0.0 
P3=C(1) 

1=1 

DO  WHILE  (I.LE.IX) 
PO-PO+DINC 

IF  ((PO.GT.SLEN).OR.(ICUR.GT.N))  THEN 
IX=I 
I=IX+1 
GO  TO  10 
END  IF  ■■' 

DO  WHILE  (P0.GT.P3)     ^Vi^i'  ... 

icuR-icuR+1  ■:     '"■     v.^;;- 

pi=P3  '■■-■'  '-'^  ..'  '^":.'"% 

P3-P3+C(ICUR)  -"■■    V- 

END  DO  ■    ,.    '■.;.  ■.:■;   '^    - 

P2=P0-P1  ■  ;  r  !-• 

D-XI  ( ICUR  )/2 . 0-RAD  ( ICUR  )     ' '' ' 

C2=C(ICUR)/2.0 

IF  (P2.GT.C2)  THEN  P2=C (ICUR  )-P2    r 

Q-C2-P2  '::■'■■ 

HE=XI(ICUR)/2.0-Q  .„,  j 

DO  J»1,IY 

H(I,J)-SQRT(4.0*HE»(XI(ICUR)-HE))/2.0-D 
END  DO 

I-I+1 
10  END  DO  '■ 

WRITE (12, 100)  ITMP 
100  FORMAT (•  •,5X, 'RANDOMIZING  SEED  NUMBER  IS  ',13/ ) 

WRITE(12,110)  DINCYT 
110  FORMATC  •,5X,'X  AND  Y  RESOLUTION  ARE  •,F11.8, 
&'  AND  ',F11.8/) 
WRITE (12, 120) 
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120  FORMATC  '.5X,'A  CLOSEST  PACKING  OF  A  SPHERICAL  PARTICLE 
&  DISTRIBUTION  IS  ASSUMED.'/)  ^-AHIICLE 

RETURN 

END        ;-■:         ..^/i-        ';;:• 


c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
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20 
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SUBROUTINE  SHADE  (IX,IY,L,X,T)  ;■  '  v 

.u.T^  SUBROUTINE  ALLOWS  THE  DETERMINATION  OF  THE  FRACTION  OF 
THE  SURFACE  WHICH  IS  SHADED  TO  BE  PERFORMED.  IT  SUMS  THE  AReI 
OF  THE  UNSHADED  PIXELS  AND  DIVIDES  BY  THE  TOTAL  Pm^  Irea 

COMMON/HARY/H(5000,2)/IARY/B(5000.2)  F 
BASE=X»TAN(T) 

IM1=IX-1  ■■:■■.. 

suM=«o.o  ; ;:  ■  i;  ';  : 

DO  J=1,IY  ''-  '^  '.'   -■ 

DO  1  =  1,  IM1  ^::; 

B(I,J)=SQRT(X»X+(H(I+1,J)-H(I,J))»»2) 
SUM=SUM+B(I,J) 

END  DO  .  "  y:. 

END  DO  ,  .       ■'  r  V    > 

F=SUM  ""  ■-,.':-;■■  '':"■■  '" 

DO  J=1,IY  '-:.'^0-;/C, ' 

DO  1=1,  IM1  :  ':  ; 

K=0  -  •-• 

IF  (K  LT.L)  THEN    ^    o/^.. 

IF  ((I+K).GT.IX)  GO  TO  20 

IF  ((H(I,J)+K*BASE).GE.H(I+K,J))  GO  TO  10 
F=F-B(I,J)^-^"*f, 

B(I,J)=0.0 

END  IF  .  'i  .  .  .  *  ". 

END  DO 
END  DO  ■i\}  ,  -■  ..  ti  \'  !. 

FL=F/SUM 
WRITE (12, 100)  FL      „  '  ;N.  ^ 

rSn^'  '•^°^''™^  SHADING  FRACTION  IS  '.FT. 5) 
END  V 


C 
C 

c 

C 

c 


SUBROUTINE  ACTHET  (IX.IY.X.T) 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


V  v 


THIS  SUBROUTINE  CALCULATES  THE  ACTUAL  AVERAGE  ANGLE  OF  PHOTO- 
ELECTRON  ESCAPE  WITH  RESPECT  TO  THE  SURFACE  MEAN.   THIS  VALUE 
DIFFERS  FROM  THE  INSTRUMENTAL  SETTING  DUE  TO  SURFACE  CONTOUR 
EFFECTS.   IN  ADDITION,  THE  AVERAGE  VALUE  OF  THE  SINE  OF  THIS 
ANGLE  IS  CALCULATED,  AS  IS  THE  RATIO  OF  THIS  AVERAGE  SINE  VALUE 
TO  THE  SINE  OF  THE  INSTRUMENTAL  SETTING.   FOR  THE  AVERAGE  SINE 
VALUE,  THE  STANDARD  DEVIATION,  RANGE,  AND  DISTRIBUTION  OVER 
[0,1]  ARE  CALCULATED  AS  WELL. 

^  *  A- 

DIMENSION  DISC  10) 

COMMON/HARY/H ( 5000 , 2 )/IARY/B ( 5000 , 2 ) , F/KARY/SAAR ( 5000 , 2 ) 

DATA  PI  /3. 141 592654/ 

P2=PI/2.0 

ISTEP=0 

A-O.O 

SA-O.O  ,   ■  ■■   '■         -  ,.: 

SASQ=o.o        .  ■  ^?":/:  ; 

SMAX— 2.0  '-■-";;:      :  :. 

SMIN-2.0  ''- 

IM1-IX-1 

DO  1-1,10  '.> 

Dis(i)=o.o       ;:..,•  *"■.■?;- ' 

END  DO  ■'   7 

DO  J-1,IY  '     ''.■', 

DO  1-1, IM1 

IF  (B(I,J).NE.O.O)  THEN   •   v:- 
R=B(I,J)/F 

PHI=ATAN((H(I,J)-H(I+1  ,J))/X) 
AP-T+PHI 

IF  (AP.GT.P2)  AP=PI-AP 
A-A+(AP»B(I,J)) 
SAP=SIN(AP) 

IF  (SAP.GT.SMAX)  SMAX=SAP 
IF  (SAP.LT.SMIN)  SMIN=SAP 
IF  ((O.O.LE.SAP).AND.(SAP.LT.O.I ))  DIS(1) 
IF  ((0.1.LE.SAP).AND.(SAP.LT.0.2))  DIS(2) 
IF  ((0.2.LE.SAP).AND.(SAP.LT.0.3))  DIS(3) 
IF  ((0.3.LE.SAP).AND.(SAP.LT.0.4))  DIS(4) 
IF  ((0.4.LE.SAP).AND.(SAP.LT.0.5))  DIS(5) 
IF  ((O.5.LE.SAP).AND.(SAP.LT.0.6))  DIS(6) 
IF  ((0.6.LE.SAP).AND.(SAP.LT.0.7))  DIS(7) 
IF  ((0.7.LE.SAP).AND.(SAP.LT.0.8))  DIS(8) 
IF  ((0.8.LE.SAP).AND.(SAP.LT.0.9))  DIS(9) 
IF  ((0.9.LE.SAP).AND.(SAP.LE.1.0))  DIS(10 
SA=SA+(SAP«B(I,J))       :  • 
SAAR(I,J)=SAP 

ISTEP-1  :'  i.  :^ 

ELSE  -^  .-'v^\  ^ 

SAAR(I,J)=0.0        "-::'--,    : 
END  IF  .»>";:'-,,    ■  ■' '• 


=DIS(1 )+R 
=DIS(2)+R 
=DIS(3)+R 
=DIS(4)+R 
=DIS(5)+R 
=DIS(6)+R 
-DIS(7)+R 
'=DIS(8)+R 
=DIS(9)+R 
)-DIS(10)+R 
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C 
C 
C 


C 
C 

c 
c 
c 

c 

c 
c 
c 
c 


END  DO 

END  DO      -  •   ■    :^  :         _ 

IF  (ISTEP.NE.O)  THEN  fe  •  ■        •■--;' 
A=A/F 
SA=SA/F 
RS=SA/SIN(T)  .     .  ^: 

WE  CONVERT  ATHETA  BACK  TO  DEGREES.      ,  ■ 

A=A*180.0/PI 
WRITE(12,100)  A 

^^\   «.in°c!?^^^'  ''5X''™E  AVERAGE  ACTUAL  THETA^  VALUE  IS  • 

WRITE(12,110)  SA 

^^^i.Fsfs^^^^'  '•^^•'™^  ^^^'^AGE  VALUE  OF  SIN  (ATHETA)  IS  • 

WRITE (12, 130)  SMIN.SMAX  "'■        ' 

130    FORMATC  '.SX.'THE  RANGE  IS  '.FS.S.'  •  F8  5/) 
WRITE  (12,140)  .POO/; 

''*°    WRmnP  Wn7nrln?'^''V?   '""^^^'^  °^  ''^   (ATHETA)') 
WRITE(12,150)  DIS(1),DIS  2  ,DIS(3).DIS(i1)  DlSfS^ 

&DIS(6).DIS(7).DIS(8).DIS(9).DIS(10         ^  ^^' 
150    FORMATC  M0(2X,F5.3)/) 

WRITE(12,160)  RS  .  ^  , 

'^\  is'?'^y^;/'^^''™'  ^^^"^  0^  sinYatheta)/sin  (THETA) 

ELSE 

WRITE(12,170) 
^^°  END^?F^^^'  ''5^' 'ACTUAL  THETA  VALUE  CANNOT  BE  COMPUTED.'/) 

RETURN 

END         •  .v  -  '  ■     :■■■  '. 


SUBROUTINE  AREA  (IX, IY,X,Y,S)         *     '  ' 

SUnlfrt   ^R^f  """Jf  n^ic°''^  ^"^   APPROXIMATE  CALCULATION  OF  THE 
SmES  OR  ?R;.MrJ.c°^?,™^^  ^^  DIVIDING  THE  SURFACE  INTO  A 
SERIES  OF  TRIANGLES,  CALCULATING  THEIR  AREAS,  AND  SUMMING. 

COMMON/HARY/H(5000,2)/IARY/B(5000.2).F/JARY/HAAR(5000,1) 
HARE=0.0 

jMi=iY-i  *:; 

iMi=ix-i     ^r^:-^''^'         :  .•  •  ;■■;  ::■■ 

RMS-SQRT  (X»X+Y»Y)     'vvV^  ; , 
DO  J=1,JM1 
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DO  1=1, IM1 

AY=ATAN(ABS(H(I,J+1  )-H(I,J))/Y) 
AX=ATAN(ABS(H(I+1,J+l)-H(I,j+i))/x) 
BY-ATAN(ABS(H(I+1,J+1)-H(I+1,j))/y) 
^,.  BX=ATAN(ABS(H(I+1,J)-H(I,J))/X) 

:::;    c=atan(abs(h(i+i  ,j+i  )-h(i,j))/rms)       .-.^ 

SA=0.5»(X/C0S(AX)+Y/C0S(AY)+RMS/C0S(C))    ■-    > 

SB=0.5»(X/COS(BX)+Y/C0S(BY)+RMS/CGS(C))    ^ 

TARE=TARE+0.5*X»Y/(C0S(AX)»C0S(AY))+0  5»X»Y/ 
&(COS(BX)»COS(BY))  -^VHx;;  uo  x  1/ 

HI=SQRT(SA»(SA-X/COS(AX))*(SA-Y/COS(AY)) 

^(sti;i/co1(cJn*''"'^''*^''"'''°'^''^^*^ 

HARE=HARE+HI 

IF  (B(I,J).NE.0.0)  THEN 

HAAR(I,J)=HI 
ELSE 

HAAR(I,J)=0.0 

'■""    END  IF  -;-  .: 

END  DO  :'-■'- 
END  DO  ;  . 

s=HARE  -.'-:';: 

RETURN 
END 


C 
C 

c 

c 

c 
c 
c 
c 
c 
c 
c 
c 


if)  I  ^ 


SUBROUTINE  IN0U(IX,IY,K) 


.rr^rZ^l^  SUBROUTINE   PERFORMS   AN  AUTOCORRELATION  ON  THE   INPUT   DATi 
eLh  VECTOR    Tr.l\''  '^^'ORMim  ONE-DIMENSIONAL  OPERA^IoL  oT 
!  pin        ^   ^^  ™^   X-DIRECTION,    AND  SUBSEQUENTLY  AVERAGING  THEM 
A  GOOD  MAXIMUM  LAG  VALUE    IS  TYPICALLY    1 0?   OF  THE  TOTAL   LENgS  OP 
oS?pu?Tt;  ,™^/"^°^°««^^ATI0N  VECTOR   Is  ALsTwRm\^'oNTTH°E 

DIMENSION  ACV(500).AC(500),PACV(1),WKAREA(1)  '  ' 

DIMENSION  W ( 5000) , CORSUM ( 500) , COVSUmS) 
COMMON/HARY/H(5000,2)  >: 

WRITE(12,100) 

100  FORMAT(.  •///.5X.. AUTOCORRELATION  VECTOR  DATA  IS  SHOWN  BELOW.-//) 

CORSUM(I)-0.0  ' 

COVSUM(I).0.0 
END  DO 

SUM=0.0  ' 
DO  J=1,IY 

AMEAN=0 . 0 

DO  I»1,IX 
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C 

c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 


W(I)=H(I,J)  ^-  -.'v-  '  ; 

AMEAN=AMEAN+H(I,  J)  '  '^  -  'vT' 

SUM=SUM+H(I,J)  -:'^:  -^t  . 

END  DO            ';"•";  ■    ,:''  r  ';  ; 

AMEAN=AMEAN/IX 
D0^M,™^°^"'^^'^'^'5'AMEAN,VAR,ACV,AC.PACV.WKAREA) 

COVSUM (M ) =C0 VSUM (M ) +ACV (M )/FLOAT (K ) 

C0RSUM(M)=C0RSUM(M)+AC(M)/FL0AT(K) 
END  DO 

END  DO 

SUM=SUM/FLOAT(IX»IY)        .,.,.. 
WRITE(12,110)  SUM 

110  FORMATC  '.SX.'THE  SURFACE  MEAN  HEIGHT  IS  '  El  2  5//) 
WRITE(12,120)  .ni^O/z; 

' '°  wTmlJs.  i'bO  ) ';:;;;?^.7"'°'°^^«^^^^^  • '  «^'  -AUTOCORRELATION.  / ) 

130  FORMATC  ',5X.I3,2(10X,E12.5))  .    .  ^   .. 
DO  1-1 ,K 

WRITE(12,130)  I.ACV(I).AC(I) 
END  DO 
RETURN  /■ 

END  v..::-:. 


'  <is^> 


SUBROUTINE  EPOP(IX,IY,XT, YT.T, PN. TN) 


ELECTRONf  WHTPH  ?I«'or^   ^°  CALCULATE   THE   OVERALL  POPULATION  OF 

IS  REPRELSpn   .f  ?  lltr""  ^"^  ^^^^'"''^  2^^^S.      THIS  POPULATION 
IS  REPRESENTED  AS  A  PERCENTAGE   OF  THE   POPULATION  FOR   THE   FTR^T 

DEGREES?.™'  ''''  ''  ''''''  ''  ''  ''''''''    (TmclLLY     89 

COMMON/JARY/PIXAR(5000,1)/KARY/DEPAL(5000,2) 

PSUM=0.0  .;■;  ,.v,^ 

DSUM=0.0  V 

FLAR=0.0  !- 

IM1=IX-1 

JM1=IY-1 

DO  J=1,JM1 

DO  I=.1,IM1  -5  ^ 

POP=POP+PIXAR(I,J)*DEPAL(I,J)  /i^,-; 

PSUM=PSUM+PIXAR(I,J)  .  .  '. 

DSUM=DSUM+DEPAL(I,J)     ' 

FLAR=FLAR+XT»YT 
END  DO 

END  DO  ) 'i   I  Tr^   "yirt'^^t- 

-  ■-f_  iWjIuf  y^wA  V-w^"  '"■%■' 


'-^  f* 


1*!  /  •  ->»  .♦»;■  >.  !«■*■ 
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IF   (PN.NE.0.0)   THEN 

PERFL«100.0*P0P/FLAR 

POP=POP/PN  -,v     ,'  V.      :• 

ELSE  ,     ..„  V,'         ..     ^^:r" 

PN=POP  ' 

PERFL=100.0*PN/FLAR 

P0P=1.0  ■       :,v.M^    -    --:■  '     -  ^  -J-   •. 

END   IF 
WRITE (12,  98)    PSUM 

98  FORMAT (•    '.SX.'THE  SUM  OF  UNSHADED  PIXEL  AREAS   IS    ',F12.5) 
WRITE (12,99)    FLAR 

99  FORMAT (•    ',5X, 'COMPARABLE   FLAT  SURFACE  AREA   IS    •,F12.5/) 
WRITE( 12,97)  DSUM 

97      FORMATC    ',5X,'THE  SUM  OF  NORMALIZED  PENETRATION  DEPTH 
&IS    •,F12.5/) 
WRITE ( 1 2, 1 00 )   P0P*1 00 . 0 , TN»57. 29578 

100  FORMATC    ',5X,'F0R  THIS  ANGLE,    THE  POPULATION   IS    ',F7.2, 
i'it   OF  THE  VALUE   AT    ',F5.2,'    DEGREES.'/) 

WRITE(12,120)    PERFL 
120   FORMATC    •,5X,'THIS  VALUE    IS    •,F7.2,'?   OF  THAT  FOR 
&  A   FLAT  SURFACE.'/) 
RETURN  ■   ■  :  ; 

END 


'1 


.*»/ 


i**  ;"%  «: 


■^ 


-Hr    *f      -'<-     -"S 


\.i    ^,-     ^       i-f^- 


A  > . 


"*  g-'-^'0'%ji 


APPENDIX  B 
NILT  PROGRAM  SOURCE  CODE 


vi  c^i  1/.M  U  ;J  <r>v 


c 
c 

C  .,     :,  ■  '  ^      '  :  -  :  ;, 

c  "  .-   ■  ■■  --^  ■■■■ 

c 

C  DEPTH  PROFILE  DETERMINATION  PROGRAM 

C 

C  PAUL  C.  MCCASLIN 

C 

C        THIS  PROGRAM  PROVIDES  FOR  INPUT  OF  A  DATA  VECTOR  OF  XPS  PEAK 

C  AREA  VALUES,  OR  RATIOS  OF  SUCH  VALUES,  AS  A  FUNCTION  OF  THE 

C  INSTRUMENTAL  SETTING  FOR  THE  PHOTOELECTRON  ESCAPE  ANGLE.   GIVEN 

C  THIS  DATA,  THE  PROGRAM  ALLOWS  FOR  INPUT  OF  AN  INITIAL  GUESS  TO 

C  THE  CONCENTRATION,  OR  CONC.  RATIO,  DEPTH  PROFILE  OF  THE  SPECIES 

C  UNDER  INVESTIGATION.   THIS  INITIAL  GUESS  MAY  BE  GENERATED  BY  THE 

C  PROGRAM  IF  ONE  OF  THE  ANALYTIC  FUNCTIONAL  FORMS  IS  CHOSEN.   FROM 

C  THIS  POINT,  THE  PROGRAM  ITERATES  TOWARD  A  MORE  EXACT  SOLUTION, 

C  ITS  APPROACH  BEING  BASED  ON  SEVERAL  USER-VARIABLE  PARAMETERS. 
C 

DIMENSION  XI(12),YI(12),C0(12,12) 
DIMENSION  X0UT(12),YIN(12),AXB(12) 
DIMENSION  SSC(12)  ■  ■ • 

REAL  K 

COMMON/ONE/ROOTLEG (12,2) 

C0MM0N/TW0/NP,S(12),YNT(12),S0LIN(12),Z(12) 
OPEN  (UNIT=12,FILE«'DYN.DAT',STATUS='0LD' ) 
OPEN  (UNIT=10,FILE='D0UT. DAT '.STATUS- 'NEW') 
DATA  PI  /3. 141 592654/ 
READ(12,90)  NFL,N0UT 
90   F0RMAT(2(5X,I1  )) 

C  IF  NFL  =  0,  THEN  THE  INITIAL  GUESS  TO  THE  DEPTH  PROFILE  IS 

C  INPUT  VIA  THE  MAIN  DATA  FILE. 

C        NFL  -  1 ,  THEN  THE  INPUT  DATA  FILE  IS  GENERATED  BY  ASSUMING 

C  AN  ANALYTIC  FORM  FOR  THE  LAPLACE  TRANSFORM  INTEGRAL. 

C 

C  IF  NOUT  -  0,  THEN  THE  ITERATION  RESULTS  WHICH  ARE  OUTPUT  ARE 

C  EVENLY  SPACED. 

C        NOUT  =  1 ,  THEN  THE  ITERATIONS  WHICH  POSSESS  DESIRED 

C  CUTOFF  VALUES  FOR  THE  SUM  OF  SQUARES  OF  ERRORS  ARE 

C  OUTPUT. 

C 
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READ(12,100)  NP      "'  ^  ■   -    •       ,'-^':\ 
100  F0RMAT(5X,I2) 
C  ■         '      • 

C     WE  READ  IN  THE  MAXIMUM  NUMBER  OF  ITERATIONS,  AND  THE  VALUE  FOR 
C     THE  COURANT  PARAMETER.   THE  GREATER  THE  COURANT  PARAMETER,  THE 
C     FASTER  THE  APPROACH  TO  THE  SOLUTION,  AND  THE  GREATER  CHANCE  OF 
C     NOISE  BEING  INTRODUCED. 
C 

READ(12,110)  NITER, COUR  \:       .      ■  -•.'•■ 

110  F0RMAT(5X,I4,5X,F7.5)  ■   •' 
IF  (NOUT.EQ.O)  THEN  "' 

C 

C     NINTER  IS  THE  INTERVAL  BETWEEN  PRINTOUTS  OF  ITERATION  RESULTS. 
C  i-i.     ./"'.»"  '•*'  *---^  i 

READ (12,1 11 )  NINTER 

111  FORMAT  ( 5X ,  14 )        '  i  i^  f'  '  r  f   f^  \ " 

SSCUT-O.O      '  v«^  -i  '^/'  i  '■-  'i  i  ^.'   i ' 
NSS-1 
ELSE      --   •   -  -  -^ '^f  f.  :'■   ■••>  "^    "  ".  ■" 

c  :  j  ^.:'  ■  - 

C     NSS  IS  THE  NUMBER  OF  SUM  OF  SQUARES  VALUES  FOR  WHICH  RESULTS  ARE 

C     EXPECTED.   THE  ACTUAL  SS  VALUES  ARE  STORED  IN  SSC(I). 

C 

READ (12, 100)  NSS 

DO  1=1,  NSS  ■  ■  '''^^>-[\ 

READ(12,112)  SSC(I)  ,  -   \  '" 

112  FORMAT(6X,E9.4) 

END  DO  ..■••■  - 

NiNTER=i  -'.:  '''^-   :: 

END  IF  -,. 

c  -■■  •  -■'  ■"'■"'  •  "  ■ 

WRITE(10,122)  ■ 

122  FORMAT  (//'  MOX, 'S ',1  3X, 'CORR.  INT.V) 

DO  1-1, NP 

READ(12,120)  XKD.YKI) 

120    F0RMAT(5X,F5.2,5X,E14.8) 

S(I)-1.0/SIN(XI(I)»PI/180.0)       <^-.' 
C 

C     ANY  FUNCTIONAL  CORRECTION  OF  THE  INTENSITY  DATA  SHOULD  BE 

C     MADE  ON  THE  FOLLOWING  PROGRAM  LINE. 

C 

YNT(I)-YI(I) 

WRITE(10,125)  S(I),YNT(I) 
125    FORMATC  ' ,  5X,F1  2.  9,  5X,E1  4.  8)     ■■';:;•:''     ,  - 
END  DO  -   -.;,  -  .. 

C 

C     FIRST,  WE  LOAD  IN  THE  PROPER  SET  OF  ROOTS  AND  WEIGHTS  FOR 
C     SETTING  UP  THE  GAUSSIAN  QUADRATURE  COEFFICENT  MATRIX. 
C  ■,  .■  .--•..  -  -,      -'„■-■  ..  ,  ^ 
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:•->*'■ 


193 
194 


C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


CALL  ROOTEN    -  '  '  '^ ''^^' .;     A-  ;  ;  ■■ 

DO  i=i,NP    ■'■"'._  /!----■;  .  ,  ,  ■  ^-'^'  .      "■/--■ 

DO  J=1,NP 

C0(I,J)=R00TLEG(J,2)*R00TLEG(J,1)»*(S(I)-1  0) 
END  DO  .  ,       '^ 

Z(I)=-L0G(R00TLEG(I,1 ))   ''  - 
IF  (NFL.EQ.O)  THEN       '-''^'.   '.. 

READ(12,193)  SOLIN(I)        "V 

F0RMAT(10X,E14.8) 

IF  (I.EQ.1)  WRITE(10,194) 

A  DA?rVECT0R:^) '''"''"''"  ''"'^^^  ^°  '"^^''^'^  "^^  "^^'^  ^^^^^0 

END  IF        '     '  "' 


&  AS 


IF 


THE  VAlJe'of'nSEl!'"'"  "'"''''''  TEST  ' FUNCTIONS  IS  GOVERNED  BY 
IT  IS  A  STEP  FUNCTION. 

WE  ASSUME  A  SLOPING  INCREASE  WITH  A  CUTOFF 
IT  IS  A  POLYNOMIAL  WITH  CUTOFF. 
AN  OVERLAYER  DISTRIBUTION  IS  ASSUMED. 
WE  USE  A  POWER  LAW  FIT. 
A  POWER  LAW  FIT  WITH  CUTOFF  IS  ASSUMED 
IT  IS  AN  EXPONENTIAL  FIT. 
WE  USE  A  LINEAR  FUNCTION  WITH  A  NON-NEGATIVE 


NSEL  = 
NSEL  = 
NSEL  = 
NSEL  - 
NSEL  = 
NSEL  = 
NSEL  - 
NSEL  = 
INTERCEPT 


1, 
2. 

3, 

^i, 
5, 
6, 
7, 
8, 


END  DO 

IF  (NFL.EQ.1)  THEN 
READ(12,100)  NSEL 
CALL  ANALYT(NSEL) 

END  IF 


•'.-^ ' 


'■o^'   •■  ..,«.S! 


DO  I  =  1,NP  ::;  •'        -'  ■'}  ':V 

YIN(I)-YNT(I)       ^  ■'■%".■■  -     •'■  i 

XOUT(I)=SOLIN(I)    -v  r  :.  v 

END  DO  -  -^-'y  -■' 

CALL  AXBSQ  (CO ,  XOUT ,  AXB ,  SS )  ' '-  '- "  ^  ^ 

CALL  WRITOUT(XOUT,AXB,SS,0)  ./^ r  ^  ■ 

IF  (N0UT.EQ.1)  SSCUT-SSCd)  "^^^v  : 

M»i  ,  ■  ■.>;"': 

DO  1  =  1,  NITER  "     '■^- 

CALL  SMOOTH (CO, YIN, XOUT, COUR) 

CALL  AXBSQ(CO,XOUT,AXB, SS)  ;   <:  -^ 

ITEST-MOD(I,NINTER)  "  ^v/  ' 

W1  =0  .         :  v.  '\ ,.       ,,  ' 

IF  ((ITEST.EQ.O.AND.NOUT.EQ.O).OR 
4(SS.LE.SSCUT.AND.NOUT.EQ.1.AND.W1.EQ.'o))  THEN 
CALL  WRITOUT(XOUT,AXB,SS,I) 
M=M+1 

IF  (N0UT.EQ.1)  THEN 

WRITE(1G,166)  SSCUT    \'-^  ^r?^^-,  . 
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166  FORMATC    '.SX, 'CUTOFF  ERROR  VALUE   IS    ',£9.^1) 

IF    (M.LE.NSS)    THEN 
SSCUT=SSC(M) 

wi=i       -    ■  :     - 

END   IF 
END   IF  ■     i 

END   IF  .'..:. 

DO  j-i,NP  ; 

soLiN(j)=xouT(j)  :;^>*^ 

YIN(J)=YNT(J) 
END  DO 

END  DO 

WRITE(10,160) 
160   FORMAT(///'    ',20X, 'FINAL  RESULTS'/) 

WRITE(10,165)  GOUR 
165   FORMATC    ',5X,'THE  COURANT   PARAMETER  IS    '.FT.S) 

CALL  WRITOUT(XOUT,AXB,SS, NITER) 

STOP 

END  :'   . 

c 
c 
c 
c 

SUBROUTINE  AXBSQ(C ,X0, AXB.SS)   '  " 
C 

C        THIS  ROUTINE  CALCULATES  THE  SUM  OF  SQUARES  OF  THE 

C  MATRIX  MULTIPLIED  BY  THE  ANSWER  VECTOR  MINUS  THE  ACTUAL 

C  SOLUTION  VECTOR,  AND  WRITES  IT  OT  THE  OUTPUT  DATA  FILE. 
C 

DIMENSION  C(12.12),X0(12).AXB(12) 

C0MM0N/TW0/N,S(12),Y(12),X(12),Z(12) 
SS-0.0 

DO  1-1, N 

suM-0.0         '  .  ■■; ■^;^;; 

DO  J-1,N 

SUM=SUM+C(I,J)*XO(J) 
END  DO 

AXB(I)-(SUM-Y(I))»*2 
SS-SS+AXB(I)  :..     . 

END  DO 

RETURN  >■  , 

END 
C  ■-■- •■■;^-.,  ■•■ 

C 

c  .  ''TJ  .  .. ^  .  .'  /  _  .-'-■■    ■ 


SUBROUTINE  WRITOUT(X,AXB,SS, ITER ) 
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C        THIS  ROUTINE  WRITES  THE  DATA  FOR  A  PARTICULAR  ITERATION 

C     ONTO  THE  OUTPUT  DATA  FILE. 

C 

DIMENSION  X(12),AXB(12) 
COMMON/TWO/  N,S(12),Y(12),S0L(12),Z(12) 
WRITE(10,140)  ITER 
140  FORMAT(//'  ', 5X, ' ITERATION  ',14/ ) 

WRITE(10,167) 
167  FORMATC  •,9X,'Z/LAMBDAM2X.'CONC.  RATIOM  IX,  •  (Ax  -  b)"2'/) 
DO  1-1, N 

WRITE(10,170)  Z(I),X(I),AXB(I) 

170   FORMATC  '.scsx.Eie.g)) 

END  DO  ."f 

WRITE(10,222)  SS  ;  '^ 

222  FORMATC  •/,4lX,'SUM  -  ',E16.9) 
RETURN 
END 
C  .   ^  . 

c  ;.'"  ■  ■•;   ■■ "-  ■■-■ 

C  \.-: 

C  J.    i 

SUBROUTINE  SMOOTH(A,B,X,XLAM)   '^•"^'■"'   .   , 
C 

C        THIS  SUBROUTINE  PERFORMS  THE  ACTUAL  DYNAMIC  PROGRAMMING. 
C 

DIMENSION  A(12,12),B(12),X(12),XIDENT(12,12),XK(12), 
&AKRN(12,12),Q(12,12),P(12),ALPHA(12,12),RH0(12) 
C0MM0N/TW0/N,S(12),Y(12),C(12),Z(12) 
C  ..  .,  '  _ 

DO  1-1, N  .  .:-':,: 

DO  J-1,N  '  ■ -vC  ^  ■•■■' 

XIDENT(I,J)=0.0 
END  DO 

XIDENT(I,I)-1.0  "<:.. 

END  DO 

c  .      ■  '       •  >.^:?':'::  ■/ 

XK(1)-0.0  •         v;        '  '^ 

DO  1-1, N 

XK(1)-XK(1)+A(I,1)»»2    •  '■   ,  • 

END  DO  '  ^.ir      V-   'v'" 

DO  1-1, N  -  •.        -' 

DO  J-1,N 

AKRN(I,J)-(A(I,1)»A(J,1))/(XLAM+XK(1)) 
END  DO 

END  DO 

c  '■>_■'"-  '^  ,' . 

DO  1-1, N  -  "  ■ 

DO  J=1,N 

Q(I,J)=XIDENT(I,J)-AKRN(I,J) 
END  DO 
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C 

c 


END  DO 

DO  1=1 ,N 

P(I)=(-XLAM*C(1 )»A(I,1))/(XLAM+XK(1)) 
END  DO 

R=XLAM«(C(1 )*»2)-((XLAM»C(1  ) )»*2  )/ (XLAM+XK(1 )) 

DO  1-1, N 

ALPHAdJ  )=A(IJ  )         ..."  .  ;;v 
END  DO  ■   ■  ^ 

RHOd  )-0.0 

DO  IXX1-2,N  /    :     1  .     '• 

DO  1=1  ,N  '<:^fr:\"  -■ 

ALPHAd.IXXI  )=0.0 

DO  J»1,N  ■  •  ' 

ALPHAd.IXXI  )=ALPHA(I,IXX1  )+Q(I, J)»A(J,IXX1  ) 

END  DO 
END  DO 

RHOdXXI  )=0.0  I    .     '  . 

DO  1=1  ,N 

RHOCIXXI )=RH0(IXX1 )+P(I)»A(I, IXX1 ) 
END  DO 

XKCIXXI )=0.0 
DO  1=1  ,N 

XKCIXXI  )=XK(IXX1  )  + (ALPHAd.IXXI  )»Ad,IXX1  )) 
END  DO 
DO  1=1, N  '  .    •     . 

DO  J»1,N  '■''■  '''■ 

AKRN  (I,  J  )=  (ALPHAd.IXXI  )*ALPHA(J,IXX1  )  )/(XLAM+XK(IXX1  )) 
END  DO 
END  DO 
DO  1=1  ,N  ■••  ^' 

DO  J=1,N  ^.  .;  -■■• 

Q(I,J)=Q(I,J)-AKRN(I,J)         ^  ■^•'■:. 
END  DO 
END  DO 
DO  1=1 ,N 

P(I)=P(I)-((XLAM»C(IXX1 )+RH0(IXX1 ))»ALPHA(I,IXX1  ))/(XLAM+ 
&XK(IXX1  ))  ^  .;' 

END  DO 

R»R+XLAM»C(IXX1  )«»2-((XLAM*CdXX1  )+RH0(IXX1  ))»»2)/(XLAM+ 
&XK(IXX1 )) 
END  DO 

DO  IXX2=1,N  .'■      '■■  '<     •/';•. 

IXX3=N+1-IXX2 

X (IXX3 )= (XLAM»C ( IXX3 )+RHO ( IXX3 ) ) / (XLAM+XK (IXX3 ) ) 
DO  1=1, N 

X(IXX3)=X(IXX3)+(ALPHA(I,IXX3)*B(I))/(XLAM+XK(IXX3)) 
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END  DO      ■.;-  ■■  " 

c 

DO  1=1 ,N 

B(I)=B(I)-X(IXX3)*A(I,IXX3) 
END  DO  ,  ,.  , 

END  DO  '  '      .  -•  ; 

RETURN       ■  -  ■■  ■ 

END  ^  ^.  -*   .  ^ 


^  r  .,',. ..«  (^'  '  \  \I  » ,.4^i   . 

SUBROUTINE  ROOTEN   ,.;  :.  -  ,.  j-.  ,^  ^,  -^  -^  ., 


C 


c 

C        THIS  ROUTINE  SERVES  THE  FUNCTION  OF  LOADING  IN  THE  PROPER 

C     SET  OF  SHIFTED  LEGENDRE  ROOTS  AND  WEIGHTS. 

C 

COMMON/ONE/ROOTLEG ( 1 2 , 2 ) 

C0MM0N/TW0/N,S(12),Y(12),S0L(12),Z(12) 

CHARACTER»9  DFILE 

DFILE=' LEGXX.DAT'  ^ 

NX-N-2  .y 

GO  T0(30, 40, 50, 60, 70, 80, 90, 100, 11 0,1 20). NX 
30   DFILE(4:5)='03' 

GO  TO  160  ' 

40  DFILE(4:5)=''04'  .   " 

GO  TO  160 
50   DFILE(4:5)-'05'  '■•..       .*:.,.;,.: 

GO  TO  160  '  ■  ■:.    :       .'■'     ■         '.: 

60  DFILE(4:5)='06'  '  ;-     ;      ^%<-       ' 

GO  TO  160  '     ,-      '    ■  '^^  -  "- 

70   DFILE(4:5)='07'  '      '  •  :       '. 

GO  TO  160  :' -'  ^'•■"     \. 

80  DFILE(4:5)»'08'  ;      -     :' 

GO  TO  160  ...     •    .     .     '•■■-.  J- •^.     !■■{■: 

90   DFILE(4:5)-'09'  -"' /i   . 

GO  TO  160  ';     -  .-  :v,:  •^  :  ■,^'      '  ,  .'^  •■■-■■  :_ 

100  DFILE(4:5)=-'10'   '  "'. 

GO   TO    160  ':yr'--  '-^    ' 

110   DFILE(4:5)-'11  •  • -^     ' 

GO  TO  160  .  ;•:;•  V        ^  :  ■''         ■- 

120   DFILE(4:5)-'12' 

160   OPEN    (UNIT-14,FILE-DFILE,STATUS='0LD') 
DO   10   1-1, N 

READ(14,200)   ROOTLEG(I, 1 ) ,R00TLEG(I,2) 
C 

C  R00TLEG(I,1)   COMPRISES   THE  ROOTS  OF  THE  SHIFTED  LEGENDRE  POLYS. 

C  R00TLEG(I,2)   COMPRISES  THE   CORRESPONDING  WEIGHTS. 

C 

200  F0RMAT(2(5X,F15.14))  1 
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10   CONTINUE 
RETURN 

END        '    '  ■''  :,^'-':^-' ■■ 

c  ,  .  '...;:-;'': 

c     "■-..■ 
c 

c     ,  •.■;.,■,, :  -..  .  , 

c 

SUBROUTINE  ANALYT(ISEL) 
C 

C        THIS  ROUTINE  ALLOWS  THE  CHOICE  OF  ONE  OF  SEVERAL  ANALYTIC 

C  FORMS  FOR  THE  LAPLACE  TRANSFORM  EQUATION.   GIVEN  THIS  CHOICE, 

C  THE  BEST  FIT  ASSUMING  THIS  FORM  IS  DETERMINED  USING  A  LEAST 

C  SQUARES  ANALYSIS. 

DIMENSION  SS(12)  j^ 

C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12),Z(12)      ^-      ,      '' 

C0MM0N/THREE/YT(12)  '  ;;  "  , 

IF    (ISEL.LT.4)   READ(12,100)    EXK        \  C'  ;■■■<■> 

100  F0RMAT(5X,E14.8)  "' "-  .:-": 

C 

C     EXK  IS  THE  EXPERIMENTAL  VALUE  OF  K.   ' -" 
C 

IF  (ISEL.EQ.6)  READ (12, 100)  EXZ    '  ---  ^■-'•»*-  - 
C 

C     EXZ  IS  THE  EXPERIMENTAL  CUTOFF  VALUE  FOR  Z/LAMBDA. 
C 

GOTO(10,20,30,40,50,60.70,80),ISEL 
10  CALL  STEFUN(EXK) 

GOTO  90  ;•  :'     ':?■  • 

20  CALL  SLOINC(EXK) 

GOTO  90 
30  CALL  POLCUT(EXK)      ....  •  ,      "- 

GOTO  90  .  J,   ;•'    ■'■-.;; 

40  CALL  OVERLAY  .  '  ' 

GOTO  90  ■'  '' 

50  CALL  POWLAW  .  ''  ' 

GOTO  90     .  ':'''■  '■  ;• 

60  CALL  POWCUT    "  .  -  '  .  ;> 

GOTO  90       .  ,^  .;  .    "  '  ;  ■ 

70  CALL  EXPONE       ''^    ••.:  .  '  "  ' 

GOTO  90  '  '  ■  .  "  V  •■: 

80  CALL  LININT  .;  '  '  ' 

90  SSQ-O.O  '  ;  "• 

DO  I»1,N 

SS(I)-(YNP(I)-YT(I))»»2 

SSQ-SSQ+SS(I)  •.   ,  •  ^ 

END  DO  ■...--■ 

RETURN 

END  '  . 
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C 
C 

c 
c 
c 

c 
c 
c 


SUBROUTINE  STEFUN(EXK) 

THIS  IS  THE  ROUTINE  ACCESSED  FOR  A  STEP  FUNCTION  DISTRIBUTION. 


C 

c 
c 
c 
c 


C0MM0N/TW0/N,S(12),]fNP(12).S0LIN(12),Z(12) 
C0MM0N/THREE/YT(12) 

DIMENSION  XL(12),YL(12) 

DO   1-1, N  .  ^  :- 

yL(I)=LOG(S(I)»YNP(I))        •    --  '^ 
XL(I)=S(I) 
END  DO 

CALL  LS ANAL (N, XL, YL.RINCEP, SLOPE, CO) 

C=EXP(RINCEP) 

THEOK=-SLOPE 

WRITEdClOO)    EXK  '  "^'^^       ' 

'°°  Soo'l'li)^'"''"''''"''''^''^*''   ^  ^^"^""^   =    ''^^^-^^ 

110  FORMATC/'  ',5X,'THE  EQUATION  IS  OF  THE  FORM') 

WRITE(10,120) 
120  FORMATC  M5X.'N(Z)  =  0  IF  0<Z<K' )      -■  • . 

WRITE(10,130)  -;,, 

130  FORMATC  ',22X,'C  IF  Z>K' )       ':  '      ''    :i:;%' 

WRITE(10,140)  C 
140  FORMATC/'  ',5X,'C  -  ',E16.9)        V  • 

WRITE(10,150)  THEOK 
150  FORMATC  •.5X,'K  =  ',E16.9) 

PERCENT=100.0»ABS(THEOK-EXK)/EXK 

WRITE (10, 160)  PERCENT 

160  FORMAT(/.  .,5X,'J  ERROR  BETWEEN  EXPTL.  AND  THEOR.  K 
atfro.2,'%'; 

DO  1=1, N 

YT(I)=C»EXP(-THEOK»S(I))/S(I) 

IF  (Z(I).LE. THEOK)  THEN      .  .  ,•  :- 

SOLIN(I)=.0.0  •:"  ' 

ELSE  ^7J   ■-'-■'^^ 

soLiN  ( I )  -c.^.«*. .,*,w  «^:  •  TS  0  0 1" '  '■'   .  ■  - . 

END  IF  ■■  .*^  ■ 

END  DO  .  ;,  .^/.  ^   , -^  ^  . 

RETURN       .i^OUmhini. 

END  :• 


«iii 


V^.fc.-  ,^i    }.^     ,^ 
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SUBROUTINE  SLOINC (EXK) 
C  ' 

C     CUtSff?  ^^  ™^  ^°"^^"^  ACCESSED  FOR  A  SLOPING  INCREASE  WITH 
C 

C0MM0N/TW0/N,S(12),yNP(12),S0LIN(12),Z(12) 
COMMON/THREE/ YT(  12) 

DIMENSION  XL(12),YL(12) 
DO  1=1, N 

YL(I)»LOG(S(I)#S(I)*yNP(I)) 
XL(I)=S(I)        ^  .,  . 

END  DO  "■  " 

CALL  LSANAL(N,XL,YL,RINCEP,SLOPE,CO) 

C=EXP(RINCEP) 

THEOK=— SLOPE 

WRITEdClOO)  EXK  '- 

100  FORMAT (//•  '.SX. 'EXPERIMENT  K  VALUE  =  •  El  6  9) 
WRITEdO.IIO)  .c-io-i^; 

110  FORMATC/'  '.SX.'THE  EQUATION  IS  OF  THE  FORM') 

WRITE(10,120)        v-l.^-.  ^  ^""  ' 

120  FORMATC  '.ISX.'NCZ)  -  o'   IF  0<Z<K' ) 

WRITE(10,130) 
130  FORMATC  ',22X,'C(Z-K)   IF  Z>K' )   -*-,  .* .^- 

WRITEdO.UO)  C  ^   ^  • 

140  FORMATC/'  '.SX.'C  =  ',E16.9)       -     .   . 

WRITE(10,150)  THEOK 
150  FORMATC  '.SX.'K  =  •,E16.9) 

PERCENT=100.0»ABS(THEOK-EXK)/EXK 
WRITE(10,160)  PERCENT 

^^°  !:S'^fyy:;'5x.'J  error  between- exptl.  and  THEOR.  K  =  '. 

otFo.2,  '>  • )  • 

DO  1=1,  n  ..,  ;=  ■•  • 

YT(I)=C»EXP(-THE0K*S(I))/S(I)#»2 
IF  (Z(I).LE.THEOK)  THEN 

SOLIN(I)=0.0 
ELSE 

SOLIN(I)=C«(Z(I)-THEOK) 
END  IF 

END  DO  .   ■  ,:         >: 

RETURN  '■<:.' 

END  /  ■  '  ;  ^'  ■ 


C 

c 
c 
c 
c 

c 
c 
c 


SUBROUTINE  POLCUT(EXK )  -y'-fi:,  - ' 

THIS  IS  THE  ROUTINE  ACCESSED  FOR  A  POLYNOMIAL  FIT  WITH  CUTOFF. 


5^.*._:   V,r  -1C    ■  .,    :••*;  f     ■■   V'       •  ^ 


C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12),Z(12) 

C0MM0N/THREE/YT(12) 

DIMENSION  XL(12),YL(12) 

DO   1-1, N 

YL(I)-LOG(EXP(-EXK*S(I))/YNP(I))      "'  •'  /  ^        -~ 
XL(I)=LOG(S(I)) 

END  DO  .  .  : 

CALL  LSANAL(N,XL,YL,RINCEP, SLOPE, CO)        .. 

C-EXP(-RINCEP) 

WRITEdO.lOO)    EXK 
100   FORMAT(//'    ',5X,' EXPERIMENTAL  K  VALUE   =    ',E16.9) 

WRITE(10,110) 
110   FORMAT (/'    •,5X,'THE  EQUATION   IS   OF  THE   FORM') 

WRITE(10,120) 
120   FORMATC    M5X,'N(X)    =  0        IF  0<Z<K' ) 

WRITE(10,130) 
130   FORMATC    ',22X,'C(Z-K)»»(MU-1)    IF   Z>K' ) 

WRITE(10,135) 
135   FORMATC    •,22X,' •) 

WRITE(10,137) 
137   FORMATC    •,2ilX,  •  GAMMA (MU)')      • 

WRITE(10,140)  C  ,;:'-     "-, 

140  FORMATC/'    ',5X,'C=    ',E16.9) 

WRITE(10,150)   SLOPE 
150   FORMATC    ',5X,'MU  =■    ',£16. 9)  "■    " 

GMU=GAMMA(  SLOPE)  '   '   ■--;:''.; 

WRITE (10, 160)  GMU  '  '':   -r    -  • 

1  60  FORMAT C  ' ,  5X , '  GAMMA  (MU )  =  ' ,  El  6 .  9 )   •  •   '■  ■      ''■  '  • 

DO  I=>1  ,N  -"  ':   v 

YT(I)=C»EXP(-EXK»S(I))/S(I)»»SLOPE 
IF  (Z(I).LE.EXK)  THEN 

SOLIN(I)=.0.0 
ELSE 

S0LIN(I)-C»(Z(I)-EXK)»*(SL0PE-1.0)/GMU 
END  IF 
END  DO 

RETURN  ,  ;^'^  :.;       ..  ,  ,  '  .,: 

END  ■  ~'-    . 


C 
C 

c 
c 
c 

c 


SUBROUTINE  OVERLAY 


C     THIS  IS  THE  ROUTINE  ACCESSED  FOR  AN  OVERLAYER  DISTRIBUTION. 
C 

C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12),Z(12) 
COMMON/THREE/ YT(1  2) 
DIMENSION  XL( 12), YL(1 2) 
REAL  K 
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D3=1.0E10 
CI =0.0 

DO  STEP=0.001,3.0,0.001  ^  \-      "  T    -: 

DO  1=1,  N  ■■■     ^/rvii 

,,    -;     ■  YL(I)=S(I)*YNP(I)  ■     '    •-        ■ 

XL(I)=1-EXP(-STEP*S(I)) 

END  DO  ;  I  r:rK^T^ 

CALL  LSANALCN, XL, YL.RINCEP, SLOPE, CO) 
IF   (ABS(CO).GT.CI)   THEN 

CI -CO  ;  ''-..  '.\:v  r\  -'--fy    ff-  ^-  '    - 

C2=SL0PE       ■ .  i  •:  r-  r  '~  P  i  ■  ''-'^'^ '  ' ^    - 

C3=RINCEP   --^  ^■*'  '^  --  »■'  -^'  ■  -i-Jv   -     J   - 
■   :^-     C4=STEP  .  v'  '■ 

END  IF 

IF  (ABS(RINCEP).LT.D3)  THEN 
DI'CO 

^   .;-;.■        D2=SL0PE  ■,-,;■; 

D3=RINCEP  •;."'.' 
D4=STEP 
END   IF  J^;.   ; 

END  DO  -'-.^^Z 
WRITEdO.IIO)  i,     • 

110   FORMAT(//'    ',5X,'THE  EQUATION   IS  OF  THE   FORM')      ^■'     ■       '    ' 

WRITE(10,120)  '^       :--,.. 

120   FORMATC    M5X,'N(Z)  -  C      IF  0<Z<K' )  :    /-^^^ 

WRITE(10,130)  ..-;:,.,-   , 

130  FORMATC  ',22X,'0   IF  Z>K' )  ^'   A'  -   v 

WRITE(10,133) 

'^^  WRmOO.'lilofc^"'  ''"'  "'  '''°''''°  '°  ™^  CORR.COEF.  GIVES') 
140  FORMATC/'  '.SX.'C  -  ',E16.9) 

WRITE(10,150)  C4 
150  FORMATC  '  ,5X, 'K  =  ,E1 6. 9)  •:■•■-     ^ 

WRITE(10,160)  CI  "  ••      ,; 

150  FORMATC  '.SX.'CORR.  COEF.  =  '  F9  6)  ■ 

WRITE(10,170)  C3  '  ^ 

'''  WRmco.iso)'"""''''^  ^'"°"''  ''  '''°^  =  ''^^^-^^        ■   ' 

'^°  WRmOolliofD™  ^^ST^^F"  ACCORDING  TO  THE   INTERCEPT  GIVES') 

WRITE  do,' 150)   D4  '  ■  "    ' 

WRITE(10,160)  D1  ■-.  ■ 

WRITE(10,170)  D3 

C=0.5»(C2+D2)  '  'r''.:' 

K=0.5*(C4+D4) 

DO  1=1, N 

YT ( I ) -C» ( 1 -EXP (-K»S (I ) ) ) /S ( I )  - V 

IF   (Z(I).GT.K)   THEN 
SOLIN(I)=0.0 

ELSE 

SOLIN(I)=C  "' 
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END  IF 
END  DO 
RETURN 
END 


C 
C 

c 
c 

c 

c 
c 
c 


SUBROUTINE  POWLAW  ,  ;    '   - 

THIS  IS  THE  ROUTINE  ACCESSED  FOR  A  POWER  LAW  DISTRIBUTION. 

C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12),Z(12)  •-.      ' 

C0MM0N/THREE/YT(12)  ,        ,v  -^- 

DIMENSION  XL(12),YL(12) 
DO   I=.1,N 

yL(I)-LOG(YNP(I))     .j    .i/  ^^?><.' 

xL(i)»LOG(s(i))       ;in''     '^'''' 

END  DO  ■'    ^'  ■~'     "■^' 


.^ 


CALL  LSANAL(N,XL,YL.RINCEF,SLOPE,CO)      '  '■  -' 

THEOK=-SLOPE  .    .,     ; 

GAM-GAMMA  (THEOK)  "  '     '' 

C=EXP(RINCEP)/GAM  '  '• 

WRITEdO.IIO)  -.    ■ 

110   FORMATC//'    '.SX.'THE   EQUATION   IS   OF  THE   FORM') 
WRITE(10,120) 

120   FORMATC    M5X,'N(Z)    =  CZ»»(K-1)')  > 

WRITE(10,140)  C 

140  FORMATC/'  '.SX.'C  -  ',E16.9) 

WRITE(1G,150)  THEOK 
150  FORMATC  ',5X,'K  =  ',£16. 9) 

DO  1-1, N 

YT ( I ) =C»GAM»S ( I ) «» ( -THEOK ) 

S0LIN(I)-C»Z(I)»»(THE0K-1.0) 
END  DO 

RETURN 

END  ... 

C  •■     ■■  ,  -•^: ■■:...•->■"■  ■:   ,     ■  -        '        '''^   " 

C     ..  .  .■  ''■     ■■-v-..:'^,. -"-  ;     -  '    '  V  ■  .  . 

c       -..■  ^■■•. -.•'■./':  :.-.■.::;.  ■.,-,^ 

C       ,.■•■■'  -   -.-■•  ,;.■  ', .  ■■  '-     )  ■'■.:'.         - 

c  '-';■■'■; 

SUBROUTINE  POWCUT(EXZ)^" 
C 

I  WITh"cUTOFF™^  "°""'*^  ACCESSED  FOR  A  POWER  LAW  DISTRIBUTION. 

C  *        .■.■■•'•■•■■•;  }.     . 


C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12),Z(12) 
C0MM0N/THREE/YT(12) 

DIMENSION  XL(12),YL(12) 
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DO  1=1, N       /  v*-^iM;''  '^  -  " 

YL(I)=LOG(YNP(I))     ^~'  '■  v.-  ' 

xL(i)=LOG(s(i)-Exz)  '  •  ,. :  r:  ./A  . 

END  DO  •  ;  'v,;^'S(; 

CALL  LSANAL(N,XL,YL,RINCEP,SLOPE,CO) 

THEOK=-SLOPE 

GAM=GAMMA(THEOK)        '     ;>;^ 

C-EXP(RINCEP)/GAM      .  -     '•  ' :.  . 

WRITEdO.IIO) 

110   FORMAT(//>     ',5X,'THE  EQUATION   IS  OF  THE  FORM') 
WRITE(10,120) 

120  FORMATC  M5X,'N(Z)  =  CZ»»(K-1  )»EXP(AZ) ' ) 

WRITE(10,140)  C 
140  FORMATC/'  ',5X,'C=  ',E16.9) 

WRITE(10,150)  THEOK  .  ^-: 

150  FORMATC  '.5X,'K  =  ',£16. 9)   i  ' 

WRITE(10,150)  EXZ 
160  FORMATC  ',5X,'A  =  ',E16.9)   .'/  • 

DO  1=1,  N  ■'-•■- 

YT(I)=C»GAM»(S(I)-EXZ)**-THEOK 

SOLIN(I)-C*EXP(EXZ»Z(I))»Z(I)»*(THE0K-1.0) 
END  DO 

RETURN 

END 


C 
C 

c 
c 
c 

c 
c 
c 


SUBROUTINE  EXPONE  :  :' - 

THIS  IS  THE  ROUTINE  ACCESSED  FOR  AN  EXPONENTIAL  DISTRIBUTION. 

C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12),Z(12) 
C0MM0N/THREE/YT(12) 

DIMENSION  XL(12),YL(12)        4  J  Q i 

DO  1=1,  N  r.    'v-;  :  •. 

YL(I)=YNP(I)»S(I) 
XL(I)-YNP(I) 
END  DO 

CALL  LSANAL(N,XL,YL,RINCEP,SLOPE,CO) 
WRITE(10,110) 

110  FORMAT(//'  •,5X,'THE  EQUATION  IS  OF  THE  FORM') 

WRITE(10,120)  . 

120   FORMATC    ',15X,'N(Z)    -  C«EXP(AZ)') 

WRITEdO.UO)   RINCEP  ■- 

140  FORMATC/'    '.SX.'C"    ',£16. 9)         ..^  -   ^-Z 

WRITEC10,150)   SLOPE 
150   FORMATC    ',5X,'A  =    ',E16.9) 

DO   1=1, N 
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YT(I)=RINCEP/(S(I)-SLOPE) 
SOLIN ( I ) =RINCEP*EXP (SLOPE»Z ( I ) ) 

END  DO 

RETURN 

END  ,! 


C 
C 

c 
c 
c 

c 


SUBROUTINE  LININT 


r  x,o»™i^  ^^  ™^  ROUTINE  ACCESSED  FOR  A  LINEAR  FUNCTION  WITH  A 

C     NON-NEGATIVE  INTERCEPT.  wxin  h 


C 


C 

c 
c 
c 
c 

c 
c 
c 


C0MM0N/TW0/N,S(12),YNP(12),S0LIN(12).Z(12) 
C0MM0N/THREE/YT(12) 
DIMENSION  XL(12),YL(12) 
DO  I-I.N 

YL(I)=S(I)»S(I)*YNP(I) 
XL(I)=S(I)       .  . 

END  DO  ■  '  ■'■■  - 

CALL  LSANAL(N,XL,YL,RINCEP,SLOPE,CO) 
WRITE(10,no) 

110   FORMAT (//•     '.SX.'THE  EQUATION   IS   OF  THE   FORM') 

WRITE(10,120) 
120   FORMATC    M5X,'N(Z)    =  A»Z   +  C,    OO  • ) 

WRITEdO.UO)   RINCEP 
140  FORMAT(/'    '.SX.'A  -    ',E16.9) 

WRITE(10,150)   SLOPE 
150   FORMATC    ',5X,'C  =    ',£16. 9) 

WRITE(10,160)   -SLOPE/RINCEP    ■     ^     *-^ 
160   FORMATC    ' ,  5X, 'Z-INTERCEPT  =    ',£16. 9) 

DO  1=1, N 

YT ( I ) =SLOPE/S ( I ) +RINCEP/S (I ) **2 

SOLIN(I)=RINCEP»Z(I)+SLOPE 
END  DO 

RETURN  :;1 

END  t  ■'■;'•■ 


SUBROUTINE  LSANAL(N,X.Y,R,S,CO ) 

OF  ?HE^NPOT'nATr"'?Sf .'  ''^'''  ''^'^'  ^^AST  SQ"A«ES  ANALYSIS 
rnJ^T  .rTn.^^.j;iJ'  .™^  ^^°P^  AN°  Y-INTERCEPT  ARE  OUTPUT.   THE 


C     CORRELATION  COEFFICIENT  IS  ALSO  AVAILABLE  FOR  OUTPUT. 


C 
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C 
C 

c 
c 
c 

c 


DIMENSION  X(12),Y(12) 
A=0.0 
B-0.0 
P=0.0 
Q=0.0 
E=0.0 
DO  1=1, N 
A=A+X(I) 
B=B+Y(I) 
P=P+X(I)»X(I) 
Q-Q+y(I)*Y(I) 
E=E+X(I)»Y(I) 
END  DO 
SX=P-A»A/N 
SY=Q-B*B/N 
XY=E-A»B/N 
S=-XY/SX 
XBAR=A/N 
YBAR=B/N 
R=»YBAR-S»XBAR 

CO=XY/SQRT(SX»SY) 

RETURN 

END 


FUNCTION  GAMMA (ARG) 


C        THIS  ROUTINE  CALCULATES  THE  GAMMA  VALUE  FOR  A  GIVEN  POSITIVE 

C  REAL  NUMBER  INPUT.   IT  MAKES  USE  OF  A  FORMULA  TAKEN  FROM 

C  DWIGHT'S  "TABLES  OF  INTEGRALS  AND  OTHER  MATHEMATICAL  DATA". 

c 


LOGICAL  FLAG 

INTEGER  WG  :*^-       . 

REAL»8  SP(500) 

DATA  EU, ERR    /O. 5772157, 1.0E-10/ 

OPEN   (UNIT=16,FILE='SP.DAT',STATUS='0LD') 
FLAG-. TRUE. 
SP(1)-0.0 
POM- 1.0 

WG=2  ::■■■■ 

ISI=-1  '";;,•  :^;:>:", 

DO  1=2,26,2  '  :^",;...'W.:r>  :■'  ;•- 

READ(16,100)   SP(1),SP(I+1) 

100  F0RMAT(2(5X,F10.8))  

END  DO  '*^".i"''-''  '.'     '  ■"'■ 

DO  1=28,500  '-"  '-'  C^ 

SP(I)=1.0  ^'         ,    S      ■ 

END  DO 
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BR=ARG-1.0 

IF    (BR.GE.1.0)   THEN 
TX=INT(BR) 

BR=BR-TX  .>  ^%  -••-»:   ^    - 

FLAG=. FALSE.  -M ^    ,  U 

END   IF  .w  .,.   Ut^\l<>^\.Xi;^ 

GP=-EU»BR 

DO  WHILE    (ABS(P0M).GT.ERR.AND.WG.LE.500) 

P0M=(-1»ISI)*SP(WG)»(BR*»WG)/FL0AT(WG) 
GP=GP+POM 

WG=WG+1 

ISI=-1»ISI 
END  DO 

GAMMA-EXP(GP)  ,,  .  v  V 

IF    (FLAG)   RETURN      '•   •      :  7     '  .' 

BR=BR+1.0 

DO  WHILE    (ABS(BR-ARG ).GT. ERR) 

GAMMA=BR»GAMMA 

BR=BR+1.0 
END  DO 
RETURN 
END 
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